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Chapitre 1 



Introduction 



Depuis de nombreuses annees, dans differents pays, les informaticiens ayant 
quelques pretentions academiques ont lutte pour etablir leur discipline de maniere 
independante. Sans vouloir dire que la lutte est terminee (certains n' ayant pas encore 
accepte que la terre n'est pas plate), on peut constater que, dans les universites 
respectees, il existe des laboratoires d'informatique independants, des diplomes 
Specialises, et -les enseignants et/0U chercheurs en informatique sont desormais 
considered comme des scientifiques a part entiere. 

Pourquoi cette lutte ? Et pourquoi en parler dans un livre sur Falgorithmique ? 
Le fait est que les informaticiens ont represente — et represented toujours — un 
enjeu economique. Comme cet enjeu a ete concretise par des moyens materiels et 
financiers mis a la disposition des enseignants et des chercheurs en informatique, 
tout un chacun a eprouve le besoin de reclamer l'6tiquette. Le tri entre les "vrais" et 
"faux" informaticiens n'est pas termine. D'ailleurs, a notre avis il ne le sera jamais, 
et c'est peut-etre heureux ainsi. 

Malheureusement, en voulant affirmer leur independance par rapport aux autres 
disciplines, les informaticiens ont perdu une partie de l'essentiel. En se concentrant 
sur les techniques non-numeriques (importantes et indispensables), ils ont perdu 
jusqu'a la notion de l'existence des nombreL reels. De meme, un peu en singeant les 
mathematiciens, qui ont montre la voie vers la tour d'ivoire, le besoin scientifique 
mais aussi psychologique de la construction dune (ou de plusieurs) theorie(s) a fait 
perdre la vraie justification de notre existence : l'ecriture de programmes qui sont 
utiles. On est done en presence de plusieurs guerres, numerique contre non- 
numdrique, theorie contre application, utilisateurs contre specialistes, vendeurs de 
vent contre professionnels serieux. Si certaines guerres peuvent etre utiles et 
salutaires, d'autres resteront toujours steriles. 
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Ce livre ne saurait bien sur en corriger les ecarts. Mais nous voulons 
temoigner de notre foi dans l'existence de l'informatique en tant que discipline 
independante, mais surtout utile. L'informaticien comme le mathematicien — meme 
si ce dernier l'a peut-etre oublie — est un esclave des autres. Sa raison d'etre est de 
rendre service, c'est-a-dire resoudre des problbmes d' application. II n'y a pas 
d'informatique academique differente de celle de l'application numerique ou de 
gestion. II n'y a pas une micro-informatique differente de l'informatique "classique". 
II y a une seule discipline, appelee a intervenir dans un tres grand nombre de 
domaines de l'activite humaine. 

Dans cette optique, la formation des informaticiens dans les universites et les 
ecoles d'ingenieurs doit se faire de maniere equilibree. Un de nos criteres est que 
si nous ne traitions pas tel sujet, nous pourrions par la suite avoir honte de nos 
Aleves ?". Le monde du travail s'attend a ce que nos eleves sachent programmer (dans 
le langage qu'utilise la compagnie concernee), qu'ils connaissent un minimum de 
methodes de resolution de problemes numeriques, et que "probabilites" et 
"statistiques" ne soient pas des mots 6sot6riques a rechercher dans le dictionnaire. Le 
cours d6crit dans ce livre essaie de prendre en compte ces considerations. Nous 
l'enseignons, dans des variantes appropriees, depuis vingt-cinq ans, a des 61eves de 
premier et de deuxieme cycle, en formation permanente, dans des diplomes 
specialises ou non. L'experience a montre que l'enseignement de ce cours de base de 
l'informatique est difficile, que personne ne possbde une recette miracle, mais que 
tout le monde (surtout ceux qui n'ont jamais ecrit un logiciel utilise par d'autres 
personnes) pense l'enseigner mieux que les autres. 

Nous presentons done ici, humblement (ce n'est pas notre tendance profonde), 
quelques idees d'un cours d'algorithmique et de programmation dont le niveau 
correspond a la premiere annee d'un deuxibme cycle pour specialistes en 
informatique, en supposant que les eleves concernes n'ont pas necessairement subi 
une formation prealable dans la matiere, mais qu'ils sachent tous un peu 
programmer Cela pose d'ailleurs un probleme particulier. Nous avons l'habitude de 
melanger des etudiants d'origines divers. Les uns peuvent avoir obtenu de bons 
l"6sultats dans un IUT d'informatique, tandis que d'autres n'ont que peu touche un 
clavier (situation de plus en plus rare, mais l'&riture ^ m programme de 10 lignes en 
BASIC peut etre consideree comme une experience vide, voire negative, en 
informatique). A la fin de la premiere annee de deuxieme cycle, il reste une 
correlation entre les etudes prealables et les resultats academiques. Cette correlation 
disparait au cours de la deuxieme annee, provoquant quelques remontees 
spectaculaires au classement. La seule solution que nous avons trouve a ce probleme 
de non-homogeneite est de l'ignorer en ce qui concerne le cours, mais d'amenager des 
binomes "mixtes" en TP. L'experience jusqu'alors est positive. 
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C'est en respectant l'idee que tous nos 616ves savent un minimum sur la 
programmation qu'a disparu le chapitre de ce livre destines aux debutants. Ainsi, les 
types de donnees de base (entiers, rdels, caracteres), leurs representations en machine 
et les instructions simples (affectations, boucles, conditions) ne sont pas definis dans 
ce volume. 

Dans sa forme actuelle, le cours dure une annee scolaire, au rythme de deux 
heures par semaine. II est necessairement accompagne de travaux diriges et de travaux 
pratiques. Chez nous, il y a trois heures de travaux diriges et quatre heures de travaux 
pratiques par semaine. De plus, la salle informatique est ouverte en libre service aux 
Ctudiants autant de temps que possible en respectant les regies el6mentaires de 
Security. Ce n'est qu'au prix de la mise a disposition d'un materiel suffisant que les 
6tudiantS peuvent reellement apprendre. 

Nous proposons de nombreux exercices et problbmes. Le sujet que nous 
attaquons necessite un investissement personnel important. Le cours doit servir a 
stimuler des efforts individuels et la realisation de projets de toutes sortes. II doit 
obligatoirement se terminer par la creation d'un logiciel en grandeur nature, de 
preference produit par un petit groupe d'eleves (deux a quatre). 

Les exemples dans ce livre sont r6dig6s dans un langage de programmation 
fictif qui ressemble a PASCAL. Comme boutade bien connue des quelques 
universites l'ayant subi, le langage dans lequel nous redigeons nos algorithmes a pris 
le nom "GRIFFGOL", qui resulte de reflexions de l'epoque de MEFTA [Cunin 1978]. 
C'est un style de programmation relativement independant d'un langage de 
programmation particulier, c'est-a-dire que nous utilisons les concepts fondamentaux 
dans la forme qui nous semble la plus appropriee. La construction des algorithmes se 
passe mieux quand on se permet une certaine liberte d'expression. Leur mise en 
ceuvre dans un ordinateur necessite une simple mise en forme en fonction du langage 
et du compilateur disponibles. 

Un corollaire est que nous disons a nos etudiants qu'ils doivent toujours 
r^pondre "oui" a la question "connaissez-vous le langage X ?", quelle que soit la 
valeur de X. En effet, sous condition que le langage soit de style algorithmique 
classique, l'apprentissage d'un langage et/ou d'un compilateur inconnu ne devrait 
durer que quelques jours. D'ailleurs, un exercice classique que nous pratiquons est de 
faire recoder un travail pratique dans l'un ou l'autre des langages a grande diffusion 
que nos etudiants ne connaissent pas. C'est ainsi qu'ils absorbent, par exemple, 
FORTRAN. En pratique, a l'heure actuelle, ils programmed dans une version de 
PASCAL, avec des prolongements en C. Le choix est remis en cause a chaque 
rentree universitaire, car nous ne sommes pas des missionnaires d'un langage 
quelconque. 
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Algortthmique et programmation 
1.1. Quelques mots sur I'environnement 

Une petite phrase ci-dessus merite qu'on s'y attarde un peu plus longtemps. 
Nous avons parle de la disponibilite de materiels, essentielle & l'apprentissage de la 
programmation, qui est une activite constructive. On ne peut apprendre qu'en 
s'exercant. Or pour s'exercer de maniere satisfaisante, l'ideal, si nous pouvons nous 
permettre le parallele, est que les ordinateurs doivent etre comme les WC : il y en a 
toujours un de libre quand on en a besoin, voire meme quand on en a envie. (Cette 
phrase a ete prononcee pour la premiere fois, a notre connaissance, par P.M. 
Woodward, du Royal Radar Establishment a Malvern, Angleterre. C'est un plaisir de 
rendre cet amical hommage a un maitre mal connu, en se rappelant qu'en 1964 il 
fallait une vision ties large pour s'exprimer ainsi.). 

Certains de nos collegues restreignent volontairement le cote experimental de 
la programmation, dans le but d'imposer, des le debut de l'apprentissage, toute la 
rigueur necessaire. C'est une reaction saine par rapport a une situation historique 
datant de l'epoque ou la programmation se faisait n'importe comment. Mais nous 
n'allons pas jusqu'a empecher nos 616ves de faire leurs betises. C'est en comparant 
des versions sauvages de programmes avec d'autres qui sont bien faites qu'ils 
comprennent reellement l'interet d'une methodologie. 

Ainsi, nous voulons que les eleves passent beaucoup de temps dans la salle des 
machines. Au debut, ils travaillent mal, mais deviennent — pour la plupart — 
raisonnables a la fin de la premiere annee. Cette approche profite d'une certaine 
fascination pour l'ordinateur (la jouissance de le dominer ?), qui s'estompe apres cette 
premiere annee. Le fait de vouloir reflechir, plutot que de se lancer immediatement 
sur la machine, est un signe majeur de maturite chez l'eleve. Cette etape ne peut etre 
franchie que parce que le materiel est toujours disponible. Un etudiant ne doit pas 
etre stresse par la longueur d'une file d'attente, ni frustre par des difficultes 
materielles. Nous notons d'ailleurs que, bien que les programmes ecrits en deuxieme 
annee soient assez importants, l'occupation des postes de travail diminue. 

L'entrainement a la programmation est une necessite pour tous les 
informaticiens, quelle que soit leur experience. Une source de stimulation pour les 
eleves est de travailler en contact avec des enseignants qui programment bien. Tant 
qu'ils sentent qu'il leur reste du chemin a faire pour arriver au niveau de rendement de 
ce modele, ils se piquent au jeu. Cela signifie que l'enseignant doit lui-meme 
continuer a programmer regulierement, comme le musicien qui continue a faire des 
gammes. Meme si nous n' avons plus le temps de produire de gros logiciels, il faut 
s'astreindre a resoudre regulierement des petits problemes. Ceux-ci peuvent, par la 
suite, contribuer au renouvellement de notre stock d'exemples et d'exercices. 
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1.2. Notes bibliographiques sommaires 

II existe deja, en francais, un certain nombre de livres [Arsac 1980], [Arsac 
19831, [Berlioux 19833, [Boussard 19831, [Courtin 1987a,b], [Gerbier 19771, 
[Gregoire 1986, 1988], [Lucas 1983a,b], [Meyer 19781, [Scholl 19841 dans le 
domaine que nous abordons ici. Un pourcentage eleve porte un nom d'auteur qui 
nous est familier, car il s'agit de collegues et sou vent amis. Cela indique que les 
reflexions concernant ce sujet sortent, en France au moins, d'un nombre limite 
d'ecoles qui, en plus, ont pris l'habitude de discuter. Le nombre de livres indique 
1' importance que chacun accorde au sujet, et les differences d'approche demontrent 
qu'il est loin d'etre epuise. Comme dans la tradition litteraire, il y aura toujours des 
idees differentes sur ce sujet. En continuant le parallble avec l'ecriture, nous 
recommandons aux eludiantS de prendre connaissance de plusieurs styles differents, 
puis de developper leur propre style en profitant des apports de chacun. 

Ayant indique quelques livres en francais, il serait deraisonnable de laisser 
l'impression que la France possede un monopole des idees. Au niveau international, 
il existe une bibliographic consequente en langue anglaise, dont voici les r6f6rences 
qui correspondent a une s61ection personnelle parmi les ouvrages les plus connus : 
[Aho 1974, 1983], [Dijkstra 1976], [Gries 1981], [Ledgard 1975], [Wirth 1976, 
1977 1. 



1.3. R emerciements 

Sur un sujet tel que le notre, il serait impensable de citer tous ceux qui ont 
influence notre reflexion. Nous nous limitons done a la mention de quelques groupes 
de personnes, en nous excusant vis-a-vis de tous les collegues que nous ne citons pas 
individuellement. 

Comme premiere influence directe, il y a eu le groupe de P.M. Woodwarda 
Malvern au debut des annees soixante. L'auteur y a fait ses premieres armes, et ses 
premiers cours, & partir de 1962, sous la direction de J.M. Foster et D.P. Jenkins, 
avec I.F. Currie, A.J. Fox et P.R. Wetherall comme compagnons de travail. Le 
foisonnement d'idees a Grenoble entre 1967 et 1975 a 6t6 d'une grande importance. 
Signalons seulement une collaboration directe avec P.Y. Cunin, P.C. Scholl et J. 
Voiron [Cunin 1978, 1980], bien que la liste aurait pu etre nettement plus longue. 
L'6cole de C. Pair a Nancy a montre la rigueur et a donne des opportunites de 
comparaison de styles. Finalement, la creation du diplome d'ingenierie informatique 
a Marseille en 1985 a provoque un effort de synthese dont le resultat est ce volume. 
De nouvelles versions du polycopie ont vu le jour a Nantes en 1990 et 1991. 
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Des rencontres ponctuelles ont aussi eu leur importance, en particulier a travers 
les &oleS organises par F.L.Bauer [Bauer 1974, 1975, 1976, 1979]. Celles-ci ont 
permis de travailler avec son equipe et de rencontrer et de se confronter avec 
E.W.Dijkstra, G.GOOS, D.Gries, J.J.Horning, P.C.Poole et W.M.Waite, parmi de 
nombreux autres collegues. 

Le feedback de generations d'etudiants nous a permis de constater que 
l'apprentissage de l'informatique n'est pas toujours aise si Ton veut atteindre un bon 
niveau. Nous remercions ces jeunes (certains ne le sont plus) pour leur apport, dont 
ils ne se sont pas toujours rendu compte (nous aussi, nous apprenons !). Emmanuel 
Gaston du DU-II et un groupe du mastere d' intelligence artificielle de Marseille, 
promotion 198889, ont corrige un certain nombre d'erreurs de francais. Nathalie 
Wurbel a aid6 en signalant des eneurs de francais et de programmation. Christian 
Paul et Annie Taftier out egalement apporte des conections au niveau du francais. Au 
cours d'un projet du DU-II, Vita Maria Giangrasso et Thieny Guzzi [Giangrasso 
1989] ont mene a bien une etude comparative de certains algorithmes. Le journal 
Jeux et Strategic, source de problemes interessants, nous a aimablement permis d'en 
utiliser dans certains exemples. 

Ce livre a ete prepare sur des micro-ordinateurs mis a notre disposition par 
diff&entS organismes, dont feu le Centre mondial d'informatique et ressources 
humaines, le CNRS, l'universite d'Aix-Marseille III, l'institut mediterraneen de 
technologie et l'universite de Nantes, que nous remercions. 

Un dernier mot sera pour mon collegue, et surtout ami, Jacek Gilewicz, 
professeur a l'universite de Toulon. Au debut, nous avions envie de preparer un livre 
combinant l'algorithmique numerique et non numerique. Pour differentes raisons, ce 
projet n'a pas vu le jour. L'utilisation du "nous" dans cette introduction represente ce 
pluriel. Je l'ai laisse dans cette version definitive, en esperant que Jacek r6digera le 
volume numerique dont nous avons besoin, en meme temps que je lui temoigne ma 
reconnaissance pour son soutien et son amitie sans faille. Ce livre lui est dedie. 

1.4. Le choix de programmes 

On apprend a ecrire des programmes en pratiquant. C'est pour cette raison que 
nous travaillons a partir d'exemples. Ceux-ci sont de plusieurs types : 

■ les classiques, qui se trouvent deja dans d'autres ouvrages, mais qui sont 
essentiels a la culture d'un informaticien, 

■ les pedagogiques, que nous avons crees ou repris comme materiel de base. Ici, 
on aurait pu en choisir d'autres, mais chaque enseignant a sa liste d'exercices, 
souvent partagee avec des collegues, 

• les amusements, qui sont la parce qu'ils nous ont fait plaisir, mais qui 
presentent neanmoins un interet pour l'6tudiant 
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Libre a chacun d'apprecier notre choix. De toute facon, il est entendu que chaque 
enseignant mettra la matiere "a sa sauce". 

Dans la mesure du possible (ou en fonction de nos connaissances), nous avons 
essaye d'attribuer la paternite des exemples, mais beaucoup d'entre eux ont des 
origines deja inconnues. Toute information suppldmentaire sera le bienvenue. 
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Chapitre 2 

Des programmes pour commencer 



2.1. Le mode d'un vecteur 

Ce probleme nous est venu de D. Gries, qui l'utilise depuis longtemps dans ses 
cours d' introduction k l'informatique. U a 6t6 repris par differents auteurs dans le cadre 
de l'analyse d'algorithmes [Arsac 1984], [Griffiths 19761. 

On considere un vecteur, disons d'entiers, dont les elements sont ordonnes. Son 
mode est la valeur de l'element qui y figure le plus grand nombre de fois. Ainsi, 
prenons le vecteur suivant : 

(1 3 3 6 7 7 7 11 13 13) 

Le mode de ce vecteur est la valeur 7, qui figure trois fois. Nous allons ecrire un 
progmmme qui prend en entree un vecteur ordonne et qui livre le mode du vecteur en 
sortie. 

2.1.1. Construction de la premiere version du programme 

Commencons par une premiere idee d'algorithme. On note que 
l'ordonnancement du vecteur fait que les differentes occurrences d'un element qui se 
r6pe-te sont contigues. II s'agit done de trouver la plus longue chaine d'elements 
identiques. On va considerer les elements & tour de role, en se rappelant de la plus 
longue chaine vue jusqu'a present. Pour chaque nouvel element, on se posera la 
question de savoir si sa lecture mene h une chaine qui est plus longue que la 
prficedente. Dans ce cas, e'est la chaine courante qui devient la'chaine la plus longue. 



ALGORITHMIQUE et programmation 

Pour garder les informations necessaries, il faut disposer des valeurs suivantes : 

■ n est le nombre d' elements dans le vecteur v, 

. i est le nombre d'eTdmentS qui ont deja etc considered, 

■ lmax est la longueur de la plus longue chaine en v[l..i], 

. m est l'index en v d'un element dans la chaine de longueur lmax : 

v[m] = mode(v[l ..i]), 

. lc est la longueur de la chaine courante (a laquelle appartient v[i]). 
On appellera cet ensemble de definitions Vital du monde ("state of the world"). 

Le programme 2.1, qui met en ceuvre ces idees, est une application du schema 
suivant : 

Initialiser 

TANTQUE NON fini 
FAIRE avancer 
FAIT 

Le programme est commente par la suite. 



DEBUT DONNEES n: entier; 

v: TABLEAU [1 ..n] DE entier; 
VAR i, lmax, m, lc: entier; 
i:=1;lmax:=1;m:=1;lc:=1; 
TANTQUE kn 
FAIRE i:=i+l ; 

SI v[i] = v[i-1 ] 

ALORS lc:=lc+1 ; 

si lolmax 

ALORS lmax:=lc; m:=i 

FINSI 
SINON lc:= 1 
FINSI 

FAIT 

FIN 

Programme 2.1. Le mode d'un vecteur 



Les donnees du probleme sont n, la longueur du vecteur, et le vecteur ordonne, 
v- Les deux sont initialises par ailleurs. Apres les declarations des objets de l'etat du 
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monde (i, lmax, m, lc), la phase d'initialisation sert a leur dormer des valeurs 
permettant de demarrer. Cette ligne correspond au traitement du premier element : la 
chaine maximum est la chaine courante, de longueur 1. Dans le TANTQUE, comme 
i est le nombre d'elements deja traites, le test de la fin est bien kn (autrement dit, il 
reste des 616ments a considerer). Pour avancer, on prend le prochain 61ement (i:=i+l), 
en se demandant si l'examen de cet element met en cause l'etat du monde. II y a deux 
cas : soit on allonge la chaine courante (v[i] = v[i-l]), soit on commence une 
nouvelle chaine. La nouvelle chaine, decrite dans la partie SINON, est de longueur 1 
(lc:=l) et ne met pas en cause lmax ou m (qui ne sont done pas modifies). Si la 
chaine courante a ete allongee (partie ALORS), on augmente lc (lc:=lc+l) avant de 
tester si la chaine courante est devenue la chaine maximale (SI lolmax). Dans ce 
cas, hnax et m sont mis a jour (hnax:=lc; m:=i). 

Nous pretendons que cette construction demontre la justesse du programme 
donne. La demonstration depend de l'acceptation de la recurrence (a partir du cas i-1 
on cr6e le cas i, le cas 1 etant traite dans 1' initialisation) et de la validite de l'etat du 
monde. En particulier, dans la boucle, avancer est fait par i:=i+l et les autres 
instructions remettent les e^mentS de l'etat du monde a jour en fonction de ce que 
Ton y trouve. 

Pour etre complet, revenons sur la confirmation du nouvel etat du monde. 
i a avance de 1, ce qui est correct, lc a pris une valeur de 1 (nouvelle chaine) ou de 
lc+1 (allongement de la chaine en cours). La chaine courante devient la chaine la 
plus longue si lolmax. Dans ce cas, lmax et m recoivent des valeurs appropriees. 
Comme tous les elements de l'etat du monde ont des valeurs correctes, le programme 
entier est juste. 

2.1.2. Remarques methodologiques 

La construction d'un programme correct depend de la logique employee par son 
createur, Cette logique ne peut s'exercer que si les bases declaratives sont saines. 
C'est ici qu'une bonne formalisation de l'6tat du monde est fondamentale. Cet etat du 
monde peut etre decrit en francais (ou en anglais ou en polonais, nous ne sommes 
pas racistes) ou en style mathematique. Le choix entre les deux (ou les quatre) styles 
est une question de gout et de type d' application. 

Mais, que l'etat du monde soit ecrit en francais ou en formules mathematiques, 
il se doit d'etre precis. De nombreux programmeurs agrementent leurs programmes 
avec des commentaires dll type "i est l'index de l'element de v". Us font souvent 
l'erreur classique de confondre le dernier element traite avec son suivant (celui qui va 
etre mite). 
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L'utilisation de tableaux et de boucles impose d'autres contraintes. En 
particulier, il faut demontrer la terminaison de chaque boucle et confirmer que l'index 
de chaque reference a un element de tableau est entre les bornes (i doit etre contenu en 

[l..n]). 

Pour la terminaison de la boucle, il n'y a pas de probleme ici. La variable de 
controle i est augmented de 1 a chaque tour de la boucle, et arrivera done a n pour 
arreter le massacre. 

Les seules references a des 616mentS de tableaux sont dans le test SI v[i]=v[i-l]. 
i commence a un et la boucle s'arrete quand i=n. Avant l'affectation i:=i+l a 
l'int^rieur de la boucle, on a done I'indgalite' suivante : 

1 ^i<n 

Apres l'augmentation de i, ceci devient : 

1 <i£n 

On voit facilement que i et i-1 sont tous les deux a l'int6rieur du domaine l..n. 

Cette solution du problbme, d'une grande simplicity est correcte, mais elle 
n'est pas la meilleure possible. II existe une variante qui est moins longue, qui 
utilise moins de variables, et qui s'execute plus vite. II est meme surprenant que 
l'amelioration echappe a la quasi-totalite des programmeurs. Nous y reviendrons dans 
un paragraphe ulterieur (cf. 7.1). 

Par ailleurs, la solution donnee comporte une faiblesse. En general, il y a au 
moins un etudiant, ou une etudiante, dans la classe qui reste suffisamment r6veill6(e) 
pour poser la question suivante : 

"Qu'est ce qui se passe si deux chaines sont de meme longueur, e'est-a-dire si 
deux modes se pr6sentent ?", 

Bonne question ! La reponse peut etre de prendre n'importe lequel des modes 
possibles, mais la faiblesse reside dans le fait que notre specification du probleme ne 
parlait pas de cette possibilite. On voit qu'il n'est pas simple de bien specifier des 
algorithmes, surtout quand il faut etre complet. 

2.2. Recherche d'un objet 

Dans ce nouveau programme, il s'agit de trouver ['emplacement qu'occupe un 
objet (disons un entier) dans un tableau. C'est une operation menee frequemment 
dans des programmes de tous types. Supposons en premier lieu que nous ne savons 
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rien du tableau, e'est-a-dire qu'il faut en examiner tous les elements pour demontrer 
l'absence eventuelle de l'objet recherche. Bien sur, le programme s'arrStera des que 
l'objet est trouv6. Par la suite, nous montrerons comment rendre l'algorithme plus 
efficace en introduisant des connaissances concernant la forme des donnees. 



2.2.1. Recherche lineaire 

L'algorithme le plus simple examine a tour de role chaque element jusqu'a la 
r£ussite de la recherche ou l'epuisement des candidats. II necessite l'etat du monde 
suivant : 

DONNEES 

n: entier, le nombre d'elements de t, 

i: TABLEAU [1 ..n] DE entier, le tableau a examiner, 

objet: entier, l'objet a trouver, 

VARIABLES 

i: entier, le nombre d'elements deja examines, 
trouve: bool, trouve = t[i]=objet. 

Le programme 2.2 est une mise en ceuvre possible. 



DONNEES n, objet: entier; 

t: TABLEAU [1 ..n] DE entier; 
PRECOND n>0 

DEBUT VAR i: entier, trouve: bool; 
i:=0; trouv§:=FAUX; 
TANTQUE i<n ET NON trouve 
FAIRE i:=i+l ; trouve := t[i]=objet 
FAIT 

POSTCOND (trouve ET objet=t[i]) OU (NON trouve ET j=n) 
FIN 

Programme 2.2. Recherche d'un objet dans un vecteur 

Le programme est encore une application du schema suivant : 

Initialisation 
TANTQUE NON fini 
FAIRE avancer 
FAIT 
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L' initialisation de l'etat du monde dit que, le programme n'ayant rien vu (i:=0), 
l'objet n'a pas encore 6t6 trouve (trOuv&=FAUX). Pour avancer, on augmente i de 1 
et on remet trouve a jour. Le processus peut terminer soit par la reussite (trouve 
devient VRAI), soit par l'inspection du dernier element (i=n). Notons que les deux 
conditions peuvent etre verifiees en meme temps, car l'objet peut se trouver en t[n]. 
La boucle termine, car i augmente a chaque tour. L' index de la seule reference a t[i] 
est juste, car 1' initialisation et la condition de la boucle font que 0<i<n a l'entree 
dans la boucle. Avec Faugmentation de la valeur de i, ceci devient 0<i<n au moment 
de la rtsfifirence, c'est-a-dire l<i<n. 

Dans ce programme, nous avons introduit deux nouvelles notions : la pre- 
condition ("precondition") et la post-condition ("postcondition"). La precondition est 

ce qui doit etre vrai pour que le programme soit executable. Ici, il faut qu'il y ait au 
moins un element dans le tableau. Notons en fait que, si cette condition n'etait pas 
verifiee, par exemple n=0, le programme dirait que l'objet est absent, mais ceci est 
un hasard. La post-condition est ce que Ton sait & la fin du programme. Elle doit 
servir de specification du resultat. 

Prenons 1' exemple d'une procedure de calcul d'une racine carr6e. On aurait la 
structure suivante : 

DONNEE x: reel; 
PRECOND x^O; 
Procedure de calcul de y=sqrt(x) 

POSTCOND y . y = x (en general, d epsilon pres) 

On sait calculer les racines carries des nombres non negatifs (pre-condition) et 
le resultat est un nombre dont le carre est la donnee d'origine (post-condition). 

Dans notre programme de recherche lineaire, la pre-condition est raisonnable, 
mais la post-condition n'est pas entierement satisfaisante. Elle comporte deux 
lacunes : elle n'est pas directement deductible du texte et le sens de la variable trouve 
n'est pas bien defini, meme si le programme est juste. Regardons maintenant le 
processus de deduction. 

Reprenons le programme avec sa pre-condition, mais sans post-condition pour 
le moment : 

PRECOND n>0; 

DEBUT VAR i: entier, trouve: bool; 

i:=0; trouv6:=FAUX; 
TANTQUE kn ET NON trouve 
FAIRE i:=i+l ; trouve := objet=t[i] 

FAIT 

FIN 
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Que sait-on sur ce programme ? Une premiere deduction, evidente, est que, 
comme le programme a quitte la boucle, la condition de continuation de la boucle est 
maintenant fausse. Ainsi, apres FAIT, on peut deduire : 

NON (kn ET NON trouve). 

Par application de la loi de Morgan, ceci devient : 

i^n OU trouve. 

Comme i n'avance que d'un pas a la fois, on peut deduire : 

i=n OU trouve. 

Les OU de l'informatique etant inclusifs, les deux clauses peuvent etre vraies 
en meme temps (cas de t[n] = objet). Les post-conditions du programme donne 
suivent logiquement, en considerant l'affectation a la variable trouve dans l'interieur 
de la boucle. 

En general, on peut souhaiter que la post-condition se retrouve a partir de la 
pre-condition, de deductions du type donne ci-dessus et des definitions de l'etat du 
monde. Mais la post-condition que nous avons jusqu'ici est incomplete. Quand 
l'objet a ete trouve, i indique bien son emplacement, mais la non-reussite est mal 
decrite. Nous n'avons pas encore explique que NON trouve veut dire que l'objet est 
absent. 

L'erreur, classique, a ete de construire les pre- et post-conditions aprbs coup. II 
aurait fallu specifier le programme avant de l'ecrire, les conditions servant de 
specification. Dans notre cas, on aurait : 

PRECOND 

Soient un objet et un vecteur t de longueur n, n>0; 

POSTCOND 
Ou t[i]=objet, 

ou il n'existe pas i, 0<i^n, tel que t[i]=objet. 

Notons que la variable trouve ne figure plus dans la post-condition. Elle sert a 
distinguer les deux cas de solution. II faut maintenant demontrer la verite de cette 
nouvelle post-condition. Quand l'objet a ete trouve, la logique precedente est 
suffisante. Pour montrer l'absence de l'objet, il nous faut une clause de plus. 
Reprenons le texte avec de nouvelles decorations : 
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DONNEES objet, n: entier; t: TABLEAU [1 ..n] DE entier; 
PRECOND n>0; 

DEBUT VAR i: entier, trouve: bool; 
i:=0; trouv6:=FAUX; 
TANTQUE i<n ET NON trouve 
FAIRE INVARIANT 0<j«£i t[j] * objet; 

i:= i+1; trouve := t[i]=objet 
FAIT 

POSTCOND t[i]=objet OU "i (0<i^n => t[i]*objet) 
FIN 

L'invariant dit que la boucle n'est executee que si l'objet n'a pas encore ete 
trouve. La post-condition peut maintenant etre demontree. En effet, a la terminaison 
de la boucle, si trouvd est FAUX, alors i=n ET t[i] n'est pas l'objet. Mais avant 
d'executer l'instruction i.— aucun t(j], 0<j<i, n'etait l'objet. Apres l'augmentation 
de i, si t[i] n'est pas l'objet, l'invariant reste confirme. Quand i=n, la deuxieme partie 
de la post-condition est correcte. 

L'utilisation d'assertions telles que les invariants de boucles permet d'arriver a 
des preuves formelles de programmes. En pratique, la preuve complete d'un 
programme de taille industrielle s'avere longue et couteuse. Mais ce n'est pas une 
raison pour Fetudiant de ne pas connaitre ces techniques. Une familiarite avec les 
bases des preuves de programmes est une des cles de l'amelioration de la performance 
d'un programmeur. II cree de la sorte des schemas mentaux qui font qu'il analyse ses 
problemes de maniere plus rigoureuse, produisant des programmes de meilleure 
qualite en commettant moins d'erreurs, 

2.2.2. Un piege 

Les programmeurs de la generation precedente n'aimaient pas les variables 

bool&nnes. Us preferaient ecrire le programme ci-dessus dans une forme a priori plus 
simple : 

i:=0; 

TANTQUE i<n ET t[i+l ]*objet 

FAIRE i:=i+l 

FAIT 

En supprimant la variable trouve, le test de reussite se fait sur le prochain 
element t[i+l]. Afin d'eviter le calcul de F index, on peut redefinir i comme F index de 
l'61ement k traiter et non plus l'element qui vient d'etre traite : 
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i:=1; 

TANTQUE \<J\ ET t[i]*objet 

FAIRE i:=i+l 

FAIT 

Malheureusement, ces deux programmes comportent une faiblesse qui ne se 
montrera qu'avec certains compilateurs. Supposons, dans la demi6re version ci- 
dessus, que l'objet n'est pas present dans le vecteur. La boucle toume pour la derniere 
fois avec i=n. L'execution de la boucle augmente la valeur de i, donnant i=n+l. On 
teste de nouveau avec TANTQUE. La premiere partie de la condition (i<n) est fausse, 
mais on peut quand meme en evaluer la deuxieme partie afin de disposer des deux 
operandes de Foperateur ET. Cette deuxieme partie comporte une reference a t[i], 
c'est-a-dire a t[n+l]. Or, cet element n'existant pas, le programme peut terminer par 
l'erreur "index en dehors des bornes". 

Cette faiblesse peut s'eliminer avec l'introduction de Foperateur ETPUIS 
("CAND"), qui ordonne les tests. Ainsi, on pomrait 6crire : 

TANTQUE ten ETPUIS t[i]*objet 

Quand i>n, c'est-a-dire que le premier operande est FAUX, on n'evalue pas le 
deuxieme. Cela depend du fait que (FAUX ET b) est toujours FAUX, quelle que soit 
la valeur de b. II existe egalement Foperateur OUALORS ("COR"). Les definitions 
de ces operateurs sont les suivantes : 

a ETPUIS b J SI a ALORS b SINON FAUX FINSI 
a OUALORS b J SI a ALORS VRAI SINON b FINSI 

Notons que ces deux definitions sont celles qui se trouvent dans beaucoup de 
livres de logique pour ET et OU, mais ces operateurs ne sont pas mis en ceuvre de 
cette fajon dans tous les compilateurs. 

2.2.3. La dichotomie 

Quand on ne sait rien sur les elements d'un tableau, pour etablir qu'une valeur 
donnee ne s'y trouve pas, il faut inspecter tous les elements, car la valeur peut 
figurer & n'importe quelle place. Maintenant, nous allons considerer un cas plus 
interessant, ou les elements du tableau sont ordonnes. C'est comme Faeces a un 
annuaire telephonique. Pour trouver le numero de M.Dupont, on ne regarde pas 
toutes les entrees de A a DUPONT, On precede par des approximations. 

Pour nous faciliter la programmation, nous allons proceder par des 
approximations simples. Dans un annuaire de 1000 pages, on regarde la page 500. 
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Si l'element recherche est alphabetiquement plus petit, on a restreint la recherche aux 
pages 1 a 499, ou, s'il est plus grand, aux pages 501 a 1000. Chaque regard coupe le 
domaine de recherches en deux. Les recherches s'arretent si l'element examine est le 
bon, ou si le domaine de recherches est devenu nul (l'objet est absent). 

Considerons une premiere version de ce programme, avec Petat du monde 
suivant : 

bas, haut: entier, SI t[i]=objet ALORS bas £ haut 
centre: entier, t[centre] est I'§l6ment d examiner 
trouve: booleen, trouve 1 t[centre]=objet 

Appliquons comme d'habitude le schema : 

Initialisation 
TANTQUE NON fini 
FAIRE avancer 
FAIT 

Le programme 2.3 est une solution possible. 



DONNEES n: entier, t: TABLEAU [1 ..n] DE entier; 
PRECOND n>1 ET (0<i<Un => t[ikt[j]j 
DEBUT VAR bas, haut, centre: entier, trouve: bool; 
bas:=1 ; haut:=n; trouv6:=FAUX; 
TANTQUE haut-bas>1 ET NON trouve 
FAIRE centre :=entier((haut+bas)/2) ; 
CHOIX t[centre]<objet: bas:=centre, 
t[centre]=objet: trouv6:=VRAI, 
t[centre]>objet: haut:=centre 
FINCHOIX 

FAIT; 

SI NON trouve 
ALORS SI t[bas]=objet 

ALORS centre:=bas; trouv6:=VRAI 
SINON SI t[haut]=objet 

ALORS centre:=haut; trouv6:=VRAI 
FINSI 

FINSI 

FINSI 

FIN 

POSTCOND trouve => t[centre]=objet, 
NON trouve =>" i, O&n, t[i]*objet 

Programme 2.3. Recherche par dichotomie 
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La pre-condition exprime le fait que nous disposons d'un vecteur ordonne d'au 
moins deux elements. L' initialisation indique que le domaine de recherches est 
t[l..n], l'objet n'etant pas encore trouve. La condition apres TANTQUE merite des 
explications. La proposition NON trouve est evidente, mais la partie avant le ET 
Test moins. Pour considerer une valeur intermediate, nous imposons qu'elle soit 
differente des deux extremes, c'est-a-dire que I'in6galit6 suivante soit vraie : 

Proposition A. bas < centre < haut 

Cette indgalite necessitant l'existence d'au moins une valeur entre bas et haut, 
on retrouve : 

Proposition B. haut ■ bas > 1 

Dans la boucle, le calcul de la valeur de centre necessite une conversion, par la 
fonction entier, du resultat de la division, qui est un nombre reel, en un nombre 
entier. On confirme facilement que centre est bien entre haut et bas en appliquant la 
proposition B ci-dessus. Quand la somme (haut + bas) est impaire, la division par 
deux donne un nombre reelle de la forme n,5. Que l'arrondi vers un entier donne n ou 
n+1 n'a aucun effet sur l'algorithme (les deux possibilites respectent la proposition 
A). Par la suite, la clause CHOIX force un choix entre les trois possibilites ouvertes 
apres comparaison de l'objet avec t[centre]. Ou l'objet a ete trouve (t[centre] = objet), 
ou le domaine de recherches est coupe en deux (bas:=centre ou haut:=centre, suivant 
la condition). 

Si l'objet est trouve, tout va bien. Si la boucle se termine sans trouver l'objet, 
haut et bas sont maintenant deux indices successifs : 

haut = bas + 1 

La dernifcre partie du programme teste si l'objet se trouve en tfhaut] ou en 
t[bas]. Ce test est irritant pour le programmeur. II n'est utile que si t[l]=objet ou 
t[n]=objet, car des que haut ou bas change de valeur, ils prennent celle de centre, 
t[centre] n'etant pas l'objet. Done le test ne sert que si Tune des variables haut OU bas 
a garde sa valeur initiale. 

En fait, le test est du a une faiblesse de specification. II faut decider si oui ou 
non t[haut] ou t[bas] peuvent contenir l'objet, et cela de maniere permanente. 
Essayons deux nouvelles versions du programme. La premiere respecte les 
initialisations de P original, ce qui veut dire que t[haut] OU t[bas] peut toujours 

contenir l'objet : 
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DONNEES n: entier, t: TABLEAU [1 ..n] DE entier; 
PRECOND n>1 ET 0<i<jm => t[i]^t[j] 
DEBUT VAR bas, haut, centre: entier, trouve: bool; 
bas:=1 ; haut:=n; trouve:=FAUX; 
TANTQUE hautebas ET NON trouve 
F Al R E centre :=entier((haut+bas)/2) ; 

CHOIX t[centre]<objet: bas:= centre + 1, 
t[centre]=objet: trouve:= VRAI, 
t[centre]>objet: haut:= centre - 1 
FINCHOIX 

FAIT 

FIN 

POSTCOND trouve => t[centre]=objet, 

NON trouve => " i, 0<\<f\, t[i];*objet 

Programme 2.4. Dichotomie, version 2 



Trois changements figurent dans cette version du programme par rapport a 
l'original. Dans le choix, quand tfcentre] n'est pas l'objet, la nouvelle valeur de bas 
(ou haut) ne doit pas se situer sur le centre, mais un pas plus loin, sur le premier 
candidat possible (centre est le dernier element rejete). En plus, dans le TANTQUE, 
au lieu de proposition B, nous avons simplement : 

haut £ bas 

Quand la proposition B est vraie, la situation n'a pas chang6 par rapport a la 

version originale. Quand haut=bas, centre prend cette meme valeur et Ton teste le 

dernier candidat. Si ce n'est pas le bon, on ajuste haut ou bas, avec le resultat : 

bas > haut 

Quand haut suit immediatement bas, centre va prendre une valeur qui est soit 
celle de haut, soit celle de bas (il n'y a pas d'espace entre les deux). Mais, grace au 
fait que bas, ouhaut, prend sa nouvelle valeur un cran plus loin que le centre, an 
prochain tour de la boucle, on aura : 

bas = haut ou bas > haut 

La boucle termine toujours. La preuve de la correction de ce programme peut 
se faire a partir des assertions suivantes : 
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Al. 0<i<bas => t[i]<objet 
A2. haut<i^n => t[i]>objet 

La condition de terminaison bas > haut montre alors l'absence de l'objet dans t. 
Cela correspond a la deduction que Ton peut faire a la sortie de la boucle : 

NON (haut > bas ET NON trouve) 

C'est-a-dire : 

bas > haut OU trouve 

bas > haut est la condition d'absence, trouve implique t[centre] = objet. 

La deuxieme bonne solution a notre probleme est de faire en sorte que ni t[bas] 
ni t[haut] ne peut etre l'objet. Considerons la version du programme 2.5. 

DONNEES n: entier, t: TABLEAU [1 ..n] DE entier; 
PRECOND n>1 ET 0<i<kn => t[i]^t[j] 
DEBUT VAR bas, haut, centre: entier, trouve: bool; 
bas:=0; haut:=n+l ; trouve:=FAUX; 
TANTQUE haut-bas>1 ET NON trouve 
FAIRE centre:=entier((haut+bas)/2); 
CHOIX t[centre]<objet: bas:=centre, 
t[centre]=objet: trouve:=VRAI, 
t[centre]>objet: haut:=centre 
FINCHOIX 

FAIT 

FIN 

POSTCOND trouve => t[centre]=objet, 

NON trouve => " i (0<kn => t[i]*objet) 

P rogramme 2.5. Dichotomie, version 3 

Dans cette version, seules les initialisations de bas et haut ont ete modifiees. 
Ces variables prennent des valeurs d'indices inexistantes. Mais ce n'est pas grave, car 
les 616mentS correspondants ne seront jamais references. La preuve de cette version 
du programme est facile. Elle est la meme que celle de la version originale, sauf pour 
rarrivee a la situation : 

haut = bas + 1 
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Dans ce cas, comme ni t[haut], ni t[bas] n'est l'objet, celui-ci n'est pas dans t, 
car il n'y a pas d'element entre t[haut] et t[bas]. 

La technique consistant a utiliser comme butee une valeur inexistante (ici en 
considerant [0..n+l] au lieu de [l..n]) est utilisee frequemment par de bons 
programmeurs. Notons qu'il n'a pas ete necessaire de creer reellement les elements 
fictifs introduits. 



2.3. De la complexite des algorithmes 

Le parallele avec l'annuaire telephonique montre que les performances de ces 
deux algorithmes (recherche lineaire et recherche dichotomique) sont tres differentes. 
Leur complexite est facile a etablir. 

Pour la recherche lineaire : 

• Si l'objet est present dans le tableau, il peut etre n'importe ou. En moyenne, 
le programme examine n/2 elements avant de trouver le bon. 

Si l'objet est absent, on examine tous les n elements. 
La complexite de l'algorithme est done de o(n). Cela veut dire que si Ton doublait le 
nombre d'61ementS, les recherches dureraient en moyenne deux fois plus longtemps. 

La dichotomie est tout autre. Ici, un doublement de la taille du tableau 
necessite un seul pas supplemental (chaque examen d'un element divise la taille du 
domaine de recherche par deux). Considerons le cas d'un tableau de 2 k elements. 
Apres un pas, on a 2 k_1 elements a considerer, apres deux pas, 2 k " 2 , .... apres k pas, 
un seul element (2°). La dichotomie a done une complexite de o(log2 (n)). 

La difference entre ces deux courbes est tres importante, surtout quand le 
nombre d'elements est grand. La table 2.1 compare le nombre maximum de pas (n) 
dans le cas de recherche lineaire avec le nombre maximum de pas avec la dichotomie. 



n dichotomie 

10 4 (2 4 = 16) 

100 7 (2 7 = 128) 

1000 10 (2 10 = 1024= Ik) 

1 000 000 20 (2 10 = 1024k = 1 048 576) 



Table 2.1. Comparaison entre la recherche lineaire et la dichotomie 
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Dans le jargon de l'informatique, on parle de la reduction d'un probleme en n 
vers un probleme en n-1 (recherche lineaire), ou vers un probleme en n/2 
(dichotomie). Par la suite, dans le chapitre 3 sur les tris, nous verrons que dans 
certains cas, on reduit un probleme en n vers deux problbmes en n/2. On appelle 
cette demiere technique I 'a it de diviser pour rejier ("divide and conquer"). 

Nous tirons deux conclusions de cette breve discussion. La premiere est que des 
connaissances minimales sur le calcul de la complexite des algorithmes sont 
n6cessaires si Ton veut pratiquer de l'informatique a un bon niveau. Le sujet, assez 
difficile, est tres etudie par des chercheurs. Ces recherches demandent surtout des 
competences elevees en mathematiques. Mais on peut faire des estimations utiles 
avec un bagage limite. La deuxieme conclusion concerne l'efficacite des programmes. 
On voit souvent des programmeurs s'echiner sur leurs programmes pour gagner une 
instruction ici ou la. Evidemment, dans certains cas precis, cela peut devenir 
necessaire, mais e'est rare. Le gain d'efficacite est limite. Mais le probleme n'est pas 
le meme au niveau des algorithmes, comme l'attestent les chiffres de la table 2.1. 
Des gains d'efficacite a travers l'algorithmique peuvent etre importants. II n'est pas 
tres sense d'optimiser un mauvais algorithme — • mieux vaut commencer avec un 
bon. L' optimisation locale peut se faire par la suite en cas de besoin. 

2.4. Resume des principes introduits 

Au cours de ce chapitre, nous avons introduit plusieurs principes importants, 
qui forment la base de notre technique de programmation. 

Nous travaillons souvent a partir d'un schema de programme ("program 

scheme"). C'est' une maquette qui indique la structure generate. Le seul schema 
utilise jusqu'ici est celui d'un processus lineaire : 

Initialiser 

TANTQUE NON fini 

FAIRE avancer , 
FAIT 

On complete le schema lineaire en utilisant un etat du monde, qui est une 

definition precise des variables. Cet etat permet de confirmer la justesse du 
programme en l'exploitant comme une liste a cocher en reponse a des questions du 
type suivant : 

• Est ce que l'etat du monde est completement initialise avant la boucle ? 

. Quel est l'effet de 1' operation avancer sur chaque element de l'etat du monde ? 

L' existence d'une definition precise des variables facilite la definition de la 
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condition de terminaison. De meme, expliciter la notion d'avancer diminue la 
probability de l'ecriture de boucles infinies. Neanmoins, pour eviter cette 
mesaventure, on demontre consciemment (et consciencieusement) la terminaison de 
chaque boucle. Le r^flexe de demonstration de validite doit aussi jouer chaque fois 
que l'on repere un element de tableau. On demontre que les indices sont 
n&essairement entre les bornes. 

Nous avons utilise les notions de pre-condition et post-condition. La pre- 
condition indique les limitations du programme, c'est-a-dire les caracteristiques des 
donnees en entree. Elle sert pour des demonstrations de correction, mais aussi pour la 
documentation d'un programme. Avec la pre-condition, un utilisateur eventuel peut 
confirmer qu'il a le droit d'appeler le programme avec les donnees dont il dispose, 

La post-condition indique ce que le monde exterieur sait apres F execution du 
programme. On doit pouvoir remplacer tout programme par n'importe quel autre qui 
respecte les memes pre-condition et post-condition, sans que 1'utilisateUT eventuel 
s'en rende compte. 

Une partie difficile du processus de la mise en oeuvre est la specification du 
programme. Mais le travail fait a ce niveau est payant. En effet, le cout dune erreur 
augmente avec le temps qu'elle reste presente. Mieux vaut passer un peu plus de 
temps en debut du processus que de payer tits cher, plus tard, l'elimination des 
eneurs. 

Pour demontrer la correction d'un programme, on utilise des assertions et des 
invuriants. Les assertions sont des formules logiques qui sont vraies aux endroits ou 
elles figurent dans le programme. Un invariant est une assertion qui est vraie a 
chaque tour d'une boucle. Une boucle est completement definie par son etat du 
monde et son invariant. Certains auteurs incluent l'etat du monde dans l'invariant. 
Les deux techniques sont equivalentes. 

2.4.1. Un aparti sur les preuves de programmes 

Les techniques resumees ci-dessus reprennent des notions emanant des travaux 
sur la preuve de programmes. Le but est d'impregner les cerveaux des etudiants de 
mfcanismes mentaux allant dans cette direction, sans passer a une approche trop 
rigoureuse pour etre soutenue dans la pratique. On doit savoir pourquoi le 
programme marche, sans avoir explicite tout le developpement mathematique. 

Pour le puriste, ou le mathematicien, cette approche n'est pas satisfaisante. II 
serait normal — dans leur monde ideal que tout programme soit accompagne d'une 
preuve formelle. Ce sont les imperatifs economiques qui font que le monde n'est pas 
ideal, surtout en acceptant les capacites et les motivations des programmeurs. 
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Le style presents est done un effort de compromis, maitrisable par les eUldiantS 
dont nous disposons tout en les guidant. Au cours de leurs travaux diriges, ils 
menent a bien au moins une preuve formelle complete afin de comprendre les outils 
sous-jacents. 

2.4.2. Le style d'ecriture 

Cet ouvrage s'adresse aux problemes algorithmiques. Aucun des programmes 
presenter ne depasse la taille d'un module dans un programme complet. La mise 
ensemble d'unites de programme pour la construction d'un produit industriel est un 
probleme aborde ailleurs (cours de genie logiciel, projets). 

Mais il ne faut pas perdre de vue cette question de modularite. L' utilisation de 
la clause DONNEES, avec les pre-conditions et post-conditions, vise, parmi d'autres 
buts, a preparer la definition de modules, avec leurs specifications, interfaces et 
corps. En pratique, dans le cours enseigne, ces notions forment la matibre de 
discussions, preparant ainsi le travail en profondeur a venir. 

2.5. Adressage disperse 

Les premiers programmes dans ce chapitre sont des exemples simples, 
introduits pour des raisons pedagogiques. Mais, en meme temps, nous avons 
examine des methodes de recherche d'un element dans un vecteur. Dans une premiere 
liste d'algorithmes de ce type, il faut introduire celui de l'adressage disperse ("hash 
code"), qui est frequemment utilise dans la gestion, dans les compilateurs ou dans 
l'intelligence artificielle. La methode a ete inventee pour accelerer des recherches de 
positions jouees dans le jeu de dames [Samuels 1959]. 

Supposons que nous voulons creer un annuaire telephonique au fur et & mesure 
de l'arrivee de numeros connus, sans faim de tri k chaque arrivee. C'est ce que l'on fait 
habituellement dans un carnet d'adresses. Dans un tel carnet, pour eviter d'avoir a 
examiner tous les noms de personnes, on les divise en 26 classes, en fonction de la 
premiere lettre du nom. Dans le carnet, on commence une nouvelle page pour chaque 
lettre. Les recherches vont plus vite parce que les comparaisons ne sont faites qu'avec 
les noms ayant la meme premiere lettre. La premiere lettre sert ici de cte ("key"). 
Une cl6 est une fonction des caracteres constituant le mot qui sert k diviser 
l'ensemble de mots possibles dans des classes. Si les noms Ctaient distribues de 
maniere egale entre les classes, on divise le nombre de comparaisons par le nombre 
de classes. Pour un carnet, on ne regarde qu'un nom sur 26. Notons que ce rendement 
n'est pas atteint en pratique, parce que, par exemple, les pages K, X, Y, . . . 
contiennent moins de noms que certaines autres. 
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L'adressage disperse est une mise en oeuvre du principe du carnet, avec quelques 
changements en raison des caracteristiques des ordinateurs. En particulier, un carnet 
comporte beaucoup de lignes vides sur les pages moins remplies. Nous decrivons 
deux versions de l'algorithme d'adressage dispersde, une premiere avec un nombre de 
Cl6s pl us petit que la taille de la memoire disponible, et une deuxieme ou le nombre 
de cl6s est le meme que le nombre de cases disponibles. 



2.5.1. Algorithme avec chainage 

Dans cet algorithme, le nombre de cl6s est plus petit que le nombre de cases 
disponibles en memoire. Considerons le carnet, ou il y a 26 cles. On cree un tableau 
du type dontl6 dans la figure 2.2. 



Index Norn Suivant 



2 




3 




4 


DUPONT 


5 




6 




7 




8 
9 




10 

11 




12 




13 




14 




15 




16 




17 




18 




19 




20 


TOTO 


21 




22 




23 




24 


X 


25 


Y 


26 




27 


DURAND 


28 


T I NT I N 


29 


DUPOND 



Figure 2.2. Table pour l'adressage disperse avec chainage 
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Dans ce tableau, les 26 premieres cases sont reservees pour le premier nom 
recu de chaque classe (en supposant que la cle est la premiere lettre). La colonne de 
droite contient 1' index de l'entree contenant le prochain nom de meme cle. Un 
successeur (suivant) d' index -1 indique que le nom dans cette entree est le dernier 
rencontr^ dans sa classe (il n'a pas de successeur). Une case vide a egalement -1 
comme successeur. Les noms comportent 8 caracteres, etant completes par des 
espaces. Une case vide comporte 8 espaces. La figure montre l'etat de la table apres 
la lecture des noms suivants : 

X, DUPONT, TOTO, Y, DURAND, TINTIN, DUPOND. 

A la reception d'un nom, on calcule sa cle i, ici la premiere lettre. Differents 
cas se presentent : 

- Le nom est la premier occurrence d'un nom avec cette cle. Alors, la case 
d'index i est vide. Le nom s'insere h cette place et le processus est termine. 

■ Un nom de cle i a deja 6t6 rencontre. On compare le nouveau nom avec celui 
d'index i dans la table. Si les deux sont identiques, le nom est trouve et le processus 
se termine. 

- Si les deux noms sont differents, il faut examiner les successeurs eventuels 
comportant la meme cle. Un successeur de valeur • 1 indique que la liste est terminee. 
On ajoute le nouveau nom £ la premiere place disponible et le processus est termine. 

• Si le successeur existe, son index est donne dans la colonne correspondante. 
Oncomparelenouveaunomavecle successeur, en se ramenant au cas precedent. 

Cela donne l'algorithme du programme 2.6, page ci-apres, appele a l'arrivee de 
chaque occurrence d'un nom. 

A la fin de cet algorithme, la variable adresse contient 1' index du nom dans la 
table. La variable pi indique la premiere case de debordement libre dans la table (avec 
la valeur de 27 au depart de l'algorithme). L'algorithme ne traite pas le probleme d'un 
debordement eventuel du tableau. 
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DONNEES cle: PROC(chaTne) -> entier; 
nom: chaTne(8); 

cars: TAB [1 ..taille] DE chaTne(8); 
succ: TAB [1 ..taille] DE entier; 
DEBUT VAR i, adresse: entier; 
trouve: bool; 
pi: entier INIT 27; 
i:=cl6(nom); 
S I cars[0 = " " 

ALORS cars[i]:=nom; % complete par des espaces % 

adresse:=i; trouv6:=VRAI 
SINON trouv6:=FAUX; 

TANTQUE NON trouve 
FAIRE SI nom = cars[i] 

ALORS adresse:=i; trouv6:=VRAI 
SINON SI succ{i] = -1 

ALORS cars[pl]:=nom; % avec des espaces % 
succ[i]:=pl; succ[pl]:=- 1; adresse.^pl.- 
trouve^VRAI; pl:=pl+1 
SINON i:=succ[i] 
FINSI 

FINSI 

FAIT 

FINSI 

FIN 

Programme 2.6. Adressage disperse 



2.5.2. Autant de des que de cases 

L'etablissement d'un chainage entre les differents noms d'une meme cle gaspille 
de la memoire. Pour eviter cette depense, on peut choisir une fonction qui donne 
autant de cles que de cases dans le tableau. En restant avec notre cle (peu realiste) de 
la premiere lettre, on peut refaire la table de la figure 2.2 pour obtenir celle de la 
figure 2.3. 
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Index Nom 



4 DUPONT 

5 DURAND 

6 DUPOND 

7 
8 
9 

10 

11 

12 

13 

14 

15 

16 

1 7 
18 
19 

2 0 TOTO 

2 1 TINTIN 

22 

23 

24 X 

25 Y 
26 

Figure 2.3. Table sans zone de debordement 



Cette figure reprend la situation de la figure 2.2. On note que les noms 
supplementaires commencant par D ont pris les places des noms commencant par E 
et F. Que se passe-t-il alors si Ton rencontre le nom ESSAI ? On comparera ESSAI 
avec DURAND (case 5), puis avec DUPOND (case 6), avant de decouvrir que la case 
7 est vide. ESSAI rentrera dans la case 7. La cle sert & etablir un point de depart des 
recherches, donnant ainsi l'algorithme du programme 2.7. 
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DONNEES cle: PROC(chame) -> entier; 
nom: chaTne(8); 
tab: TAB [ 1 ..taille] DE chaTne(8); 
DEBUT VAR i, adresse: entier; 
trouve: bool; 
i:=cle(nom); trouve:=FAUX; 
TANTQUE NON trouve 
FAIRE SI tab[i] = " " % Case vide % 
ALORS trouve :=VRAI; tab[i]:=nom 
SINON SI tab[i]=nom 

ALORS trouve :=VRAI 

SINON SI i=taille % Derniere case du tableau % 
ALORS i:=l % Recommencer en haut % 
SINON i:=i+l 
FINSI; 

SI i=cle(nom) ALORS table pleine FINSI 

FINSI 

FINSI 

FAIT 

FIN 

Programme 2.7. Adressage disperse^ version 2 



La table doit etre circulaire (on recommence a regarder en haut si Ton arrive a la 
fin). Si la table est pleine, la recherche d'un nouveau nom fait tout le tour du 
tableau. Sauf si la table est pleine et le nom est absent, a la fin du programme, i 
COntient l'index du nom dans le tableau. 



2.5.3. Choix de cle et efficacite 

Jusqu'alors, nous avons utilise comme cle la premiere lettre du nom. Cette cle 
n'est pas la meilleure dans la plupart des cas. Le choix de la bonne cle depend des 
caracteristiques des noms que Ton va rencontrer. Dans les compilateurs, on utilise 
souvent comme cle la somme des representations internes des caracteres (code ASCII 
ou EBCDIC) modulo la taille de la table. La taille choisie est habituellement line 
puissance de deux afin de calculer le modulus par decalage sur les bits. Ce choix est 
assez bon pour les identificateurs dans des programmes. II permet d'utiliser le 
deuxi&me algorithme, sans chainage. 

Le deuxieme algorithme est assez efficace tant que la table ne se remplit pas. 
Son rendement devient mauvais si la table est presque pleine. Le premier algorithme 
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garde ses qualites jusqu'au remplissage du tableau, au prix d'une occupation de 
memoire superieure. 

Rappelons que les deux algorithmes necessitent une fonction qui distribuent 
bien les noms parmi les cles. Si cette distribution est mauvaise, les deux 
algorithmes se rapprochent de 1'algorithme de recherche lineaire. 

Notons finalement que les algorithmes d' adressage disperse ne diminuent pas la 
complexite theorique des recherches par rapport a celle des recherches lineaires. 
L' amelioration de cette methode reside dans la diminution de la constante dans le 
formule. Toute la famille est d'ordre o(n) pour la recherche d'un nom. 



2.6. Exercices 

1. Avant de passer au chapitre qui donne la solution, chercher 1' amelioration du 
programme de recherche d'un mode de vecteur. Ce nouveau programme est plus court 
et plus efficace que l'ancien. En particulier, il comporte une variable de moins et un 
test (SI ALORS . . . ) de moins. 

2. Considerons une nouvelle version du programme de dichotomie, ecrit dans le style 
"sans booleens" : 

DEBUT VAR bas, haut, centre: entier; 

bas:=1 ; haut:=n; centre :=entier((n+1)/2); 
TANTQUE haut>bas ET t[centre]*objet 
FAIRE SI t[centre]<objet 

ALORS bas:=centre 
SINON haut:=centre 
FINSI; 

centre :=entier((haut+bas)/2) 

FAIT; 

SI t[centre]=objet 

ALORS ... % c'est trouve % 

SINON ... % absent % 

FINSI 

FIN 

Cette version, qui est souvent proposee par des 61&ves, contient un piege. 
Lequel ? Comme mise sur la voie, on considerera le problfcme de la terminaison de la 
boucle. 
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3. Mise en ceuvre dans Fordinateur des algorithmes d'adressage disperse. On prendra 
des textes de programmes afin de disposer de suites de noms (identificateurs). En 
prenant des programmes de grande taille, on peut mesurer l'efficacite des differents 
algorithmes de recherche d'objets (lineaire, dichotomie, adressage disperse) en 
dressant des graphiques du temps d'execution contre le nombre de noms lus. On 
examinera aussi rinfluence du choix de la fonction de calcul des cles sur les 
algorithmes d'adressage disperse. 
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Les tris 



Ce n'est pas la premiere fois qu'un livre sur 1' algorithmique et la 
programmation aborde le sujet des tris ("sorting"), loin de la. Mais le sujet est 
essentiel — on ne peut pas s'appeler informaticien sans avoir quelques connaissances 
des algorithmes de base. En plus, la matiere est un excellent terrain d'entrainement. 
C'est done sans honte que nous abordons ce chapitre, meme si d'illustres 
pr6d6cesseurs ont trace la route. Parmi ceux-ci, accordons une mention particuliere a 
[Knuth 1973], pour un volume tout entier consacre aux recherches et aux tris. 

Trier un vecteur, c'est l'ordonner. On peut trier des entiers, des reels, des 
chaines de caracteres, . . .11 suffit qu'il existe une relation d'ordre entre chaque paire 
d'616ments, e'est-a-dire que, si a et b sont des elements, une et une seule des relations 
suivantes est vraie : 

a<b a = b a>b 

Les axiomes habituels sont v6rifi6s : 

a<b <=> fc»a 

(a<b ET b<c) => a<c 

a=b <=> b=a 

On peut egalement decreter que les elements du vecteur sont tous differents, 
evitant ainsi de considerer l'6galit6. Cette restriction n'apporte pas grand-chose pour 
la programmation. 



ALGORITHMIQUE et programmation 
3.1. Recherche du plus petit element 

Un premier algorithme, que nous ne recommandons pas pour des applications 
pratiques, consiste en la recherche du plus petit element restant. Ainsi, au premier 
tour, on recherche t[min], le plus petit element en t[l..n]. 11 devient t[l] dans le 
vecteur trie. Le plus simple est d'echanger t[ 1] avec t[min], ce dernier arrivant ainsi & 
sa place definitive. Le problfeme en n est reduit a un probleme en n-1, car il reste a 
trier le vecteur t[2..n]. Au deuxibme tour on recherche le plus petit element en t[2..n] 
pour l'echanger avec t[2], et ainsi de suite. On aurait pu prendre l'element le plus 
grand, en l'echangeant avec t[n] . . . 

Le programme 3.1 montre un schema d' algorithme applique a un vecteur 
d'entiers : 



DONNEES n: entier; 

t: TABLEAU [1 ..n] DE entier; 
PRECOND n>0; 

DEBUT MONDE i: entier, t[1 ..i] est trie; 
i:=0; 

TANTQUE i < n-l 

FAIRE i:=i+l ; 

trouver j tel que t[j] est le plus petit element dans t[i..n] 
6changer(t[i], tffl) 

FAIT 

FIN 

POSTCOND t[1 ..n] est trie, cad (0<i<j<n => t[i]<t[j]) 
Programme 3.1. ScMmatiu tri par recherchedu plus petit element 

La seule instruction necessitant une explication est la condition aprfcs 

TANTQUE. 11 est necessaire de trier n-l elements, le dernier, t[n], etant 

necessairement a sa place a la fin (tOUS les elements qui le precedent sont plus petits 
que lui par construction, sauf cas d' elements egaux). 

Reste a coder le contenu de la boucle dans ce texte. Ce sera une nouvelle 
boucle (programme 3.2). 
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MONDE k: entier, I'index du dernier element consid6r6 dans la 
recherche du plus petit; 

min: entier, t[min] est le plus petit element en t[i..k]; 
k:=i; min:=i; 
TANTQUE k<n 
FAIRE k:=k+l ; 

SI t[k]<t[min] 

ALORS min:=k 

FINSI 

FAIT 

echange(t[i],t[min]) 
Programme 3.2. Boucle interne 



Cette boucle ne necessite pas d'explication. Notons neanmoins que le 
programme accepte l'egalite. Si le plus petit element existe en plusieurs 
exemplaires, on prend le premier arrive, le restant (ou les restants) etant considere(s) 
au tour suivant. Le programme 3.3, complet, utilise les mondes decrits 
prec&temment 

DONNEES n: entier; 

t: TABLEAU [1 ..n] DE entier; 
PRECOND n>0; 

DEBUT VAR i, k, min, temp: entier; 
i:=0; 

TANTQUE kn-1 

FAIRE i:=i+l ; k:=i; min:=i; 
TANTQUE k<n 
FAIRE k:=k+l ; 
SI t[k]<t[min] 
ALORS min:=k 
FINSI 

FAIT; 

temp:=t[i]; t[i]:=t[min]; t[min]:=temp 

FAIT 

FIN 

POSTCOND t[l ..n] est trie 
Programme 3.3. Programme complet 
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L'echange a &6 programme de maniere evidente, avec une variable temporaire. 
La complexite de cet algorithme est simple a calculer. Au premier tour, n-1 
comparaisons sont necessaires pour trouver le plus petit element en t(l..n). Au 
deuxieme tour il en faut n-2, puis n-3, . . . Le nombre de comparaisons est done : 

1 + 2 + 3 + . . . + (n-l) = n*(n-1)/2 

Comme les autres operations (il y a n-l echanges) se font moins souvent, 
Falgorithme est d'ordre o(n 2 ). Le nombre de comparaisons est independant du 
contenu du vecteur. 

3.2. Tri par insertion 

Dans le programme du paragraphe precedent, apres i tours, le vecteur t[l..i] est 
trie, mais en plus, ses i premiers elements sont deja a leur place, e'est-a-dire : 

0<j<k<i => t[j]<t[k] 
i<k<n => t[i]<t[k] 

La premiere assertion dit que les i premiers elements sont tries, la deuxieme dit 
que les elements apres t[i] sont plus grands ou egaux & ceux deja tries, t[i] etant le 
plus grand des elements tries. 

Une nouvelle version du tri ne necessite pas cette deuxieme condition. Apres i 
tours, les i premiers elements sont tries. On considere t[i+l]. Cet element va etre 
insere & sa place en t[l..i+l], ce qui implique, en general, la recopie un cran 
plus loin de chaque element de t[l..i] qui est plus grand que t[i+l] (voir le 
programme 3.4). 
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DONNEES n: entier; 

t: TABLEAU [1 ..n] DE entier; 
PRECOND n>0; 
DEBUT VAR i, j, temp: entier; 
arret: bool; 

MONDE i: t[l ..i] est trie; 

i:=1; 

TANTQUE i<n 

FAIRE i:=i+l ; 

MONDE j est la position du trou, 

temp = t[i], 

arret Ht[j-1]<temp; 
temp:=t[i]; j:=i; arret:= t[j-1]<temp; 
TANTQUE j>1 ET NON arret 
FAIRE SI t[j-1]>temp 

ALORS t[j]:=t[j-1]; j:=j-1 

sinon arret:=VRAI 
FINSI 

FAIT; 

t[j]:=temp 

FAIT 

FIN 

FOSTCOND t[1 ..n] est trie 
Programme 3.4. Tri par insertion 

Au tour i, ce programme fait "remonter" t[i] a sa place, ou plutot il fait 
descendre les elements de t[l..i-l] qui sont plus grands que t[i], afin de lui laisser sa 
place. On a done l'idee du "trou". On enleve l'element a considerer du vecteur, en le 
mettant dans la variable temporaire temp. Cela laisse un trou. On regarde le 
pr6d6cesseur du trou. S'il est plus grand que temp, il descend, e'est-a-dire le trou 
monte. Si le predecesseur n'est pas plus grand, ou s'il n'y a pas de predecesseur 
(l'element considere est le plus petit vu jusqu'ici), la recherche s'arrete, l'element en 
temp trouvant sa place dans le trou. 

Pour confirmer la validite des indices, notons que la boucle externe impose 
i<n. Apres i:=i+l, on a i<n, Dans la boucle interne, j>l (condition) et j<i 
(initialisation, avec j:=j-l dans la boucle). On deduit : 

1<j<n, done t[j] et t[j-1] existent. 

Les boucles se terminent, car i et j avancent, Fun vers n (i:=i+l), I'autre vers 1 

(N-D- 
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La preuve de ce programme est faite en demontrant que le nouvel dement 
retrouve bien sa place, ce qui revient a demontrer qu'a la fm de la boucle interne on 
a : 

0<k<j => t[k]<temp 
j<k<i => temp<t[k] 

Comme les elements Ctaient deja ordonnes, et comme le nouveau est 
compatible avec l'ordonnancement par les deux assertions ci-dessus, l'invariant de la 
boucle externe (t[ l..i] est trie) reste vrai. 

On pourrait supprimer la variable booleenne arret en utilisant l'operateur 
ETPUIS (programme 3.5). 



DONNEES n: entier; 

t: TABLEAU [1 ..n] DE entier; 
PRECOND n>0; 
DEBUT VAR i, j, temp: entier; 

MONDE i: t[l ..i] est trie; 

TANTQUE i<n 

FAIRE i:=i+l ; 

MONDE j est la position du trou, 
temp:=t[i]; j:=i; 

TANTQUE j>1 ETPUIS t[j-1]>temp 
FAIRE tD]:=t[j-1]; j:=j"1 

FAIT; 

t[j]:=temp 

FAIT 

FIN 

POSTCOND t[ 1 ..n] est trie 
Programme 3.5. Version avec ETPUIS 



On peut aussi faire la meme chose avec une butee (voir exercice a la fin du 
chapitre). La complexite de cet algorithme depend du nombre d' elements qui 
"descendent" a chaque tour. Avec une distribution aleatoire, la place de l'element 
COnsid6r6 va etre en moyenne au milieu des autres. Ainsi, le nombre d'elements & 
d£placer est : 

(1 + 2 + 3 + . . . + n-l)/2 = n*(n-1)/4 
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On est toujours en ordre o(n 2 ), mais cet algorithme est moins mauvais que le 
pr6c6dent (sans pour autant etre recommande dans le cas g6n6ral). En effet, si le 
tableau est d6ja trie ou presque, le premier algorithme fait le meme nombre de 
comparaisons que dans le cas d'une distribution aleatoire, tandis que dans le second, 
on constate n fois de suite qu'aucun mouvement n'est necessaire, descendant ainsi a 
l'ordre o(n). Cet algorithme est done bon pour des vecteurs que Ton sait d6ja tries, ou 
presque. 

3.3. Tri par bulles 

Bien connu aussi, ce tri a une mauvaise reputation du point de vue de 
l'efficacid, reputation qui n'est pas tout a fait justifiee. A chaque parcours du vecteur, 
on compare successivement chaque paire de voisins. Si leur ordre n'est pas le bon, 
on les dchange. Une des raisons de la mauvaise reputation de F algorithme est que 
certains enseignants montrent la version du programme 3.6, qui est effectivement 
particulierement inefficace. 

DONNEES n: entier; 

t: TABLEAU [1 ..n] DE entier; 
PRECOND n>0; 
DEBUT VAR fin, i: entier; 

MONDE fin: t[l ..fin] reste a trier, 

fin :=n ; 

TANTQUE fin>1 

FAIRE MONDE i: on va comparer t[i] et tp+0; 
i:=l ; 

TANTQUE Win 

FAIRE ASSERTION 0<j<i => t[j]<t[i]; 

SI t[i+l]ct[i] 

ALORS echange(t[i],t[i+l]) 

FINSI; 

i:=i+l 

FAIT; 

ASSERTION 0<j<fin => t[j]<t[fin]; 

fin:=fin-1 

FAIT 

FIN 

POSTCOND t[l ..n] est trie. 
Programme 3.6. Tri par bulles primitif 
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Cette version de l'algorithme est demontrable ft partir des assertions donnees 
dans le texte. A la fin de chaque tour, l'616ment "le plus lourd" est passe a la fin de la 
zone considdrge. C'est une version moins efficace de l'algorithme du §3.1, en 
trouvant le plus grand element au lieu du plus petit. 

II y a deux facons d'ameTiorer cette version : 

. Considerons le vecteur suivant : 

(23456789101) 

A chaque passage dans le vecteur, le seul echange se fera entre la valeur 1 et son 
pr£d6cesseur immediat. Nous proposons de faire des passages alternativement de 
gauche ft droite, puis de droite ft gauche. Dans ce cas, le deuxieme passage ramenerait 
le 1 au debut du vecteur, qui serait done trie apres deux passages. 

• Mais cela ne suffit pas, car ralgorithme ne se rend pas compte que le vecteur 
est trie. Considerons le vecteur suivant : 

(21345678910) 

On voit qu'il sera tri6 ft la fin du premier passage. Ce qu'il faut ajouter au programme 
est la realisation du fait que, etant donne qu'aucun echange n'a eu lieu, pendant ce 
premier passage, ft partir de la valeur 3, t[2..10] est trie et les 616mentS sont deja ft 
leurs places definitives. La preuve de cette affirmation vient du fait que le test 
demontre que ces 616ments sont ordonnes deux ft deux. L'ordonnancement etant 
transitif (a^b ET b<C => tic), ils sont done tous ordonnes. En plus, l'assertion 
montre que le dernier objet ft prendre sa place (ici t[2]) est le plus grand de tous ceux 
vu jusqu'alors. II en resulte que Ton peut r6duire Fespace ft trier plus rapidement. 

Ces ameliorations, qui prennent en compte le travail accompli "en route", 
e'est-a-dire les echanges intermediaires, donnent lieu a ralgorithme du programme 
3.7, 
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DONNEES n: entier, 

t: TABLEAU [1 ..n] DE entier; 
PRECOND n>0; 
DEBUT VAR bas, haut, i: entier; 

MONDE bas, haut: t[bas..haut] reste a trier; 

INVARIANT t[l ..bas], t[haut..n] sont ordonnes, 
t[l ..bas-1], t[haut+l ..n] sont places; 

bas:=1 ; haut:=n; 

TANTQUE bas<haut 

FAIRE MONDE i: t[i..haut] a comparer deux a deux, 

der: dernier objet bouge; 
INVARIANT t[der..i] est ordonne; 
i:=bas; der:=i; 
TANTQUE khaut 
FAIRE SI t[i+1]<t[i] 

ALORS <§change(t{i],t[i+1]); der:=i 

FINSI; 

i:=i+l 

FAIT; 

MONDE i: t[bas..i] a comparer, 

der: dernier objet bouge; 
INVARIANT t[i..der] est ordonne; 
haut:=der; i:=haut; 
TANTQUE bas<i 
FAIRE SI t[i]<t[i-1] 

ALORS echange(t[i],t[i-1]); der:=i 

FINSI; 

i:=i-1 

FAIT; 
bas:=der 

FAIT 

FIN 

POSTCOND t[1..n] est trie 

Programme 3.7. Tripar bulles normal 

Dans ce texte, ordonne veut dire que le vecteur indique est trie, e'est-a-dire que 
ses elements sont dans l'ordre, mais qu'ils pourront changer de place dans le vecteur 
final en fonction d'autres elements ft inserer. Un vecteur place est non seulement 
ordonne, mais ces elements occupent les places qu'ils auront dans le vecteur final. 
Ainsi : 
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t[i..j] est ordonne <=> (i3«l<j => t[k]<t[l]) 
t[Ljj est plac6 <=> (tp.-fl est ordonne ET 

(0<k<i,i<l<j=>t[k]<t[l])ET 

(i<l<j, j<m<n => t[l]<t[m])) 

Les invariants permettent la demonstration directe de la justesse du 
programme. 

On peut comparer la complexite de cet algorithme avec celle de l'algorithme 
par insertion. Les deux menent a un meme nombre d'eehanges, que Ton calcule de la 
mani&re suivante : en prenant les 616mentS deux a deux, le nombre d'echanges est le 
nombre de fois oil une paire d'elements n'est pas ordonnee. Mais le nombre de 
comparaisons n'est pas le meme pour les deux algorithmes. Malheureusement, ni 
l'un, ni Fautre n'est le meilleur dans tous les cas. Les deux algorithmes sont du 
meme ordre de complexit6, et ils partagent les proprietes d'etre bons dans le cas d'un 
vecteur bien conditionne et d'etre mauvais pour un vecteur mal conditionne. 

3.4. Diviser pour regner 

Les trois algorithmes de tri etudies jusqu'ici sont tOUS d'une complexite d'ordre 

o(n 2 ). Dans chaque cas, un passage du vecteur reduit un problfime en n a. un 
probl&me en (n-1). Comme pour la dichotomie par rapport a la recherche lineaire, il 
est possible de faire mieux en r&hlisant un probleme en n a deux problemes en n/2. 

Le terme generalement employe dans ces cas est celui de diviser pour regner. On 
divise le vecteur en deux moities, on trie chaque moiti4 et on remet les deux moities 
ensemble. 

II existe deux groupes d'algorithmes dans cette categorie : les tris par partition 
et les tris par fusion. Le tri par partition permet de travailler "sur place", c'est-a-dire 
en gardant les 616mentS dans le vecteur en evitant de copier le vecteur entier dans un 
nouvel emplacement dans la memoire. Le hi par fusion, utilise du temps heroique 
des programmes de gestion sur support de bandes magnetiques, necessite la creation 
d'un deuxieme vecteur dans la memoire (eventuellement secondaire), vecteur qui 
recoit les elements dans le nouvel ordre. 

Pour des raisons evidentes, nous accordons plus d'importance, dans cet 
ouvrage, au tri par partition. Le tri par fusion est souvent pris comme exemple dans 
les cours de programmation avec le langage PROLOG. 

3.4.1. Diviser pour regner avec partition 

Dans le cas d'une distribution aleatoire de valeurs, surtout si n est grand, 
l'algorithme de diviser pour regner avec partition est le meilleur de ceux present& 
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dans cet ouvrage. Le principe est de partitionner les elements en deux classes en les 
comparant avec un element dit pivot. Tous les elements plus petits que le pivot vont 
se trouver a sa gauche, les autres se trouvant a sa droite. Considerons le programme 
3.8, qui met en ceuvre cette classification par rapport au pivot t[ 1], 

DONNEES n: entier, 

t: TABLEAU [1 ..n] DE entier; 
PRECOND n>0; 

DEBUT VAR pivot, COmp, bas, haut: entier; 
MONDE pivot: base de la comparaison, 

COmp: valeur a comparer au pivot, 

bas: index du "trou" gauche, 

haut: index du "trou" droit; 
pivot:=t[l]; comp:=t[n]; bas:=l; haut:=n; 
TANTQUE bas<haut 
FAIRE SI comp<pivot 

ALORS t[bas]:=comp; bas:=bas+l ; comp:=t[bas] 

SINON t[haut]:=comp; haut:=haut-1 ; comp:=t[haut] 

FINSI 

FAIT; 

t[bas]:=pivot 

FIN 

POSTCOND (ki<bas<j<n => t[i]<t[bas]<t[j], bas=haut, t[bas]=pivot. 
Programme 3.8. Partition d'un vecteur 

On CT6q deux trous dans le vecteur, en extrayant t[l], le pivot, et t[n], un 
element de comparaison. Si l'616ment de comparaison est plus petit que le pivot, il 
est place dans le trou gauche, sinon dans le trou droit. Un nouveau trou est ct€6 a 
cote de celui que Ton vient de remplir, en extrayant comme Pigment de comparaison 
le voisin du trou rempli. Le processus continue jusqu'a la rencontre des deux trous. 
On insere le pivot dans le trou (unique) qui resulte de cette rencontre, sachant que 
tous les elements plus petits que le pivot sont a sa gauche, les autres etant a sa 
droite. Notons que cette formulation permet l'existence de valeurs Cgales, une valeur 
egale au pivot 6tant mise a sa droite. 

La division du vecteur en deux parties a reduit le probleme de hi en deux SOUS- 
problemes. II nous reste a trier les elements a gauche du pivot, puis ceux a sa droite, 
le pivot 6tant & sa place (il ne bougera plus). L'algorithme est applique 
recursivement, c'est-a-dire que chaque moitie du tableau est de nouveau divisee en 
deux par comparaison avec un pivot a lui, et ainsi de suite. Apres avoir divise par 
deux un certain nombre de fois, il reste au plus un seul 616ment dans chaque classe 
(certaines sont vides), qui est automatiquement a sa place s'il existe. 
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Pour mener a bien cette operation recursive, nous avons besoin de parameter 
le programme ci-dessus, qui va devenir une procedure dans le programme 3.9. 

PROC tri(i, j); 

GLOBAL n: entier, t: TABLEAU [1 ..n] DE entier; 
SPEC i, j: entier; 

trier t[i..j] par la methode de diviser pour regner; 
PRECOND 0<i,j<n; 

DEBUT VAR pivot, comp, bas, haut: entier; 
SI j>i 

ALORS pivot:=t[i]; comp:=t[j]; bas:=i; haut:=j; 
TANTQUE bas<haut 
FAIRE SI comp<pivot 

ALORS t[bas]:=comp; bas:=bas+l; comp:=t[bas] 
SINON t[haut]:=comp; haut:=haut-1 I comp:=t[haut] 
FINSI 

FAIT; 

t[bas]:=pivot; tri(i, bas-l); tri(bas+1, j) 

FINSI 

FINPROC 

Programme 3.9. Insertion dans une procedure recursive 

Le terme GLOBAL indique que la variable n est declaree a l'exterieur de la 
procedure (dans le programme englobant). Ainsi, la valeur de n est la meme pour 
chaque appel de la procedure, tandis qu'il existe un nouvel exemplaire de i et de j 
pour chaque appel (chacun possede le sien). Le tableau t existe egalement en un seul 
exemplaire, manipule par chacun des appels de la procedure. 

Cette procedure mene a bien le tri complet. Elle comporte le programme ecrit 
pnScedemment, avec les changements necessaires pour trier t[i..j] au lieu de t[l..n]. 
Une fois que la division en deux zones a eu lieu, il reste a les trier, l'une apres 
Fautre, par deux nouveaux appels de tri avec des parametres appropries : 

tri(i, bas-l) ET tti(bas+l , j) 

La procedure teste par i<j qu'il y a au moins deux elements a trier. Ce test 
confiie que les appels recursifs ne peuvent pas constituer un ensemble infini, car le 
nombre d'elements a trier diminue avec chaque appel. Evidemment, on trie le vecteur 
complet par un appel de la procedure de la forme suivante : 
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tri(1,n) 

La specification dune procedure (SPEC) prend la place de la clause DONNEES 
d'un programme. On y trouve la definition du jeu de parambtres, avec une indication 
de ce que la procedure doit faire. Cette indication est la post-condition de la 
procedure. 

L'algorithme a une complexite theorique d'ordre o(n*log 2 (n)), comme le 
montre le raisonnement suivant : 

Le cas parfait de cet algorithme arrive quand n=2'-l et le pivot divise toujours la 
zone en deux parties de longueurs egales. A la fin du premier tour, on a un pivot 
place et deux zones de taille (2'-l-l)/2 = 2'"' -L On voit que le nombre de 
comparaisons dans chaque zone est sa taille moins 1 (le pivot). En remultipliant par 
le nombre de zones, on deduit que le nombre total de comparaisons par tour est 
successivement : 

n-l, n-3, n-7, n-15, . . . 

Le nombre de tours est i, c'est-a-dire log 2 (n+l). On obtient, pour des donnees 
bien conditionnees, la complexite theorique annoncee. 

Ce calcul suppose que le pivot divise chaque fois son monde en deux parties de 
failles comparables. Considerons maintenant un vecteur deja ordonne. Comme nous 
prenons le premier element pour pivot, un tour va require le probleme en n dans un 
probleme en 0 et un autre en n-l. Le processus est de nouveau d'ordre o(n^). Ce 
fdsultat, surprenant, montre que diviser pour regner, tout en etant le meilleur 
algorithme dans le cas d'une distribution aleatoire avec une grande valeur de n, est le 
plus mauvais pour un vecteur deja trie (ou presque trie). On peut toujours prendre le 
pivot au centre du vecteur pour ameliorer quelque peu l'algorithme dans ce cas. 

3.42. Solution sans appel ricursif 

Le dernier programme ci-dessus comporte deux appels recursifs de la proceclure. 
Si, pour une raison quelconque, on voulait enlever la recursivite, on utiliserait une 
pile ("stack" ou "LIFO • Last In, First Out"). Cette derniere sert a memoriser, 
pendant le tri d'une zone, les limites des zones restant a trier. Ainsi, on trie t[i..j], 
avec, au debut : 

i=1 ET j=n 

Apres un tour, avec p l'index de la position finale du pivot, on a la situation 
suivante : 
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t[l ..p-1] est a trier 
t[p] est plac6 

t[p+l ..n] est a trier. 

Les bornes 1 et p-1 sont mises sur la pile et le prochain tour redemarre avec : 
i=p+1 ET j = n 

Apres chaque partition, une paire de bornes est mise sur la pile, I'autre etant 
trait6e de suite. Apres un certain nombre de partitions, la taille de la zone & traiter est 
1 (ou 0), c'est-a-dire qu'elle est terminee. Dans ce cas, on recommence avec la 
premiere paire de bornes sur la pile, et ainsi de suite. Cette technique donne lieu au 
programme 3.10. 

DONNEES n: entier, 

t: TABLEAU [ 1 ..n] DE entier, 
PRECOND 0<i,j<n; 

DEBUT VAR i, j, pivot, comp, bas, haut, pi: entier, 
fini: bool, 

pile: TABLEAU [1 ..taille] DE entier; 
j:=1;j:=n;fini:=FAUX;pl:=1; 

MONDE fini: t[l ..n] est tri6 

pi: index de la premiere case libre dans la pile 
t[i..j] est a trier 

les zones indiquees dans la pile sont a trier 
TANTQUE NON fini 
FAIRE MONDE comme d'habitude 

TANTQUE j>i 

FAIRE pivot:=tp]; comp:=t[j]; bas:=i; haut:=j; 

TANTQUE bas<haut ^ 
FAIRE SI comp<pivot 

ALORS t[bas]:=comp; bas:=bas+1 ; comp:=t[bas] 
SINON t[haut]:=comp; haut:=haut-1; comp:=t[haut] 
FINSI 

FAIT; 

t[bas]:=pivot; pile[pl]:=i; pile[pl+1]:=bas-1 ; pl:=pl+2; i:=bas+1 

FAIT; 

fini:= pl=1; 
SI NON fini 

ALORS pl:=pl-2; i:=pile[pl]; j:=pile[pl+1] 
FINSI 

FAIT 

FIN 

Programme 3.10. Version avec pile 
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On reviendra sur les piles au cours du chapitre 4, en les utilisant souvent par la 

suite. 

3.4.3. Quelques commentaires sur la recursivite 

Dans le programme de recherche lineaire, ou dans les (ris primitifs, la base des 
algorithmes a £t& la reduction d'un problbme en n vers un probleme en (n- 1). Cette 
Induction est triviale a mettre en oeuvre; a chaque tour d'une boucle, on execute 
l'affectation : 

n:=n-1 

et le tour est joue. 

En passant a la recherche dichotomique, il s'agit de reduire un probleme en n 
vers un probleme en n/2. On cherche un objet dans le domaine d'indices [g..d]. Par 
rapport a un 616ment mediane d'index c on reduit le domaine soit a [g..c], soit a 
[c.d]. La mise en oeuvre est egalement triviale, s'agissant de l'une des deux 
affectations suivantes : 

d:=c OU g:=c 

Notons neanmoins que ces utilisations de l'affectation n'existent que pour des 
raisons de mise en oeuvre dans un ordinateur. Le raisonnement sous-jacent est de type 
recurrent. Par exemple, pour la recherche lineaire, on aurait pu ecrire la fonction 
suivante : 

chercher(objet, t[l ..n]): 
SI t[n] = objet 
ALORS trouve 

SINON chercher(objet, t[l ..n-1]) 
FINSI 

Cette ecriture fonctionnelle peut s'exprimer par l'utilisation d'une procedure 
r&ursive ou par une boucle avec affectation. Toute boucle peut se reecrire en forme 
de procedure recursive de maniere directe. 

Considerons maintenant le tri par diviser pour r6gner avec partition. Ici, on 
r6dllit un problbme en n vers deux problemes en n/2. Une forme reduite de la 
fonction de tri est la suivante : 
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tri(t[g..d]): 
SI d>g 

ALORS partition(t[g..d], c); tri(t[g..c-1]); tri(t[c+l ..d]) 
FINSI 

Une simple affectation a chaque tour d'une boucle ne suffit plus, car, si 
l'affectation relance le premier des appels de tri resultant, les parametres du deuxibme 
appel doivent etre stockes quelque part afin de revenir dessus a la fin du premier 
appel. Notons que chaque appel provoque, a son tour, deux nouveaux appels, jusqu'a 
l'arriv6e de vecteurs de longueur 1 (ou 0). 

On peut envisager la cascade d'appels dans la forme d'un arbre (figure 3.1). 



tri(g..d) 




tri(g.-cl) tri(d..d) 




tri(g..c2) tri(c2..d) 




tri(g..c3) tri(c3..c2) 
Figure 3.1. Appels aprespartition 

Cette figure montre les appels en cours apres trois partitions, les indices ayant 
tit simplifies pour alleger le dessin (suppression des +1 et -1). Apres chaque 
partition, on reprend la branche gauche, en laissant la branche droite en suspens. Le 
successeur gauche est de nouveau partitionne, et ainsi de suite. 

La modelisation par procedure recursive est done une mise en ceuvre directe et 
simple de l'algorithme. Par la suite, nous avons monh"6 une mise en ceuvre avec une 
pile. La pile a servi a se rappeler ce qui reste a faire a un moment donne. Elle 
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contient, a chaque instant, les paires d' indices conespondant aux appels laisses en 
suspens. Notons que le programme du paragraphe preeddent trie la partie droite du 
vecteur, en mettant les indices de la partie gauche en suspens sur la pile. Dans la 
figure 3.1, l'ordre est inverse. En fait, l'ordre est indifferent; a la rigueur, on pomrait 
lancer les deux appels recursifs en parallele sur deux processeurs independents dans le 
cadre d'un materiel multiprocesseur. 

L'utilisation d'une pile est un moyen general pour mettre en ceuvre la 
recursivite. Si la pile est programmee explicitement, comme dans le paragraphe 
precedent, on y stocke les informations permettant de savoir ou Ton en est dans la 
cascade d'appels. On verra d'autres exemples dans le chapitre sur la marche arriere. 

En fait, lors de l'utilisation d'une procedure recursive, un compilateur engendre 
des instructions qui gerent une pile. Y sont gardees les valeurs des parametres de 
l'appel en cours et des variables locales a chaque niveau d'appel. Ce sujet est couvert 
avec plus de details dans des cours de compilation tels que [Cunin 1980]. 

On peut se demander quels sont les schemas recurrents qui permettent une 
traduction facile vers une boucle avec affectation directe de variables. Un probleme 
qui se reduit de maniere recmrente en deux sous-problemes, ou plus, ne permet pas 
une telle traduction, car il faut toujours garder la trace des appels en suspens. Dans 
un schema a un seul appel, la traduction vers une boucle est directe dans le cas d'une 
rtcursivitt terminals Dans ce cas, aucun calcul ne reste & executer apres la fin de 
l'execution de l'appel recursive, comme dans le cas de la recherche lineaire. Le 
problbme sera reconsidere dans le chapitre 7, qui traite de la transformation de 
programmes. 

3.4.4. Deux pivots 

Comme la division d'un problbme en deux sous-problemes — chacun de la 
moitie du cout du premier — repr6sente un grand gain d'efficacite, on peut se 
demander si Fidee de diviser un probleme en trois represente encore une amelioration. 
C'est purement au titre d'une speculation intellectuelle que la question se pose pour 
le tri par diviser pour regner avec partition. L'auteur ne connait pas de publication ni 
d'application de l'i(J6e. 

La mise en ceuvre d'un tel algorithme (programme 3.11) necessite done deux 
pivots, avec un element de comparison. Avec les deux pivots, on partage les 
elements en trois ensembles : ceux qui sont plus petits que le plus petit des pivots, 
ceux qui sont plus grands que le plus grand des pivots et ceux qui sont entre les 
deux. On trie par la suite les trois ensembles par trois appels recursifs. 
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DEBUT DONNEES n: entier; 

t: TAB [1 ..n] DE entier; x 
PROC tri(min, max: entier): \ 
VAR tg, td, pi, p2, t1, t2, t3: entier; i 
si max>min 

ALORS pi :=t[min]; t2:=min+1; p2:=t[t2]; tl :=min; 
test:=t[max]; t3:=max; 
TANTQUE t2<t3 
FAIRE SI test>t2 

ALORS t[t3]:=test; t3:=t3-1 ; test:=t[t3] 
SINON SI tesktl 

ALORS t[tl]:=test; tl :=tl+l; t[t2]:=t[t1] 
SINON t[t2]:=test 
* FINSI; 

t2:=t2+1; test:=t[t2] 
FINSI 

FAIT; 

t[t1]:=p1;t[t2]:=p2; 

tri(min, tg-1); tri(tg+1, td-1); tri(td+l, max) 

FINSI 
FINPROC; 
tri(1,n) 

FIN 

Programme 3.11 . Diviser pour regner avec deux pivots 



La seule difficult dans ce programme reside dans le maniement des trous et les 
ensembles pendant la partition. II y a trois trous : deux pivots et un 616ment de 
comparaison. Les trois trous definissent quatre zones dans le vecteur : trois 
ensembles a trier et une quatribme zone qui contient les elements en attente de 
partition. La figure 3.2 montre la distribution de trous et de zones pendant 
l'operation de partition. 



n 



t[i]<P1 P1<t[i]<P2 non traites t[i]>P2 

Figure 3.2. Vecteur en cours de partition 
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Dans cette figure, PI est la valeur du plus petit des deux pivots, P2 celle du 
plus grand. Les elements les plus petits occupent la premiere zone dans le vecteur, 
les plus grands occupant la demiere (quatrieme) zone. Les 616mentS ayant des valeurs 
interm6diaires occupent la deuxieme zone, laissant la troisieme pour les elements 
n'ayant pas encore 6t6 COnsidereS. 

Introduire un element dans la quatribme zone ne pose aucun probleme; il 
pousse le trou droit un cran a gauche, le dernier element non traite etant pris comme 
element de comparaison. De meme, on peut introduire un element a la fin de la 
deuxifcme zone, repoussant le trou central un cran a droite et en considerant le 
premier des elements non traites. Pour introduire un element dans la premiere zone, 
il faut "mordre" sur la deuxibme. Le premier 61ement de la deuxieme zone passe a la 
fin de cette zone, en poussant le trou central. Le premier trou peut ainsi avancer d'un 
cran, laissant de la place pour l'element a inserer. On considere le premier element 
non traite\ ejecte par l'avancee du trou central. 

Cette methode apporte-t-elle des ameliorations ? En fheorie, oui; la complexite 
est d'ordre o(n*log3(n)) au lieu de o(n*log 2 (n)). En pratique, il faut des vecteurs de 
tres grande taille pour s'en apercevoir. 

3.4.5. Tri par fusion 

Le tri par fusion date d'une epoque ou les memoires centrales etaient petites 
par rapport aux fichiers a trier, ces derniers etant stockes sur des bandes magnetiques. 
Le principe est simple : on coupe le fichier en deux moities, on trie chaque moitie, 
puis on fusionne les deux mditids triees en intercalant les valeurs de l'une et de 
l'autre dans le bon ordre. L'operation est repetee recursivement autant de fois que 
necessaire. L'operation de fusion necessite une deuxieme copie du fichier, ce qui 
double l'espace memoire occupe. En pratique, la fusion se faisait avec deux bandes 
magnetiques en entree, chacune avec une moitie triee du vecteur, et une troisieme 
bande en sortie pour recevoir le tout. 

La recursivite ne va pas jusqu'au bout; on divise le vecteur par deux jusqu'a ce 
que le r6sultat de la division puisse tenir dans la memoire principale. On peut des 
lors trier cette zone du vecteur en memoire avec un algorithme du type deja VU. Tout 
l'art de la programmation de ce type de situation comportant des bandes magnetiques 
consistait a bien organiser les informations sur les bandes afin d'eviter de couteux 
allers et retours pour rechercher le prochain bloc. 

En ignorant l'existence des bandes magnetiques, l'algorithme peut travailler de 
la manifcre suivante : 

■ On divise le vecteur tl autant de fois que necessaire pour que les zones soient 
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de longueur au plus 1. 

■ On fusionne des zones de longueur 1, par paires, pour creer des zones de 
longueur 2 dans une nouvelle copie t2 du vecteur. 

. Les zones de longueur 2 sont fusiondes pour creer des zones de longueur 4, 
de nouveau en tl. 

. On continue de la sorte jusqu'a la fin du processus. 

Prenons comme exemple un vecteur de 8 elements, trie a l'envers. Les etats 
successifs de fusion sont les suivants : 



t1: (87654321) Vecteur de depart 

t2: (78563412) Fusion d'elements creant des paires 

t1: (56781234) Fusion par paires 

t2: (12343678) Fusion par groupes de quatre 



Le programme 3.12 montre une mise en ceuvre correspondante. 
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DONNEES n: entier, t: TAB [1 „n] DE entier; 
DEBUTVAR I, pi, p2, p3, il, i2, i: entier; 

MONDE on fusionne des paires de zones de longueur I en zones de 
longueur 2*1; 

pi est I'index du premier Element de la premiere zone, p2 I'index 
du ctebut de la seconde zone, p3 de la zone suivante (n+l s'il 
n'y en a pas); 

i elements ont 6te recopi6s vers le nouvel exemplaire du 
vecteur, il est I'index du premier element non recopi£ de la 
premiere zone, i2 de la deuxieme; 

l:=l; 

TANTQUE l<n 
FAIRE i:=0; pi :=1; 
ANTQUE j<n 

faire p2:=p1+l; p3:=min(p2+l, n+1); il :=p1; i2:=p2; 

ANTQUE i<p3-1 
FAIRE i:=i+l; 
SI M=p2 

alors t2[i]:=t1[i2]; i2:=i2+1 
sinon si i2=p3 OUALORS !1[i1]<t1[i2] 
ALORS t2[i]:=t1[i1];il :=i1 +1 
sinon t2[i]:=t1[i2]; i2:=i2+1 
FINSI 

FINSI 

FAIT; v 

Pi :=p3 ' 

FAIT; 
l:=2*l; 
SI l<n 

ALORS pi :=1;i:=0; 
TANTQUE kn 

faire p2:=l+1;p3:=min(p2+l,n+1); M :=pi ;i2:=p2; 
TANTQUE i<p3-1 
FAIRE i:=i+1; 
si i1=p2 

alors t1[i]:=t[i2]; i2:=i2+1 
sinon si i2=p3 OUALORS t2[i1]<t2[i2] 
alors H[i]:=t2[i1]; ii:=ii+i 
sinon t1[i]:=t2[i2]; i2:=i2+1 
FINSI 

FINSI 

FAIT; 
p1:=p3 
FINSI; 

l:=2*l 

FAIT 

FIN 



Programme 3.12. Triparfusion 
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3.5. Resume de la complexite des algorithmes 

En fonction des discussions precedentes, on peut degager des raisons menant a 
un choix sense d'algorithme de tri pour une situation donnee. On aboutit aux 
conclusions suivantes : 

. La methode de recherche du plus petit element n'est jamais bonne. On ne 

l'utilisera pas 

. Pour de grands vecteurs dont les elements sont distribues aleatoirement, la 
methode de diviser pour regner avec partition est la meilleure. Elle est mauvaise dans 
le cas d'un vecteur deja trie, ou presque trie. Pour limiter les degats dans ce cas, il 
vaut mieux prendre le pivot au milieu de la zone, ce qui complique legerement le 
programme (voir exercice). 

- Pour des vecteurs presque tries, les deux m&hodes d' insertion ou par bulles 
sont de bonnes candidates, & condition d'utiliser la version optimisee pour la 
deuxieme. L'une ou l'autre peut etre la meilleure, en fonction des propri6t6s 
particul&res de l'ordonnancement approximatif deja existant. 

3.6. Exercices 

1. Dans le tri par insertion, l'utilisation de ETPUIS permet de supprimer la variable 
arret. Dans un langage sans ETPUIS, comment arriver au meme resultat par 
Fintroduction d'une butee ? 

2. Donner une regie permettant de calculer le nombre de comparaisons necessaires 
danslecas: 

• d'un tri par insertion, 
« d'un tri par bulles. 

3. Dans la derniere version du tri par diviser pour regner, on fera la modification 
consistant a prendre comme pivot l'etement au milieu de la zone a trier. 

4. Un bon exercice au niveau dune classe est de comparer les differents tris, mis en 
machine par differents eleves. On mesurera, par l'horloge et en comptant les 
operations, le cout de l'execution de chaque tri, en essayant une serie de vecteurs 
differents. On variera la distribution d'elements (tries, presque tries, aleatoire, tries a 
Ten vers, . ..) et la longueur du vecteur. L'etablissement de courbes pour chaque 
methode permet d'en confirmer la complexite theorique. 
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Des structures de donnees 



Les langages de programmation classiques permettent l'utilisation de variables 
simples (entiers, reels, . . .) et de tableaux, voire, eventuellement d'autres types de 
donnees composees (records, . ..). Dans le dernier exemple du chapitre 3 (tri par 
diviser pour regner), nous avons eu besoin d'une structure de donnees particuliere, la 
pile, qui n'existe pas directement dans ces langages. En effet, dans la plupart des 
langages de programmation, on ne dispose pas d'un type "pile". 

Un certain nombre de structures de donnees de ce type reviennent constamment 
dans les programmes. Dans ce chapitre, nous allons decrire les principals structures, 
avec les moyens de les representer dans des langages existants. Avec le 
ddveloppement de nouveaux langages, on peut s'attendre a voir ces objets devenir des 
types standards. En premier lieu, nous allons considerer les piles, les files ("queue" 
ou "FIFO - First In, First Out"), les arbres ("tree"), les treillis ("lattice") et les 
graphes ("graph"). 

4.1. Les piles 

Une pile est un ensemble ordonne d'objets de meme type (entiers, reels, . ..). 
C'est comme si Ton gardait une pile de livres sur son bureau. On peut poser un 
nouveau livre sur la pile, ou reprendre le livre qui est en haut de la pile. Extraire un 
livre du milieu est tellement difficile que nous renoncons a cet exercice. 

On dispose d'un vecteur dont les elements sont du type approprie. II existe deux 
operations fondamentales : poser un objet (empiler, ou "push") et en retirer un 
(depiler, ou "pull"). Le vecteur sert a stacker les objets poses. Pour savoir combien 
d'objets sont dans la pile a un moment donne, on utilise un pointeur de niveau. Le 
pointeur sert aussi a retrouver, dans le vecteur, le dernier objet depose. 



ALG ORITH M IQUE ET PROGRAMMATION 

Le programme 4.1 utilise l'index de la premiere case libre (pi) dans une pile 
d'entiers. 



DONNEES taille: entier, 

pile: TABLEAU [1 ..taille] DE entier; 
VAR pi: entier INIT 1; 
PROCEDURE empiler(x: entier); 
PRECOND pl<taille; 

pile[pl]:=x; pl:=pl+1 
FINPROC; 

PROCEDURE depiler(x:entier); 
PRECONQpl>1; 

pl:=pl-1 ; x:=pile[pl] 
FINPROC 

Programme 4.1. Mise en ceuvre d'une pile 



On voit que le nombre d'objets dans la pile a un moment donne est pl-1. Les 
procedures supposent que la pile ne deborde pas (PRECOND plltaille) et que l'appel 
de depiler n'a pas lieu avec une pile vide (PRECOND pl>l). En pratique, on teste 
ces conditions par des SI ALORS speeifiques en debut de procedure. 

Dans le chapitre precedent, nous avons programme directement la pile sans 
utiliser les procedures empiler et depiler. La taille de la pile etait supposee suffisante. 
Une pile vide etait le signe de terminaison de l'algorithme. 

La pile donnee ici est une pile d'entiers. Le meme schema de programme peut 
servir pour la construction d'une pile de reels, de booleens, . . . (il suffit de changer le 
type des elements du tableau pile). 

4.2. Les files 

Une file est une structure qui ressemble a une pile, a la difference pres que Ton 
retire 1' element le plus ancien au lieu du plus recent. Une telle structure est 
particulierement utile, par exemple, dans un systeme d' exploitation, pour la gestion 
des files d'attente. On pourrait imaginer la paire de procedures du programme 4.2. 
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PROGRAMME DISCUTABLE 

DONNEES taille: entier; file: TABLEAU [1 ..taille] DE entier; 
VAR ancien: entier INIT 1; libre: entier INIT 1; 
PROCEDURE mettre(x: entier); 
PRECOND libre<taille; 

file[libre]:=x; libre :=libre+1 
FINPROC; 

PROCEDURE enlever(x: entier); 
PRECOND ancien<libre; 

x:=file[ancien]; ancien :=ancien+l 
FINPROC 

Programme 4.2. Une file di scuta ble 

Pour mettre un objet dans la file, il faut qu'il y ait une place libre (PRECOND 
libreltaille). Pour en enlever, il faut qu'il y ait au moins un objet dans la file 
(PRECOND ancienclibre). Les deux indices ancien et libre "cement" les objets 
restants dans la file, qui se trouvent en file[ancien..libre-l]. 

Pourquoi avoir dit que le programme ci-dessus est discutable ? Tout 
simplement parce qu'il ne reutilise pas l'espace dans la file. Une fois arrive au bout 
(libre=taille+l), on ne peut plus y mettre de nouveaux elements, meme si le retrait 
d'autres elements a libere de la place. II faut rendre la file circulaire, en testant 
autrement la presence d'au moins un element (programme 4.3). 

DONNEES taille: entier, file: TABLEAU [Maille] DE entier; 
VAR n: entier INIT 0; ancien: entier INIT 1; libre: entier INIT 1; 
PROCEDURE mettre(x: entier); 

PRECOND n<taille; 
file[libre]:=x; 
SI libre=taille 
alors libre:=l 

SINON libre:=libre+l 
FINS1 
FINPROC; 

PROCEDURE enlever(x: entier); 
PRECOND n>0; 
x:=file[ancien]; 

SI ancien=taille 
ALORS ancien:=l 
SINON ancien:=ancien+l 
FINS1 
FINPROC 

Programme 4.3. Une file circulaire 
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La nouvelle variable n indique le nombre d' elements actuellement presents dans 
la file. Les pi€-conditions de non-pl6nitude et de non-vide dependent done de la valeur 
de n. Les variables libre et ancien avancent chaque fois de 1, en recommencant au 
debut de la file une fois arrivees a la fin. Dans un systeme d'exploitation, on parle 
souvent d'un producteur ("producer") et d'un consommateur ("consumer") a la place 
de mettre et enlever (voir [Griffiths 19881). 

4.3. Les arbres 

L'arbre est une structure fondamentale dans 1' algorithmique. Nous en avons deja 
vu un dans l'algorithme de tri par diviser pour regner. Le pivot sert a diviser le 
vecteur en trois zones : 

• une partie gauche, 
■ le pivot, 

• une partie droite. 

Les parties gauche et droite sont de nouveau divisees, chacune en trois 
nouvelles zones, et ainsi de suite. La figure 4.1 presente ces divisions en forme 
d'arbre. 



t[1 ..n] 




tf1-P-1] t[p] t[p+1..n] 




t[1 -pl-l] t[p1] t[p1+L.p-1] t[p+1..p2-1] t[p2] t[p2+1..n] 
Figure 41. Arbre du tri diviserpour regner 

Notons que les informaticiens ont pris l'habitude d'inverser les arbres : la 
situation initiale, la ratine ("root"), est en haut, et les 616mentS terminaux, les 
feu i 1 1 es ("leaf, leaves"), en bas. Des arbres australiens, en quelque sorte ... On parle 
done d'algorithmes qui descendent de la racine jusqu'aux feuilles. 
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Formellement, un arbre est une structure composee de noeuds ("node") et de 
branches (arc, "branch"). Une branche menant du nceud nl au noeud n2 fait que n2 est 
un successeur ("successor") de nl, nl etant le pridicessew ("predecessor") de n2. 

Chaque noeud, a F exception de la racine, possede un et un seul predecesseur. Un nceud 
peut avoir un nombre quelconque de successeurs. Un nceud sans successeur est une 
feuille. 

Suivre un chemin dans l'arbre se fait en suivant des branches successives. 
Ainsi, si n2 est un successeur de nl et n3 est un successeur de n2, alors il existe un 
chemin de nl a n3 (par l'intermediaire de n2). Nous limitons nos arbres & des arbres 
connexes, c'est-Mire que tout nceud n est accessible a partir d'une unique racine. 
Accessible veut dire qu'il existe un chemin de la racine de l'arbre jusqu'au nceud n. 
Notons que l'unicite des predecesseurs fait que ce chemin est unique. 

4.3.1. Arbres binaires et arbres n-aires 

La definition donnee ci-dessus permet a un nceud d'un arbre d' avoir un nombre 
quelconque de successeurs. En pratique, les programmeurs se limitent souvent a 
l'utilisation d'arbres binaires, dans lesquels un nceud a au plus deux successeurs. A 
priori, cela pourrait constituer une limitation du pouvoir d' expression, mais en fait 
ce n'est pas le cas. On demontre que-tout arbre n-aire peut etre represente sous la 
forme d'un arbre binaire, sans perte d' information. 

La forme binaire d'un arbre n-aire s'appelle un diagramme en forme de vigne 

("vine diagram"). Au lieu de garder, pour chaque noeud, des pointeurs vers ses 
successeurs immediats, on garde un pointeur vers son premier successeur (fils aine) 
et un deuxieme vers le prochain successeur de son pr6decesseur (frere cadet). La figure 
4.2 en donne un exemple. 

On voit sur la figure 4.2 les couples de pointeurs represented dans des boites k 
deux cases. Une case barree indique l'absence de successeur. Cette facon de decrire un 
arbre correspond d'assez pres a sa representation physique dans la memoire de 

rordinateur. 



71 



Algortthmique et programmation 



A 




Figure 4.2. Transformation d'arbre n-aire en arbre binaire 

4.3.2. Representation d'arbres 

Pour illustrer la representation physique des arbres binaires, reprenons les 
boites de la figure 4.2, en les mettant dans une table (en pratique on utilise deux 
vecteurs). Les noms des noeuds deviennent des entiers qui servent d' index dans ses 
vecteurs. En reprenant l'arbre de la figure 4.2, avec a=l, b=2, et ainsi de suite, on 
obtient les vecteurs de la figure 4.3. 



nom 


index 


succg 


succd 


A 


1 


2 


0 


B 


2 


5 


3 


c 


3 


0 


4 


D 


4 


8 


0 


E 


5 


0 


6 


F 


6 


0 


7 


G 


7 


0 


0 


H 


8 


0 


9 


1 


9 


0 


10 


J 


10 


0 


0 



Figure 4.3. Representation de l'arbre de la figure 4.2 
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Les deux premieres colonnes de cette figure ne sont donnees que pour faciliter 
la lecture. Dans la realite, seuls sont conserves les vecteurs succg (successeur 
gauche) et succd (successeur droit). 

Etant donne la possibility de transformer un arbre n-aire en arbre binaire, cette 
representation est toujours suffisante. La structure avec deux boites par nceud date de 
LISP [McCarthy 19601. Si Ton voulait conserver la forme originale, avec n 
successeurs possibles, la table de la figure 4.3 comporterait autant de colonnes que le 
nombre maximum de successeurs. 



4.3.3. Parcours d'arbres 

L'algorithme classique de parcours d'un arbre utilise un double appel recursif 
(programme 4.4). 

DONNEES succg, succd: TABLEAU [1 ..taille] DE entier; 
PROC parcours(n: entier); 
SI succg(n)*0 

ALORS parcours(succg[n]) 
FINSI; 

SI succd(n)*0 

ALORS parcours(succd[n]) 

FINSI 
FINPROC 

Programme 44. Parcours d'un arbre binaire 

Le parcours complet d'un arbre commence par un appel de 

parcours(racine) 

Le parcours d'un nceud implique le parcours de tous ses successeurs. Les tests 
sur les successeurs etablissent leur existence. 

Pour eviter de tester successivement l'existence d'un successeur gauche, puis 
celle d'un successeur droit, on peut avoir recours a la notion de butee. Un successeur 
inexistant a l'index 0. II suffit de creer un nceud d'index 0 comme butee (programme 
4.5). 
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DONNEES succg, succd: TABLEAU [OJaille] DE entier; 
PROC parcours(n: entier); 
SI n#0 

ALORS parcours(succg[n]); 
parcours(succd[n]) 

FINSI 

FINPROC 

Programme 4.5. Utilisation d'une butte 

Chaque fois qu'un noeud n'a pas de successeur (gauche ou droit), le programme 
appelle parcours(O). Le test au debut de la procedure fait que 1' execution ne va pas 
plus loin. 

Pour eliminer la recursivite, on peut utiliser line pile, dans laquelle on garde le 
souvenir des nceuds qui restent a parcouru (programme 4.6). 

DONNEES succg, succd: TABLEAU [Ctaille] DE entier; 

racine: entier; 
DEBUT VAR n: entier INIT racine, pi: entier INIT 1, 
fini: tool INIT FAUX, 
pile: TABLEAU [1 ..tpile] DE entier; 
MONDE pi: premiere case libre dans pile; 
TANTQUE NON fini 
FAIRE SI n=0 

ALORS fini := p|=1 ; 
SI NON fini 
ALORS depiler(n) 
FINSI 

SINON empiler(succd[n]); n:=succg[n] 
FINSI 

FAIT 

FIN j 
Programme 4.6. Parcours avec pile 



43.4. Parcours prefixe et post-fixe 

Considerons Farbre binaire complet a trois niveaux de profondeur qui se trouve 
en figure 4.4. Dans un arbre binaire complet, tOUS les noeuds autres que les feuilles 
ont exactement deux successeurs. De plus, tout chemin de la racine a une feuille est 
de la tneme longueur (trois dans notre cas). 
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1 




3 6 10 13 



A A A. A 

4 5 7 8 11 12 14 15 

Figure 44. Parcours en profondeur d'abord (ordre prtfixt) 

La numerotation des noeuds de cet arbre correspond a l'ordre de visite des noeuds 
dans un parcours dit "de profondeur d'abord" ("depth first search"). Nous verrons par 
la suite que cet ordre est celui de revaluation dite prifixie ("prefixed"). II est possible 
d'imaginer un parcours de Farbre en largeur d'abord ("breadth first search", figure 
4.5), mais le programme est plus difficile a ecrire (voir exercice en fin de chapitre). 




Algorithmique et programmation 

Revenons sur le parcours en profondeur d'abord, en decorant le programme avec 
une mention pour chaque visite d'un nceud (programme 4.7). 

PROCEDURE parcours(n); 
SI n>0 

ALORS visite-l ; 
parcours(sg[n]); 

visite-2; 
parcours(sd[n]); 
visite3 
FINSI 
FINPROC 

Programme 4.7. Triple visite d'un nceud 

Ces numeros de visite correspondent au fait que, dans un parcours d'arbre, on 
visite chaque nceud trois fois (figure 4.6), une fois en descendant a gauche, une fois 
ayant termine le parcours des successeurs gauche et avant de parcourir les SUCCesseurs 
droits, et une derniere fois en remontant a droite a la fin. 




Figure 4.6. Triple visite des nceuds d'un arbre 

Supposons qu'au cours du parcours en profondeur d'abord de l'arbre dans le 
programme 4.7, on imprimait le numero de nceud & la place de visite-l. Alors l'ordre 
d'impression des nceuds serait celui de la figure 4.4, ordre dit prefixe. Imprimer les 
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numeros des noeuds & la place de visite-2 donnerait l'ordre infixe ("infixed") (figure 
4.7), celui donne par visite-3 etant post-fix^ ("postfixed") (figure 4.8). Ces differentS 
ordres ont leur importance lors du traitement des expressions dans un compilateur. 




13 15 



Figure 4.7. Parcours en ordre infixi 




/\ /\ /\ As 

i 5 8 9 11 12 



9 11 



Figure 4.8. P a rco u rs en ordre post-fix^ 
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4.4. Les treillis 

Un treillis ressemble a un arbre, mais en permettant a un nceud d'avoir plus 
d'un pr6d£cesseur. Les branches peuvent done se rejoindre. Neanmoins, il n'est pas 
possible de boucler, e'est-a-dire que tous les chemins vont vers l'avant, en terminant 
sur une feuille. Nous retrouverons le probteme des boucles dans le paragraphe sur les 
graphes. La figure 4.9 donne un exemple de treillis. 



Racine 




Figure 4.9. Un treillis 

Notons, dans cette figure, que tous les chemins qui commencent a'la racine 
terminent sur Fune des feuilles Fl ou F2 dans un nombre fini de pas. 

On peut parcourir un treillis comme un arbre, par exemple avec un des 
programmes du paragraphe precedent. Ce parcours visitera certains noeuds plusieurs 
fois, car ils ont plus d'un predecesseur. En fait, chaque nceud sera visite autant de fois 
qu'il y a de chemins qui menent de la racine jusqu'a lui. II se peut que l'existence de 
visites multiples soit genante, pour des raisons d'efficacite ou de r6p6tition 
intempestive d'operations. Pour visiter chaque noeud une et une seule fois, on 
applique l'algorithme de parcours d'un graphe (voir ci-dessous). En fonction des 
besoins immediats, on traite done un treillis comme un arbre ou, le plus souvent, 
comme un graphe. 
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4.5. Les graphes 

Un graphe est toujours une collection de noeuds et de branches, mais sans 
limitations. Ceci permet des boucles, e'est-a-dire qu'il peut exister un chemin qui 
mene du nceud n jusqu'a lui-meme. Le reseau du metro parisien est un exemple de 
graphe, ou Ton peut revenir au point de depart d'un voyage. On voit qu'un treillis est 
un cas particulier du graphe, un arbre etant un cas particulier du treillis. 

Pour des raisons pragmatiques n'ayant rien a faire avec la theorie, les 
informaticiens limitent souvent les graphes qu'ils manipulent. Ainsi, on prendra 
souvent la forme de vigne, afin de ne parler que de graphes binaires. De meme, 
certains programmes ne travaillent que sur des graphes possedant une racine, e'est-a- 
dire que seront considered comme etant noeuds du graphe, les noeuds qui sont 
accessibles a partir de la racine. La notion de racine dans un graphe n'a pas de 
justification theorique; elle correspond a une volonte de simplifier les programmes de 
traitement. 

Dans le cas general, parcourir un graphe binaire a partir de sa racine avec 
l'algorithme donne pour les arbres menerait a une boucle, car l'algorithme suit tous 
les chemins possibles. Pour visiter chaque nceud une fois, sans boucler, on a recours 
a la technique de marquage. Le programme garde un drapeau booleen pour chaque 
noeud, indiquant si le noeud a d6j& ete visit6. On ne visite pas deux fois un noeud 
(programme 4.8). 

DONNEES taille: entier; 

succg, succd: TABLEAU [CLtaille] DE entier; 
racine: entier; 

DEBUT VAR marque: TABLEAU [Ctaille] DE bool INIT FAUX; 
PROCEDURE visiter(n: entier); 
SI NON marque[n] 
ALORS marque[n]:=VRAI; 

visiter(succg[n]); visiter(succd[n]) 
FINSI 
FINPROC; 
visiter(racine) 

FIN 

Programme 4.8. Visite d'un graphe 



Ce programme suppose l'existence d'un noeud de butee 0, qui est marque. 
Ainsi, le test de marquage et le test d' existence ne font qu'un. Notons que la 
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representation d'un graphe sur un couple de vecteurs succg et succd est la meme que 
celle d'un arbre. On peut enlever la recursivite de cet algorithme en utilisant une pile, 
exactement comme dans le cas d'un arbre. Par la suite, nous montrerons une methode 
de parcours d'un graphe sans recursivite mais aussi sans pile. 

En gdndral, on parcourt un treillis comme un graphe, c'est-a-dire en marquant 
les nceuds visites. Cela evite les visites r6p6t6es aux noeuds qui sont sur plusieurs 
chemins. 



4.51. Algorithme de Schorr-Waite 

L'algorithme classique de parcours d'un graphe, comme celui d'un arbre, met en 
jeu, soit un systeme d'appels recursifs, soit une pile. Une application frequente de 
l'algorithme se trouve dans des interpreteurs pour des langages tels que LISP, ou, 
quand la memoire affectee est pleine, il faut examiner les donnees en cours afin de 
recup6rer de l'espace dans la memoire principale. Pendant l'ex^Clltion d'un 
programme, on cree des noeuds et des arcs, on en supprime et on les modifie. Apres 
un certain temps, la partie de la memoire affectee au graphe est pleine. Mais, suite 
aux differentes manipulations, certains noeuds ne sont plus utiles. L'utilite d'un noeud 
est definie par son accessibilite. Tout noeud accessible a partir de la racine est utile (il 
peut toujours servir). Un noeud qui n'est pas accessible ne peut plus servir. Pour 
recuperer l'espace occupe par les noeuds inaccessibles, on commence par le marquage 
de tous les nceuds accessibles. Par la suite, dans une deuxieme phase, on examine 
tout l'espace utilisable, en recuperant les cases non marquees. Ce processus s'appelle 
le ramasse miettes ("garbage collection"). 

Le probleme vient du fait que le ramassage des miettes intervient au moment 
oil Ton constate que la memoire est pleine. On ne souhaite pas alors lancer une 
procedure recursive, etant donne que son execution necessite l'ouverture d'une pile de 
taille non previsible. Une solution de ce probleme existe [Schorr 19671. II s'agit de 
reutiliser les pointeurs pour mettre en ceuvre la notion de predecesseur. 

Considerons le parcours classique d'un graphe. La figure 4.6 a deja montre ce 
parcours sur trois nceuds d'un arbre. Chaque noeud est visite trois fois. Pour se 
rappeler ou il en est, l'algorithme de Schon-Waite garde dans la marque le nombre de 
visites au noeud deja effectuees: 

marque-0 noeud non visite 
marque=1 en train de visiter la descendance gauche 
marque=2 descendance gauche terminee, visite de la descendance 
droite en cours 

marque=3 le marquage de tous les successeurs du noeud a ete effectue 
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Avec cette representation de la marque, l'algorithme classique prend la forme du 
programme 4.9. 

PROC parcours(n); 
SI m[sg[n]]=0 
ALORS m[n]:=l ; 

parcours(sg[n]) 
FINSI; 

SI m[sdln]]=0 
ALORS m[n]:=2; 

parcours(sd[n]) 
FINSI; 
m[n]:=3 
FINPROC 

Programme 4.9. Marquage avec troisvaleurs 



En enlevant la recursivite de cette procedure, on peut deduire le programme 

4.10. 



DEBUT n:=racine; 

TANTQUE NON fini 
FAIRE CAS m[n] 
DANS 

0: m[n]:=1; 

SI m[sg[n]]=0 
ALORS n:=sg[n] 
FINSI, 

1 : m[n]:=2; 

SI m[sd[n]]=0 
ALORS n:=sd[n] 
FINSI, 

2: m[n]:=3; n:=pred[n] 
FINCAS 

FAIT 

FIN 

Programme 4.10. Version sans recursivite 



81 



Algorithmique et programmation 



Dans ce programme, la condition de terminaison est que tous les successeurs de 
la racine ont 6t6 marques. Dans ce cas, l'algorithme remonte a la racine pour la 
dernifere fois, c'est-a-dire que m[racine]=3. Le programme utilise toujours la notion de 
butee (un successeur absent pointe vers le nceud fictif 0, qui est marque), mais il 
teste si le successeur eventuel (droite ou gauche) est marque' avant d'y aller. 

La difficulte de cette version du parcours vient du fait que Ton ne sait pas 
"remonter" vers le predecesseur (n:=pred[n]) apres avoir passe la valeur de la marque a 
3. La structure des donnees ne comporte que des pointeurs vers les successeurs de 
chaque noeud. Pour r&SOUdre le probleme, Schon et Waite ont propose de renverser le 
pointeur vers le ^successeur suivi afin d'indiquer, de maniere temporaire, le 
predecesseur du noeud. Pour ce faire, on garde, a tout moment, dans une variable p, la 
valeur prec6dente de n, le noeud courant. Considerons la situation quand l'algorithme 
arrive pour la premiere fois sur le nceud 1 de la figure 5.6 : 

n=l, marque[1]=0, p=predecesseur(l), sg[1]=2, sd[1]=3. 

Appelons les cases qui indiquent les successeurs d'un nceud sg[n] et sd[n]. 
L'algorithme va proceder au marquage de la descendance gauche, n va done prendre la 
valeur de Sg[n], p prenant la valeur de n (il suit toujours n avec un temps de retard). 
Comme n a pris la valeur du successeur gauche, la case qui contient 1' index de ce 
successeur fait maintenant double emploi. On le reutilise temporairement pour se 
rappeler du predecesseur (l'ancienne valeur de p). On arrive a la situation suivante : 

n=2, marque[l]=l , p=1, sg[l]=predecesseur(l), sd[1]=3. 

A la fin du marquage de la descendance gauche, c'est-Mire quand marque[2] passe a 3, 
l'algorithme remonte sur le nceud 1 dans l'etat suivant, avec p toujours juste 
derriere n : 

n=l, marque[1]=1, p=2, sg[l]=predecesseur(l), sd[1]=3. 

On va maintenant repartir a droite, ayant restaure le successeur gauche et en 
conservant le predecesseur du noeud 1 dans sd[ 1] pendant la visite : 

n=3, marque[1]=2, p=1, sg[1]=2, sd[l]=predecesseur(l). 

En remontant de la droite, on aura : 

n=1, marque[1]=2, p=3, sg[1]=2, sd[1]=pr6d6cesseur(1). 

n va remonter vers le predecesseur, en laissant le nceud 1 dans son etat de 
depart: 
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n=pr6d<§cesseur{1), marque[1]=3, p=l, sg[1]=2, sg[1]=3. 

On retrouvera un resume de ces etats dans la figure 5.6. 

Avec cette introduction, nous arrivons au programme 4.11. 



DONNEES Un graphe code en forme de vecteurs comme avant 
DEBUT VAR n, p: entier; 

n:=racine; p:=0; 

TANTQUE m[racine]<3 

FAIRE CAS m[n] 
DANS 

0 : m[n]:=1 ; 

SI m[sg[n]]=0 

ALORS % La voie est libre a gauche % 

cycle(n, sg[n], p) 
SINON % La voie n'est pas libre, mais on fait comme si Ton y 
etait alle et comme si Ton en revenait % 

cycle(sg[n], p) 
FINSI, 

1 : m[n]:=2; 

SI m[sd[n]]=0 

ALORS % La voie est libre a droite % 

cycle(n, sd[n], sg[n], P ) 
SINON % Comme si Ton revenait de droite % 

cycle(sg[n], p, sd[n]) 
FINSI, 

2: m[n]:=3; cycle(n, sd[n], p) 
FINCAS 

FAIT 

FIN 

Programme 4.11 . Mise en oeuvre du predecesseur 



Le programme comporte deux points particuliers : 1' instruction cycle et ce qui 
se fait quand le chemin n'est pas libre, a droite ou a gauche. 

L'instruction cycle correspond a une affectation multiple. Ainsi, considerons le 
cycle suivant : 

cycle(n, sg[n], P ) 
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II correspond & 1' affectation multiple : 
(n, sg[n], p) := (sg[n], p, n) 

Pour executer une affectation multiple, on calcule les valeurs a droite, puis on 

les affecte, en parallble, dans les variables indiquees k gauche. Cette facon d'ecrire 

evite de considerer des effets de bord qui pourraient resulter de Fordonnancement des 

affectations individuelles. Par exemple, la sequence suivante d' affectations donnerait 
un resultat incorrect : 

n:=sg[n]; 
sg[n]:=p; - 
p: = n 

La premiere affectation a n fait que les deux autres affectations reperent la 
nouvelle valeur de n quand c'est Fancienne qui est recherchee. II faudrait sauvegarder 
l'ancienne valeur de n : 

xn:=n; 
n:=sg[xn]; 
sg[xn]:=p; 
p:=xn 

D'autres sequences de code donnent le meme resultat. 

Quand le successeur prevu est marque, soit parce qu'il n'existe pas, soit parce 
qu'il est en cours de traitement sur un autre chemin, le programme ne doit pas le 
prendre en compte. Mais les pointeurs doivent etre mis dans un etat qui correspond a 
F augmentation de la valeur de m[n]. L'algorithme met a jour les pointeurs comme si 
n revenait du chemin qui mene au successeur marque, bien qu'il n'y soit jamais alle. 

4.5.2. Amelioration de l'algorithme de Schorr-Waite 

Cette presentation vient de [Griffiths 1979], en reponse a celle de [Gries 1979]. 
C'est ce dernier qui a transmis l'idee de l'astuce qui permet de raccourcir le 
programme ci-dessus. Nous ne connaissons pas l'identite de l'inventeur de l'astuce. 

Considerons les etats par lesquels passent les noeuds au cours du traitement 
(figure 4.10). 
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m=0 
m=3 



m = l 



m=2 




Figure 4.10. Etats dans Schorr-Waite 

La nouvelle id& est de garder le pointeur vers le pred6cesseur dans la case du 
successeur que Ton ne prend pas. La figure 4.11 montre les etats resultants, dans 
lesquels le pointeur vers le descendant restant doit obligatoirement changer de place. 



m=0 
m=3 



m = l 



m=2 




Figure 4.11 . Etats dans Schorr-Waite ame'liore' 



85 



Algoriti imique ET programmation 



Cette id6e, a priori un peu bizarre, donne lieu au programme 4.12. 



DONNEES Un graphe code en forme de vecteurs comme avant 
DEBUT VAR n, p: entier; 

n:=racine; p:=0; 

TANTQUE m[racine]<3 

FAIRE CAS m[n] 
DANS 

0: m[n]:=l ; 
SI m[sg[n]]=0 

ALORS cycle(n, sg[n], sd[n], p) 
SINON cycle(sg[n], sd[n], p) 
FINSI, 
1 : m[n]:=2; 
SI m[sg[n]]=0 

ALORS cycle(n, sg[n], sd[n], p) 
SINON cycle(sg[n], sd[n], p) 
FINSI, 

2: m[n]:=3; cycle(n, sg[n], sd[n], p) 
FINCAS 

FAIT 

FIN 

Programme 4.12. Schorr-Waite ameliore 

On voit que, dans cette version du programme, les memes cycles se repetent. 
En regroupant, d'une part, les cycles identiques et, d'autre part, les augmentations des 
valeurs de m[n], on Obtient le programme 4.13. 

DEBUT VAR n, p: entier; 
n:=racine; p:=0; 
TANTQUE m[racine]<3 
FAIRE m[n]:=m[n]+l ; 

SI m[n]=3 OUALORS m[sg[n]]=0 

ALORS cycle(n, sg[n], sd[n], p) 

SINON cycle(sg[n],sd[n], p) 

FINSI 

FAIT 

FIN 

Programme 4.13. Version condensie 
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Cette version finale du programme, tres condensee, est assez agreable du point 
de vue de l'esthetique, mais elle a demande un travail considerable. C'est un bon 
exercice de style. En plus, le programme resultant peut etre repris tel quel et installe 
dans un ramasse-miettes performant. 



4.5.3. Representation de graphes sur une matrice binaire 

On peut aussi representer un graphe sur une matrice binaire, ce qui facilite le 
traitement de graphes n-aires. Un 1 a l'intersection de la ligne Y et de la colonne X 
veut dire que le nceud X est un successeur du nceud Y. Considerons le graphe de la 
figure 4.12. 



A 




F 



Figure 4.12. Un graphe 

Sa representation en forme de matrice binaire se trouve sur la figure 4.13. 





A 


B 


c 


D 


E 


F 


A 


0 


1 


1 


0 


0 


0 


B 


0 


0 


1 


0 


1 


0 


c 


0 


0 


0 


1 


0 


0 


D 


0 


1 


0 


0 


0 


0 


E 


0 


0 


0 


0 


0 


1 


F 


0 


0 


0 


0 


0 


0 



Figure 4.13. Representation sur une matrice binaire 
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Un 1 represente un arc qui mene du noeud indique par la ligne a celui de la 
colonne, un 0 indiquant l'absence d'un tel arc. Cette reprdsentation est plus compacte 
que la table utilis6e auparavant dans le cas d'un graphe permettant de multiples 
SUCCesseurs (graphes n-aires). Elle se prete aussi a certaines classes de manipulations, 
comme dans le paragraphe suivant. 



4.5.4. Fermeture transitive 

La fermeture d'une relation est une operation bien connue des mathematiciens. 
Elle rtJSlllte de la transitivite de l'operation consid6r£e. Par exemple, cette propriete 
nous permet de deduite: 

si a<b et txc 
ALORS a<c 
FINSI 

Intuitivement, on dit que "les amis des amis sont des amis". Considerons la 
notion de successeur dans un graphe. Dans la matrice du paragraphe precedent, un 1 
dans la position (x,y) indique qu'il existe un arc entre le noeud y et le noeud x (x est 
successeur de y). Maintenant, si x est un successeur de y, et si y est un successeur de 
z, on peut deduire que x est un successeur de z (en passant par y), c'est-a-dire qu'il 
existe un chemin entre z et x. La fermeture transitive d'un graphe donne sa 
connectivity. Pour chaque noeud n, on obtient la liste complete des noeuds qui sont 
accessibles h partir de n. 

La figure 4.14 montre la fermeture transitive du graphe de la figure 4.13. On 
voit, dans la premiere ligne, qu'il est possible, en partant de A, d'atteindre les noeuds 
B, C, D, E, F. En revanche (premi6re colonne), A n'a pas de predecesseur. On voit 
bien que A est la racine d'un graphe connexe. De la meme facon, dans la derniere 
ligne, on voit que F n'a pas de successeur, c'est-a-dire que F est une feuille. 



A B C D E 

A 0 1 1 11 

B 0 1 1 1 1 

c 0 1 1 1 1 

D 0 1 1 1 1 

E 0 0 0 0 0 

F 0 0 0 0 0 0 

Figure 414. Fermeture transitive du graphe de la figure 4.13 
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Les noeuds B, C, D prdsentent la particularite d'etre leurs propres successeurs 
(un 1 sur la diagonale principale). C'est une indication de la presence dune boucle 
(on peut aller de B jusqu'a B). Le graphe n'est done pas un treillis et, a fortiori, n'est 
pas un arbre. 

L'algorithme direct de calcul dune fermeture transitive suit la definition : "SI B 
est un successeur de A, ALORS tous les successeurs de B sont egalement 
successeurs de A". Considerons la ligne A de la figure 4.13. Elle comporte des 1 
dans les positions B et C. Ainsi, les successeurs de B sont des successeurs de A, 
comme le sont les successeurs de C. Pour les successeurs de B, on prend la ligne B 
en l'ajoutant a la ligne A par l'operation OU inclusive. On fait de meme pour chaque 
1 de la ligne A, que le 1 soit d'origine ou qu'il resulte dune addition. On ne reprend 
pas deux fois un 1 dans la meme position. Pour la ligne A, cela donne les operations 
de la figure 4.15. 



ligne A au depart : 

addition de la ligne B : 

addition de la ligne C 

addition de la ligne D 

addition de la ligne E : 

addition de la ligne F : 



A 

0 
0 
0 
0 
0 
0 



B 


C D 


E 


F 


1 


1 0 


0 


0 


1 


1 0 


1 


0 


1 


1 1 


1 


0 


1 


1 1 


1 


0 


1 


1 1 


1 


1 


1 


1.1 


1 


1 



(sans changement) 
(sans changement) 



Figure 4.15. Fermeture transitive, une ligne 

Aucun nouvel 1 n'ayant ete introduit, on a termine avec la ligne A. La meme 
operation a lieu pour chaque ligne du graphe. Cet algorithme n'est pas optimal. II a 
ete ameliore & plusieurs reprises [Warshall 1962], [Griffiths 1969]. 

4.6. Ordres partiels et totaux 

Les quatre structures, les files, les arbres, les treillis et les graphes, ont des 
proprietes d'ordonnancement importantes. Dans une file, il existe un ordre total, 
c'est-a-dire que pour chaque couple (x, y) d' elements, une et une seule des relations 
suivantes est vraie (pos est la position de l'element donne) : 

pos(x) < pos(y) 
ou pos(x) > pos(y) 
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Cette propri6t6 fait qu'il n'existe qu'un seul ordre d'ecriture dans lequel les 
616ments se trouvent chacun a leur place. Elle est evidente dans le cas d'une file (voir 
le cas des entiers positifs). 

Dans un arbre ou dans un treillis, il n'existe que des ordres partiels. Un ordre 
doit respecter la regie suivante : 

SI x est successeur de y 
ALORS pos(x) > pos(y) 
FINSI 

Consid6ronS Farbre binaire trivial de la figure 4.16. 



1 




2 3 

Figure 416. Arbre binaire simple 

La regie d'ordonnancement est respectee paries deux sequences suivantes : 

123 
132 

La figure 4.4 (prefixe ou en profondeur), comme la figure 4.5 (en largeur), 
presentent des ordres partiels dans un arbre. Les ordres des figures 4.7 (infixe) et 4.8 
(post-fix6) ne respectent pas la regie. Les arbres et les treillis permettent toujours les 
ordonnancements partiels des elements les constituant (il existe un seul ordre, e'est-a- 
die un ordre total, si et settlement si aucun nceud n'a plus d'un successeur direct). 

L'ordonnancement n'est plus possible pour un graphe qui comporte une boucle. 
En effet, comme la boucle indique qu'il existe au moins un element qui est son 
propre successeur, il n'est plus possible de respecter la regie d'ordonnancement par la 
relation de succession. 

4.7. Exercices 

telARQUE. — La notion de pointeur n'a pas ete introduite. II faut done gerer 
des tables d'index. 

1. On considere des families de type un peu particulier. II y a un seul parent, avec au 
plus deux enfants. Un enfant peut lui-meme etre parent d'autres enfants, et ainsi de 
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suite. Le tout consfitue done un arbre binaire dont les nceuds sont decorSs (dans le 
champ val) avec le nom de la personne representee. Un nom est une suite de lettres 
de 8 caract&res au plus. La racine de l'arbre s'appelle "adam" (ou "eve", si Ton 
pr6ffere). On ecrira les procedures recursives suivantes : 

. PROC trouveKnom, n) -> (bool, entier). 

La variable nom contient le nom d'une personne. Le booleen retourne en 
resultat sera VRAI si et seulement si la personne nommee est presente dans l'arbre 
ou le sous-arbre de racine n. Si le nom est present Lender retourne en resultat 
contiendra l'index du nceud comportant le nom. 

■ PROC parenUp, enf, n) -> (bool, entier). 

Comme pour trouver, mais on cherche un couple parent - enfant de noms p et 
enf. L'entier resultant contiendra, le cas echeant, l'index du parent. 

. PROC cousin_germain(a, b, n) -> (bool, entier, entier). 

Recherche de cousins germains de noms a et b. 

. PROC a!eul(a, enf, n) -> (bool, entier, entier). 

Recherche des liens entre l'aieul de nom a et son descendant enf. On imprimera 
(dans le bon ordre) les noms des personnes dans la descendance (a est parent de b, b 

est parent de c i est parent de enf), ou un message approprie si enf n'est pas un 

descendant de a. 

2. Par la suite, on peut se poser les memes questions avec des families plus 
naturelles, ou chaque enfant a un pere et une mere, avec des mariages stables et 
monogames. On garde le controle des naissances (pas plus de deux enfants par 
famille). On engendre un treillis. 

On batira un systeme conversationnel qui accepte des informations des types 
suivants : 

■ a se marie avec b, 

■ naissance de F enfant e du couple m, f. 

L'arbre familial est un fichier sur disquette. Les procedures de recherche peuvent 
prendre la forme de questions (interface ?) comme : 
- Qui est le pere de x ? 
• Qui est la grand-mere maternelle de y ? 

■ Des questions de l'exercice 1. 

3. Dans un arbre binaire complet a n niveaux, combien y a-t-il de noeuds ? 
Combien de feuilles ? 

4. Ecriture du programme de parcours d'un arbre binaire en largeur d'abord. 
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Chapitre 5 

Recurrence et recursivite 



Que ce soit pour la resolution de problemes, pour la specification 
d'algorithmes, voire meme pour leur mise en ceuvre, la recurrence ("induction") est 
un outil mathematique essentiel. Tout informaticien se doit de la maitriser. 
Heureusement, la recurrence est une technique simple, en ddpit de sa puissance. 

Dans les chapitres precedents, nous avons deja vu des programmes bases sur la 
recurrence. II est d'ailleurs difficile d'imaginer des programmes qui n'y font pas appel, 
car la simple boucle en est une forme. L'informaticien parle souvent de reduire un 
probleme en n vers un probleme en n-1. Dans d'autres cas (diviser pour regner ou 
parcours de graphes, par exemple), on parle de la reduction d'un probleme en n vers 
deux problemes en n/2. II s'agit maintenant de rendre le lecteur plus conscient du 
mecanisme et de sa puissance. 

Au niveau de la programmation, nous avons deja vu que la recurrence peut se 
traduire par la re'cursivite' ("recursion") ou par des iterations, avec ou sans pile. Le 
choix de methode depend de la complexite du probleme considere, de la puissance du 
langage de programmation utilise et de l'efficacite requise. 

Comme d'habitude, nous abordons les techniques a travers des exemples. Ceux- 
ci sont souvent bien connus, car c'est un sujet sur lequel beaucoup de collegues ont 

deja travaille. 

5.1. L'exemple type - les tours d'Hanoi 

Si ce probleme est utilise comme exemple dans tous les cours de ce type, ce 
n'est pas un hasard. II est meme arrive au niveau des journaux de vulgarisation 
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[Bezert 1985]. La source la plus ancienne que nous connaissons est de E. Lucas, 
sous le pseudonyme de M. Claus, en 1883. Les tours d'Hanoi donnent lieu a une 
solution elegante par la recurrence, une analyse iterative demandant plus de travail. 

Le titre du probleme vient de l'histoire racontee habituellement, qui est celle de 
moines bouddhistes en Asie du sud-est qui egrenent le temps en transferant des 
disques, tous de tallies differentes, sur un jeu de trois piquets (voir figure 5.1). Dans 
cette figure il n'y a que cinq disques, mais la tradition veut que les moines jouent 
avec 64. L'histoire est une invention du dix-neuvibme siecle, Lucas la placant a 
BdnareS (en Inde). On ne sait pas comment elle s'est resituee a Hanoi . . . 



YSSSSSSSSSSM WSSSSSSSSSSSM %f/S/SSSSSSSA 



B 



Figure 5.1. Position de depart des tours d'Hanoi 



Le jeu consiste a transferer la pile de disques du piquet A vers le piquet B, 
en utilisant C comme piquet de manoeuvre, tout en respectant les deux regies 
suivantes : 

- un disque ne peut pas etre pose sur plus petit que lui, 
• on ne deplace qu'un disque a la fois. 

La solution la plus simple vient de la reponse a la question suivante : "si je 
savais transferer n-1 disques, saurais-je transferer n ?". La reponse est oui, car 
autrement nous n'aurions pas pose la question. Ainsi, pour transferer n disques de A 
a B, on commence par le transfert de n-1 disques de A a C, suivi du deplacement du 
dernier disque (le plus gros) de A a B, suivi du transfert des n-1 disques de C a B. 
Cela donne la procedure du programme 5.1. 
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PROCEDURE hanoi(n, a, b); 
SI n=l 

ALORS d6placement(a, b); 
SINON hanoi(n-1, a, 6-a-b); 
deplacement(a, b); 
hanoi(n-1, 6-a-b, b) 
FINSI 
FINPROC 

Programme 5.1. Tows d'Hanoi, version de base 

Dans cette procedure, nous avons suppose que les piquets sont numerates 1.2, 
3. Deplacer un disque se fait par la procedure deplacement. Pour n=l le deplacement 
est imm£diat, autrement on applique l'algorithme decrit ci-dessus, 6-a-b etant le 
numero du troisieme piquet (6=1+2+3, done en soustrayant deux des trois 
possibility de 6 on obtient la troisieme). 

La figure 5.2 montre les appels en cascade de la procedure hanoi (notee h) pour 
le transfert de trois disques. Le deplacement d'un disque est note d. 



h(3,1,2) 




h(2,1,3) d(1 ,2) h(2,3,2) 




h(1 ,1,2) d(1,3) h(1 ,2,3) h(1.3,1) d(3,2) h(1 ,1,2) 



d(1,2) d(2,3) d(3,1) d(1,2) 

Figure 5.2. Arbre d'appels pour trois disques 

L'arbre se lit de la maniere suivante : 

• Pour transferer 3 disques du piquet 1 au piquet 2 (racine de l'arbre), on 
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transfere deux disques de 1 & 3, un disque de 1 & 2, puis deux disques de 3 h 2 
(successeurs de la racine). 

■ Les successeurs de type h sont decomposes de la meme maniere, en 
appliquant recursivement la procedure. 

• Les feuilles de l'arbre sont des deplacements d'un seul disque. 

En reprenant les feuilles de l'arbre dans l'ordre, de gauche h droite, on deduit les 
deplacements simples, dans l'ordre, pour le cas de trois disques : 

d(1,2) d(1,3) d(2,3) d(1,2) d(3,1) d(3,2) d(1,2) 

Les sceptiques peuvent toujours essayer avec trois pieces de monnaie ! 

5.1.1. Co6t de I'algorithme 

Se pose maintenant la question de savoir combien de deplacements individuels 
sont necessaires pour transferer n disques. En considerant le programme du 
paragraphe precedent, on voit que le transfert de n disques coute deux fois le prix du 
transfert de n-1 disques plus 1 deplacement, par le fait que : 

h(n)= h ( n- 1) + 1 t h(n-l) 

On obtient les equations de r6currence suivantes, avec c representant le cout de 
l'operation, mesure en nombre de deplacements individuels : 

<V l 

c n = 2 * c n _i + 1 pourn>1 

Pour voir la solution de ces equations de recmrence, considerons quelques cas : 

n c n 

1 1 

2 3 

3 7 

4 15 



On voit intuitivement que 
c n =2 n -1 

En substituant cette valeur dans l'equation, la solution est confirmee. 
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5.1.2. Une analyse plus poussee 

La solution recursive au problbme des tours d'Hanoi est satisfaisante pour 
l'esprit et marche dans l'ordinateur. Mais elle comporte l'inconvenient de ne pas etre 
tlfcs pratique pour le joueur humain. Mettons nous un instant a la place des moines, 
qui deplacent leurs disques a la main. Quel disque faut-il deplacer a un moment 
donne ? L'algorithme recursif n&essite la mise a jour constante d'une pile pour 
prendre en compte les appels en cours, mais il ne serait pas tr6s pratique pour les 
moines de gerer une telle pile. 

En fait, un peu de reflexion nous permet de faire mieux, au moins pour jouer a 
la main. Pour obtenir des regies simples pour le joueur humain, nous avons besoin 
de quelques th&Mfcmes. 

Theoreme 1. 

Dans une position quelconque, il y a au plus trois deplacements possibles. 
Preuve. 

Supposons qu'aucun piquet n'est vide. Comparons les tallies des trois disques 
en haut des piquets. Le plus gros des trois ne peut pas bouger, car il faudrait le poser 
sur plus petit que lui. Le moyen peut se poser sur le plus gros (une possibilite), 
mais pas sur le plus petit (qui est plus petit que lui). Le petit peut passer sur l'un ou 
F autre des deux autres piquets, ce qui donne trois possibilites en tout. Maintenant, si 
un des piquets est vide, il joue le role du plus gros disque, qui ne peut pas bouger de 
toute facon, et l'argument ne change pas. Si tous les disques sont sur le meme piquet 
(deux piquets vides), c'est-a-dire dans la position de depart, il n'y a que deux 
possibilites : le plus petit peut aller sur l'un ou l'autre des deux piquets restants. 

Theoreme 2. 

Le disque moyen ne doit pas bouger deux fois de suite. 

Preuve. 

Sur les piquets de depart et d'arrivee, les autres disques sont tous plus gros que 
lui. Sur le troisieme piquet, il y a le plus petit. Apres un deplacement, le moyen 
reste done le moyen. Son seul deplacement possible serait de retourner d'ou il vient, 
ce qui serait parfaitement inutile. 
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Theoreme 3. 

Le plus petit ne doit pas bouger deux fois de suite. 

Preuve. 

Ce n'est pas la peine, car il aurait pu y aller en un coup. 
Theoreme 4. 

Tous les emplacements d' index impair sont des deplacements du plus petit 
disque, les deplacements pairs etant du moyen. 

Preuve. 

C'est le petit qui commence dans la position de depart, deplacement 1. Par la 
suite, le moyen et le plus petit alternent, d'apres les th6oremes 1-3. 

Nous avons demontre que pour savoir qui bouge, il suffit de savoir si le dernier 
transfert etait du petit disque ou du moyen. Le prochain deplacement concerne done 
Fautre disque mobile. Comme les deplacements du disque moyen sont imposes (une 
seule possibilite), ils ne posent pas de problbme. Reste a trouver une r6gle pour les 
deplacements du plus petit. 

Theoreme 5. 

Le plus petit disque parcourt les piquets dans un ordre fixe, e'est-a-dire qu'il fait 
des cercles. 

Preuve. 

Par recurrence sur n. Supposons que la propriete soit vraie pour n- 1, e'est-a-dire 
que le petit disque tourne, par exemple, dans 1' ordre a b c a b c . . . On sait que 
(formule F) : 

n>1 => h(n, a, b) = h(n-1, a, c); d(a,b); h(n-1, c, b) 

Comme le d(a,b) ne concerne pas le plus petit disque, si les deux appels de 
h(n-l, . ..) font tourner le petit disque dans le meme sens, alors la propriete est vraie 



98 



Algorithmique et programmation 



pour n (le tour continue). Comme la propri6td est vraie au niveau de trois disques 
(voir les deplacements resultants de la figure 5.2), elle est vraie pour n disques, n>3. 

Cette preuve est juste, mais elle pose sou vent des probltmes aux eleves en 
raison du changement apparent de la direction du cercle entre le niveau n et le niveau 
n-1. C'est un faux probleme, mais il indique une autre propriete interessante. 
Supposons que le transfert complet souhaite est de A a B. Alors, si le nombre total 
de disques est impair, le premier deplacement (necessairement du plus petit) est 
egalement de A a B. Mais, si le nombre total de disques est pair, le premier 
deplacement est de A a C (Fautre possibilite). Avec un nombre impair de disques, le 
cercle est 

abcabc... 

Avec un nombre pair, il est 
acbacb... 

Si n est pair, n-1 est impair, et inversement. Mais, en regardant les parametres 
des appels de n et de n-1 dans le formule F ci-dessus, on se rend compte de 
l'inversion de la direction. 

Nous pouvons done deduire la regie pour les moines, qui doivent transferer, 
rappelons le, 64 disques. Pour effectuer un transfert de A a B, le premier deplacement 
est de A a C (64 est un nombre pair, pour un nombre impair de disques le 
deplacement aurait ete de A a B). Par la suite, il suffit de se rappeler si le coup qui 
vient d'avoir lieu etait un deplacement du plus petit disque, ou non : 

• Si le deplacement precedent etait du plus petit disque, on deplace le moyen 
(une seule possibilite). 

- Si le deplacement precedent etait du disque moyen, on deplace le plus petit 
vers le prochain piquet dans le cycle ACBACB... 



5.2. Des permutations 

La generation des permutations de n objets est un exercice d' algorithmique et de 
programmation qui presente un interet pedagogique certain tout en menant a des 
programmes utiles. La litterature scientifique comporte de nombreux papiers sur le 
sujet. 

Nous abordons ces algorithmes d'un point de vue personnel, en fonction de 
decouvertes aleatoires. Le premier programme date d'une des ecoles de FL. Bauer. 
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5.2.1. Permutations par ^changes de voisins 

La specification de ce probleme nous a ete donnee dans les annees 70, par un 
collegue. Malheureusement, nous ne nous rappelons plus de qui il s'agissait. Si 
jamais il lit ces lignes, ce serait avec plaisir que nous lui attribuerons son bien . . . Le 
probleme est d'engendrer toutes les permutations des n elements d'un vecteur, en 
appliquant les seules regies suivantes. 

• Entre deux permutations successives, la seule difference nSsulte de l'echange de 
deux voisins. Notons qu'en general v[i] a deux voisins, v[i-l] et v[i+l], mais v[l] et 
v[n] n'ont chacun qu'un seul voisin, ce qui veut dire que le vecteur n'est pas circulaire 
et, en particulier, v[l] et v[n] ne sont pas des voisins. 

• En numerotant chacun des objets par l'index de sa position de depart (l..n), 
l'echange qui a lieu entre deux permutations successives concerne Felement le plus 
grand possible, sans reproduire une permutation deja produite. 

Ces deux regies impose un ordre unique de production des permutations, meme 
si cet ordre n'est pas evident a premiere vue. Afin de bien cerner le probleme, 
considerons F ordre impose par cet algorithme pour les permutations de quatre entiers 
(figure 5.3). 



lero 


Permutation 


Pivot 


Monde 




1234 


4 


123 


2 


1243 


4 




3 


1423 


4 




4 


4123 


3 




5 


4132 


4 


132 


6 


1432 


4 




7 


1342 


4 




8 


1324 


3 




9 


3124 


4 


312 


10 


3 142 


4 




11 


3412 


4 




12 


43 12 


2 




13 


4321 


4 


321 


14 


3421 


4 




15 


3241 


4 




16 


3214 


3 




17 


2314 


4 


231 


18 


2341 


4 




19 


243 1 


4 




20 


4231 


3 




21 


42 13 


4 


213 


22 


2413 


4 




23 


2 143 


4 




24 


2 134 







Figure 5.3. Permutations de quatre entiers 
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Sur cette figure, on voit que le 4 "traverse" les autres objets de droite a gauche, 
puis retourne au point de depart par une traversee de gauche a droite. Au moment de 
son changement de direction, le 4 ne peut pas revenir immediatement, car l'echange 
reproduirait la permutation pr6c6dente. Chaque fois que le 4 termine une traversee, et 
done va changer de direction, F algorithme effectue un echange dans le monde des 
trois". Cela veut dire que les trois objets restants menent leur vie au niveau inferieur. 
Le 4 traverse chaque permutation du monde des trois avant que 1' algorithme en 
produise la prochaine. Dans le monde des trois, le trois traverse chaque permutation 
du monde des deux. Ainsi, on justifie 1' algorithme par la formule de recurrence 
suivante : l'objet n traverse complement chaque permutation des n-1 objets 
restants. On engendre ainsi toutes les permutations de n objets. La demibre colonne 
de la figure 5.3 indique les permutations des trois objets qui sont engendrees a chaque 
changement de direction du quatrieme. Le sens du mot pivot deviendra clair par la 
suite. 



52.2. Le programme 

Dans ce programme, nous travaillons sur Fentier i, qui est l'index de l'echange 
considered L'echange i transforme la permutation i en la permutation i+1. La valeur 
de i varie done de 1 a factoriel(n)-l, ou n est le nombre d'objets, l'echange d'index 
factoriel(n) reproduisant la position de depart. Les objets sont represented par les 
entiers (l..n). Le vecteur v contient, a chaque instant, la derniere permutation 
produite. 

Avec les conventions ci-dessus, nous allons considerer les questions suivantes : 
• A l'echange d'index i, quel est l'objet qui doit s'echanger avec un voisin plus 
petit que lui ? On appelle cet objet le pivot. 

■ Dans quelle direction le pivot doit-il bouger ? 
. Ou se trouve le pivot dans v ? 

La figure 5 . 3 indique le pivot dans le cas des permutations de quatre elements. 
Reprenons cette sequence (figure 5.4). 

12345678910 1112 13 14151617 18 19202 122 2324 
4 4 4 3 4 4 4 3 4 4 4 2 4 4 4 3 4 4 4 3 4 4 4 fin 

Figure 5.4. Valeurs de i et du pivot dans permutations(4) 

Dans cette figure, la Premiere ligne indique les index des echanges successifs, la 
deuxieme le pivot pour l'echange. On voit que le quatre prend trois echanges pour 
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traverser la permutation courante du monde des trois, suivi d'un echange dans ce 
monde des trois. On voit que chaque fois que 4 n'est pas un diviseur de i, c'est le 4 
qui est le pivot. Considerons maintenant le monde des trois. Les echanges le 
concernant sont les 4, 8, 12, 16 et 20, tOUS les i multiples de 4. H suffit de diviser 
ces valeurs par 4 pour retrouver la sequence 1, 2, 3, 4, 5. De maniere recurrente, dans 
le monde des trois, le pivot est le 3, sauf si 3 divise le nouvel i (l'ancien i/4). Pour 
le cas ggndral, on obtient le programme 5.2 pour le calcul du pivot. 

DEBUT nn:=n; ni:=i; 

TANTQUE nn DIVISE ni 
FAIRE ni:=ni/nn; nn:=nn-1 
FAIT 

% nn est le pivot pour I'echange i % 

FIN 

Programme 5.2. Calcul du pivot 

Le fait de recopier n et i dans les variables temporaires nn et ni evite de les 
ddtruiie. L'operateur DIVISE donne un resultat VRAI si son deuxieme operande est 
un multiple exact du premier, FAUX autrement. Dans la plupart de langages de 
programmation, il faudrait 6crire un test du type suivant : 

si nn*entier(ni/nn)=ni alors . . . 

Examinons maintenant le probleme de savoir oil en est le pivot dans ses 
traversees du monde inf6rieiir. A la fin de la boucle ci-dessus, nn est le pivot et ni 
l'index de I'echange dans le monde des nn. Considerons les valeurs de p et q, 
CalCulte comme suit : 

p:=entier(ni/nn); q:=reste(ni/nn); 

Notons que q>0, car la division laisse necessairement un reste (condition de 
terminaison de la boucle). On demon tre que p est le nombre de traversees deja 
effectuees par nn, q etant l'index de I'echange dans la traversee actuelle. Cela resulte 
du fait qu'une traversee necessite nn-1 echanges avec nn comme pivot. 

Avec cette information, on peut deduire les faits suivants. 

• Si p est pair, le pivot procede de droite a gauche, et si p est impair, de gauche 
a droite. Effectivement, la premiere traversee, avec p=0, a lieu de droite a gauche, la 
deuxieme, avec p=l, de gauche a droite et ainsi de suite. 

• Le pivot a effectue q-1 echanges dans cette traversee. II a done avance de q-1 
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positions a partir de son point de depart. 

• Le point de depart du pivot est completement a gauche ou completement a 

droite de l'ensemble des elements du monde des nn, qui se trouvent dans des 
positions successives (tout element d'ordre superieur est en fin de traversee). 

Le point de depart naturel dune traversee est la position 1 (traversee de gauche a 
droite) ou nn (de droite a gauche). Mais considerons, dans la figure 5.3, I'echange 
numero 4, qui est le premier du monde des 3. Le pivot, dans ce cas le 3, n'est pas en 
position 3, mais en position 4. La raison est que le 4 a termine une traversee de 
droite a gauche sans revenir. On voit qu'il faut deplacer le point de depart naturel par 
le nombre d' elements d'ordre superieur qui ont accompli un nombre impair de 
traversees. Tout cela nous donne le programme 5.3. 

DEBUT DONNEES n: entier; 

VAR i, ni, nn, p, q, d, pos: entier; 
VAR perm: TABLEAU [ 1 ..n] DE entier; 
POUR i:= 1 JUSQUA n FAIRE perm[i]:=i FAIT; 
POUR i:= 1 JUSQUA fact(n)-1 
FAIRE ni:=i; nn:=n;d:=0; 
TANTQUE nn DIVISE ni 
FAIRE ni:=ni/nn;nn:=nn-1; 

SI impair(ni) ALORS d:=d+l FINSI 

FAIT; 

p:=entier(ni/nn); q:=reste(ni/nn); 
SI pair(p) 

ALORS pos:=nn+d-q+l ; echanger(perm[pos], perm[pos-1]) 
SINON pos:=q+d; echanger(perm[pos], perm[pos+l ]) 
FINSI 

FAIT 

FIN 

Programme 5.3. Permutationspar echanges 

Dans ce programme, d est le deplacement cree par les traversees d'elements 
d'ordre superieur. Apres la division ni/nn, la nouvelle valeur de ni represente le 
nombre de traversees deja effectuees par l'ancien m. Si ce nombre est impair, nn est 
reste a gauche, provoquant une augmentation de la valeur de d. Le calcul de pos, la 
position du pivot, prend en compte les faits deduits ci-dessus. 

Le calcul de pos, avec le cumul de d, peut etre evite en gardant en memoire le 
vecteur inverse du vecteur perm. L'element perm[i] indique l'objet qui se trouve en 
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position i. L'element inv[i] va indiquer la position de l'objet i en perm. C'est done 
une representation duale de la relation biunivoque 

position <-> objet 

On obtient le programme 5.4. 



DEBUT DONNEES n: entier; 

VAR i, ni, nn, p, q, pos, voisin: entier; 
VAR perm: TABLEAU [1 ..n] DE entier; 
POUR i:= 1 JUSQOa n 
FAIRE perm[i]:=i; invp]:=i 
FAIT; 

POUR i:= 1 JUSQUA fact(n)-1 
FAIRE ni:=i;nn:=n;d:=0; 

TANTQUE nn DIVISE ni 

FAIRE ni:=ni/nn;nn:=nn-1 

FAIT; 

p:=entier(ni/nn); q:=reste(ni/nn); pos:=inv(nn); 

SI pair(p) 

ALORS voisin:=perm[pos-1]; perm[pos-1]:=nn; 

inv[nn]:=pos-1 ; inv[voisin]:=inv[voisin]+l 
SINON voisin :=perm[pos+l]; perm[pos+1]:=nn; 

inv[nn]:=pos+1 ; inv[voisin]:=inv[voisin]-1 
FINSI; 

perm[pos]:=voisin 

FAIT 

FIN 

Programme 5.4. Permutations par echanges, version 2 



Dans cette version, pos est la position du pivot, voisin l'objet avec lequel il va 
s'echanger. Les deux vecteurs sont mis a jour de maniere evidente. 

L'idee de garder la representation inverse d'une relation revient SOUvent dans les 
programmes, comme dans les bases de donnees. Considerons le fichier de 
proprietaires de voitures tenu a jour par une prefecture. On peut, par exemple, le trier 
dans l'ordre alphabetique des noms des proprietaires. Dans ce cas, connaitre le 
numero d'immatriculation de DUPONT est une operation rapide (avec un algorithme 
de recherche dichotomique, par exemple). Mais decouvrir le nom du proprietaire du 
vehicule 9999XX44 devient difficile, car il faut inspecter chaque entree dans le 
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fichier. On voit que la representation choisie doit etre une fonction de l'utilisation 
que Ton va en faire. La double representation est souvent necessaire, e'est-a-dire que 
la prefecture garderait deux fichiers, l'un trie en fonction des noms des proprietaires, 
1'autre par des numeros d'immatriculation. 

5.2.3. Relation avec les tours d'Hanoi 

La technique d'utiliser l'index du coup a jouer marche aussi dans les tours 
d'Hanoi. Supposons que les disques sont numerates de 1 (le plus grand) a n (le plus 
petit). Pour savoir quel est l'index du disque qui doit bouger, on peut utiliser un 
programme semblable a celui que nous avons deduit ci-dessus pour les permutations 
(programme 5.5). 

DEBUT nn:=n; ni:=i; 

TANTQUE 2 DIVISE ni 
FAIRE ni:=ni/2;nn:=nn-1 

FAIT; 

% C'est le disque nn qui bouge % 
FIN 

Programme 5.5. Encore les tours d'Hanoi 

Ce calcul depend du fait que le plus petit disque bouge une fois sur deux, que le 
second bouge une fois sur deux des coups restants, et ainsi de suite. 

L'idee de travailler sur les proprietes numeriques de l'index du coup s'av6re 
interessante dans un ensemble d' applications. 

5.2.4. Une variante 

En relaxant la regie concernant l'echange de voisins, on peut produire un 
deuxieme programme, similaire au premier, qui engendre toutes les permutations de 
n objets dans un autre ordre. Nous allons travailler sur la recurrence suivante : le 
pivot n traverse de droite a gauche chaque permutation du monde des n-1. Pour quatre 
objets, les permutations sont engendrees dans l'ordre de la figure 5.5. 
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Figure 5.5. Nouvelles permutations de quatre objets 



Considerons le programme 5.6. 
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DEBUT DONNEES n: entier; 
VAR i, j, nn, p, q: entier; 
perm: TABLEAU [1 ..n] DE entier; 
i:=0; 

TANTQUE kfact(n) 

FAIRE nn:=n; p:=i; 

POUR j:=1 JUSQUA n FAIRE perm[j]:=0 FAIT; 
TANTQUE nn>0 

FAIRE % placer I'objet nn % 

q:=reste(p/nn); p:=entier(p/nn); j:=n; 
TANTQUE q>0 

FAIRE % sauter q positions vides % 
SI perm[j]=0 ALORS q:=q-1 FINSI; 
j:=j-1 

FAIT, 

TANTQUE perm[j]>0 

FAIRE % trouver la prochaine position vide % 
j:=j-1 

FAIT; 

perm[j]:=nn; nn:=nn-1 

FAIT; 

% perm contient la permutation d'index i % 
i:=i+1 

FAIT 

FIN 

Programme 5.6. Traversee unidirectionnelle 



Pour chaque valeur de i, le programme calcule la position de chaque objet nn 
(l<nn<n) dans le vecteur perm. Par la meme logique que celle de 1'algorithme 
precedent, on prend la partie entiere (en p) et le reste (en q) de la division de p (une 
copie de i au depart) par nn. La partie entiere indique combien de traversees 
completes ont ete effectuees par nn. Le reste indique le nombre d'avancees deja 
effectuees par nn dans la traversee actuelle. Les index i vont de 0 ^ fact(n)-l pour 
faciliter les calculs. 

La valeur de q permet de calculer la position de I'objet nn dans perm. L'objet nn 
a deja avance q fois. En commencant a droite, Falgorithme passe sur q positions 
vides dans perm. L'objet nn va occuper la prochaine position vide. La boucle est 
repetee pour chaque objet n, n-1, n-2, .... 1. La nouvelle valeur de p est l'index de la 
permutation dans le monde des nn-1 . 
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Ce nouvel algorithme a un avantage important par rapport au precedent : il 
calcule chaque permutation a partir des valeurs de i et. de n, sans consulter les 
permutations precedentes. U est done capable de calculer la permutation i sans 
engendrer ses predecesseurs. Cette possibilite sert souvent, car de nombreuses 
applications necessitent de tirer au sort une permutation. Apres le tirage au sort d'un 
entier al6atoire entre 1 et factoriel(n), l'algorithme peut engendrer la permutation 
correspondante. 

5.2.5 Une version recursive 

Ce n'est pas par gout de la difficulte que les programmes pour les permutations 
presenter en premier lieu necessitent une certaine reflection, mais justement pour 
apprendre a raisonner. Evidemment, on peut travailler directement, en ecrivant une 
procedure recursive de maniere "naturelle" (programme 5.7). 

VAR pos: TAB [0..max-1]; 
PROC perm(i, n): 

VAR p, q, j: entier; 

SI n=l 

ALORS pos[1]:=1 

SINON p:=entier(i/n); q:=reste(i/n); perm(p, n-l); 

POUR j:=n-1 PAR -1 JUSQUA q+1 FAIRE pos[j]:=pOS[j-1] FATT; 
pos[q]:=n 
FINSI 
FINPROC 

Programme 5.7. Version recursive 



Dans cette procedure, i est l'index de la permutation a engendrer, 
0<i<faCtoriel(n-l). n est le nombre d'elements a placer dans pos[0..n-l]. La 
numerotation a partir de 0 facilite l'utilisation du systeme de partie entiere et reste, 
deja etudie. 

Ainsi, s'il existe plusieurs elements a placer (n>l), on place n-l Clements par 
l'appel recursif, puis on insere le n-ibme a sa place en decalant les q derniers 
elements chacun d'une case vers la droite. Dans le monde des n-l, l'index de la 
permutation concernee est p (entier(i/n)), 

Cet algorithme, appele a repetition pour les valeurs de i de 0 a factoriel(n)-l, 
engendre les permutations dans l'ordre alphabetique. Cet ordre, classique, conespond 
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a la definition suivante : Les elements dans chaque permutation sont considered 
comme des caracteres (numeriques) d'un mot. Ces mots sont organises en ordre 
alphabetique (ordre numerique si Ton considere les Clements comme les chiffres d'un 
nombre). 

5.3. Exercices 

Sur les tours d'Hanoi 

1. En supposant que les moines de Hanoi veuillent transferer le jeu complet de 64 
disques, en prenant une seconde pour le deplacement individuel d'un disque, combien 
leur faudrait-il de temps ? La legende veut que la terminaison du transfert signale la 
fin du monde. 

2. Etablir une courbe de complexity de l'algorithme des tours d'Hanoi en mesurant le 
temps n6cessaire au calcul pour differents nombres de disques. L'exercice sert a faire 
prendre conscience des problemes d'explosion combinatoire. 

3. Avec 64 disques, et en numerotant les disques dj (le plus petit), d 2 , . . . . d$4 0 e 
plus gros), dans quel ordre d 2 parcourt-il les piquets et, en general, quelle est la regie 
pourdj ? 

4. Considerons les index des deplacements simples i j (le premier coup), i 2 , . . . Quel 
est le numero du disque qui bouge au coup ij ? 

Sur les permutations 

5. La possibilite de tirer au sort une permutation est utile dans de nombreuses 
applications. On reprendra la methode de generation des permutations par echange de 
voisins (§5.2.2), en la modifiant afin de pouvoir calculer la permutation numero i 
sans en engendrer les i-1 precedentes. 

6. Ecrire un programme qui engendre les permutations dans l'ordre alphabetique (la 
valeur numerique de Tender produit en COnCatenant, dans l'ordre, les objets de la 
permutation i, est plus petite que celle de la permutation i+1) sans utiliser une 
procedure recursive. Le programme possedera la propriete de la question precedente 
(engendrer une permutation sans engendrer ses predecesseurs). 

7. On considere les nombres pouvant etre formes a partir des chiffres 1, 2, .... 9, 
chaque nombre comportant une et une seule occurrence de chaque chiffre. En 
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organisant ces nombres dans l'ordre numerique, quelle est la valeur du lOO.OOOeme 

nombre ? (Concours de jeux mafhematiques de France, 1986). 

8. Considerons l'ensemble de valeurs duales des permutations du §5.2.2. (vecteur 
inv). Nous l'avons utilise pour faciliter la recherche de la position courante de chaque 
objet. Est ce que cet ensemble forme un nouvel ordonnancement des permutations de 
n objets ? 

9. Un theoreme bien connu des mathematiciens est que toute permutation de n objets 
est un produit unique de cycles. On demontrera le theoreme en l'accompagnant d'un 
programme qui identifie les cycles concerned pour une permutation donnee. 

10. Le code de Gray est bien connu dans le monde de la transmission d' informations. 
II est deTini de la maniere suivante : 

On considere l'ensemble de caracteres formes d'un nombre fixe de bits. Comme 

d'habitude, n bits permettent de disposer de 2" caracti?res. Dans le code de Gray, deux 
caracteres qui se suivent ne different que par la valeur d'un seul bit. 

On ecrira un programme qui engendre les caracteres successifs du code de Gray sur n 
bits, sachant que le bit qui change de valeur a chaque fois est celui du poids le plus 
faible possible. Par exemple, pour n=3, on aura l'ordre suivant : 

000 
001 

on 

010 
110 

111 

101 
100 
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Chapitre 6 

La marche arriere 



Dans ce chapitre nous abordons une nouvelle categorie d'algorithmes : ceux qui 
mettent en jeu la marche arriere ("backtracking"). II s'agit de resoudre des problemes 
non deterministes. Un problbme qui n'est pas deterministe ne peut pas etre resolu de 
maniere directe. Considerons la situation de quelqu'un qui arrive a un carrefour dans 
un labyrinthe. II ne sait pas si la bonne decision est de tourner a gauche ou a droite. 
La solution est d'essayer une des possibilites ouvertes, en la poursuivant aussi 
longtemps que possible. Si la decision mene a un blocage, il faut revenir sur ses pas 
afin d'essayer une autre possibility. C'est cette idee de revenir sur ses pas qui donne 
lieu a l'expression "marche arriere". 



6.1. La souris et le fromage 

Nous considerons comme application de la marche arriere un probleme de 
labyrinthe artificiel. Sur un echiquier classique (8 cases par 8), il y a un morceau de 
fromage que veut manger une souris. Les coordonnees des cases qu'ils occupent sont 
[xf, yf] et [xs, ys]. Le probleme est d'etablir un chemin que peut prendre la souris 
pour aller jusqu'au fromage, sachant qu'elle ne peut avancer que d'une case a la fois, 
soit horizontalement, soit verticalement, et qu'il existe des barrieres infranchissables 
entre certaines cases voisines. Le probleme est ancien, mais la premiere utilisation 
ayant un rapport avec l'informatique que nous connaissons est [Zemanek 1971]. H a 
ete utilise par un des auteurs comme exercice de programmation depuis une ecole 
d'ete" en 1972 et a servi comme support dans un cours d' intelligence artificielle 
[Griffiths 1986, 19871. 



Algortthmique et programmation 

Nous allons done ecrire un programme mettant en ceuvre la marche arrive pour 
la recherche d'un chemin, sans essay er de trouver le chemin le plus court. Le 
programme presente les caracteristiques suivantes : 

- A chaque tour de la boucle principale, il considere, s'il en existe, un pas en 
avant, sinon il en fait un en arriere. 

• Le pas en avant est dans la premiere direction disponible (qui n'a pas encore 
ete essayee) a partir de la case courante. 

« L'ordre de disponibilite est : nord, est, sud, ouest, qui sont representes par les 
entiers 1, 2, 3, 4. 

. En arrivant sur une case, le programme la marque. Les marques ne sont 
jamais enlevees. La souris ne repasse jamais en marche avant dans une case deja 
marquee, car l'algorithme essaie toutes les possibility a partir de chaque case visitee, 
et il ne faut pas qu'il boucle. 

• Faire marche arriere consiste en un retour vers le predecesseur de la case 
courante, suivi d'un essai dans la prochaine direction disponible a partir de cette case 
(s'il en existe). 

- La marche aniere resulte du constat que Ton a essaye les quatre directions a 
partir d'une case sans trouver un chemin. 

■ Pour faire un pas en arriere, il faut savoir d'ou on est venu. Ainsi, a chaque 
decision, la direction prise est sauvegardee sur une pile. A chaque retour en arriere, 
on revient sur la decision la plus recente, qui est en tete de la pile. Si la situation est 
toujours bloquee, on prends encore le pr6decesseur, et ainsi de suite. 

Le programme 6.1 est donne en entier pour faciliter la presentation, les 
commentaires etant regroupes a la fin. 
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1 DEBUT DONNEES barriere: TAB [1..8, 1..8, 1-4] DE tool; 



2 xs, ys, xf, yf: ©ntisrj 

3 VAR marque: TAB [1..8, 1..8] DE bool INIT FAUX; 

4 pile: TAB [1..64] DE entier; 

5 x, y, xn, yn, dir, pi: entier; 

6 x:=xs; y:=ys; marque[x, y]:=VRAI; dir:=0; pl:=1 ; 

7 TANTQUE NON (x=xf ET y=yf) ET p|>0 

8 FAIRE SI dir<4 

9 ALORS dir:= dir+l ; 

10 SI NON barriere[x, y, dir] 

11 ALORS xn:= CAS dir DANS (x, x+1, x, x-l) FINCAS; 

1 2 yn:= CAS dir DANS (y+1 , y, y-l, y) FINCAS; 

13 SI NON marque[xn, yn] 

14 ALORS x:=xn; y:=yn; marque[x, y]:=VRAI; 

1 5 pile[pl]:=dir; pl:=pl+l ; dir:=0 
FINSI 

FINSI 

1 6 SINON pl:=pl-1 ; 

1 7 SI pl>0 

18 ALORS dir:=pile[pl]; 

19 X> CAS dir DANS (x, x-l, x, x+1) FINCAS; 

20 y:= CAS dir DANS (y-l, y, y+1 , y) FINCAS 
FINSI 

FINSI 

FAIT 



FIN 

Programme 6.1. La souris et lefromage 

Les num6ros de ligne servent aux commentaires. 

1. barriere[i, j, dir] est VRAI si et seulement s'il y a une barriere entre la case [i, j] et 
sa voisine en direction dir. Le tableau est initialise ailleurs. Le bord de l'echiquier est 
entoure de barrieres pour eviter de faire reference a des cases inexistantes. 

2. [xs, ys] est la position de depart de la souris, [xf, yfj celle du fromage. 

3. marquefi, j] est VRAI si et seulement si la souris a deja visite la case [i, j]. Ce 
tableau est initialise a FAUX. 

4. Une pile de 64 Clements suffit. 

5. [x, y] est la position courante de la souris. [xn, yn] est une position visee par la 
souris. dir represente le nombre de directions que la souris a deja essayees a partir de 
la case [x, y]. 
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6. L'initialisation met [x, y] & la position de depart [xs, ys], qui est marquee. Aucune 
direction n'a ete essay^e. La pile est vide. 

7. On boucle tant que le fromage n'est pas atteint et qu'il existe des chemins k 

essayer. 

8. Existe-t-il encore une direction a essayer ? 

9. Alors on prend la prochaine. 

10. Si la route n'est pas barree, on calcule les coordonnees [xn, yn] de la case voisine 
en direction dir. 

11,12. dir=l veut dire vers le nord, c'est-a-dire augmenter y de 1, et ainsi de suite. 

13. Dans 1' absence de marque sur [xn, yn], on peut y aller. Si la case est marqu6e, on 
essaiera la prochaine valeur de dir au prochain tour de la boucle. 

14. Le mouvement a lieu, avec marquage. 

15. La decision est enregistree sur la pile. Pour la nouvelle case le nombre de 
directions d6}& essayees est 0. 

16. Comme il n'existe plus de directions & essayer, il faut faire marche arriere. 

17. Si la pile n'est pas vide, on prend la direction prise pour arriver sur la case 
courante, qui se trouve en haut de la pile (18). pl=0 arrive quand tous les chemins 
possibles ont ete essayes k partir de [xs, ys] et le fromage n'a pas ete atteint. II n'y a 
done pas de chemin possible. 

19,20. On calcule, en inversant les lignes 11 et 12, les coordonnees du pred6cesseur, 
et le programme fait un nouveau tour de la boucle. 

A la fin du programme, si un chemin existe, la pile contient toutes les 
decisions prises sur lesquelles l'algorithme n'est pas revenu. Ainsi la pile contient 
une solution du probleme. II n'est pas interessant de programmer ici une boucle 
d'impression, mais le principe de la disponibilite des resultats acquis est valable pour 
les algorithmes de marche arriere en general. Pendant l'execution du programme, la 
pile contient & chaque moment les decisions en cours d'essai. 

6.1.1. Version recursive 

Dans le cas d'une recherche complete sur toutes les possibilites, la marche 
arriere est l'equivalent d'un parcours d'arbre ou de graphe (voir §6.1,2). Dans le 
probleme considere ici, nous sommes en presence d'un graphe, car il serait possible 
de revenir sur la case que l'algorithme vient de quitter, et Ton pourrait boucler. Cela 
est une autre explication du besoin de marquage, pour empecher le programme de 
suivre les boucles. Comme dans tous les algorithmes de ce type, on peut ecrire le 
programme dans la forme d'une procedure recursive (programme 6.2). 
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DEBUT DONNEES xs, ys, xf, yf: entier; 
barre: TAB [1..8, 1..8, 1..4] DE bool; 
VAR marque: TAB [1..8, 1..8] DE bool INIT FAUX; 

existe_chemin: bool; 
PROC joindre(x, y) -> bool; 
VAR trouve: bool; 

dir: entier; 
trouve:= x=xf ET y=yf; 
SI NON marqu6[x, y] 
ALORS marqu6[x, y]:=VRAI;dir:=0; 
TANTQUE NON trouve ET dir<4 
FAIRE dir:=dir+l ; 
SI NON barr<S[x, y, dir] 

ALORS trouve:=joindre(CAS dir DANS x, x+1, x, x-l FINCAS, 
CAS dir DANS y+1 , y, y-l, y FINCAS) 

FINSI 

FAIT; 

St trouve ALORS imprimer(x, y) FINSI 
FINSI; 

RETOURNER trouve 
FINPROC; 

existe_chemin:=joindre(xs, ys) 

FIN 

P rogramme 6.2. Forme recursive 



Le lecteur ayant tout suivi jusqu'ici n'aura pas de probleme pour comprendre ce 
programme. 



6.1.2. Marche arriere, arbres et graphes 

La marche arrfere est un moyen d' explorer systematiquement un arbre de 
possibilites. Evidemment, F exploration systematique ne termine que si le nombre de 
possibility est fini. II faut aussi noter l'existence de l'explosion combinatoire, ou le 
nombre de possibilites, bien que fini, est tellement grand que l'exploration est en 
dehors des capacites des ordinateurs. 

La figure 6.1 montre quelques branches de F arbre de possibility pour la souris. 
Commencant sur la case [xs, ys], elle a quatre directions disponibles (en ignorant des 
barrieres et les limites de Fechiquier). Considerons la direction "est", menant a la 
case [xs+1, ys]. A partir de cette nouvelle case, quatre directions sont encore 
imaginables, dressees sur la figure. 
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[xs, ys] 




[xs, ys+1] [xs+1,ys] [xs,ys-1] [xs-1,ys] 




[xs+1,ys+1] [xs+2, ys] [xs+1 , ys-1] [xs, ys] 
Figure 6.1. Arbre de possibility pour la souris 

Or, une de ces directions mbne la souris de retour a sa case de depart. Si Ton 
continuait a dresser un arbre, il serait de taille infinie, car on recommence 
indefiniment avec la meme case. Pour representer le labyrinthe par une structure 
formelle, il faut passer a un graphe, avec un noeud par case. Les arcs sont toujours 
bidirectionnels, reliant des cases voisines qui ne sont pas separees par des barrieres. 
Notons que ce n'est pas seulement le fait de pouvoir retourner directement vers la 
case precedente qui impose un graphe; on peut aussi arriver, par exemple, a la case 
[xs+1, ys+1] en deux pas de deux facons differentes (par [xs, ys+1] et par [xs+1, 

ys]). 

Le marquage des cases du labyrinthe correspond done exactement au marquage 
des nceuds dans le parcours d'un graphe. Dans le cas du labyrinthe, le graphe est le 
graphe de connectivite des cases, representees par des nceuds. Le programme 6.2 a la 
forme d'un parcours de graphe n-aire. 

6.2. Les huit reines 

II s'agit encore d'un exercice archi-classique, mais sa presence ici est toujours 
justifiee par des raisons pedagogiques. II faut placer huit reines sur un echiquier de 
telle facon qu'aucune reine ne soit en prise par rapport a une autre. Deux reines sont 
mutuellement en prise si elles se situent sur la meme ligne, la meme colonne ou la 
meme diagonale. Le probteme est ancien [Bezzel 18481 [Gauss 1850], mais il a ete 
remis au gout du jour avec l'arrivee des ordinateurs. De nombreux articles modemes 
lui sont consacres, jusqu'a dans des revues populaires [Paclet 19841. II y a 92 
solutions, dont evidemment un grand nombre dues a la symetrie. 
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Dans le programme 6.3, np est un entier qui indique le nombre de reines 
actuellement placees. En marche avant, le programme considere la prochaine case [x, 
y]. Le test x>8 au debut de la boucle sert a passer a la ligne suivante (y:=y+l) en cas 
de besoin. Si la case [x, y] n'est pas couverte par une dame deja placee, on peut y 
placer une nouvelle dame. Si la case n'est pas libre, la prochaine case sera consideree 
au prochain tour de la boucle (x:=x+l a la derniere ligne). 



DEBUT VAR np: entier INIT 0; % Le nombre de reines placees % 
x, y: entier; % coordonnees d'une case a essayer % 
x:=1;y:=1; 
TANTQUE np c 8 

FAIRE SI x > 8 ALORS x:=1 ; y:=y+1 FINSI; 
Sly<9 

ALORS SI libre(x, y) 
ALORS np:=np+l ; empiler(x, y); 
FINSI 

SINON d<§piler(x, y); np:=np-1 

FINSI; 

x:=x+1 

FAIT 
FIN 

Programme 6.3. Les huit reines 



L'indication de la necessite de faire marche arriere est qu'il ne reste plus de case 
a essayer (y>8). Dans ce cas, on depile x et y pour revenir dans l'etat precedent, en 
diminuant le compteur np. Le prochain tour de la boucle considerera la prochaine 
case. Le processus s'arrete quand np arrive a 8, c'est-Mire quand solution est trouvee. 
Comme nous savons qu'une telle solution existe, il n'est pas necessaire de considerer 
le cas ou la marche arriere ne peut pas avoir lieu. 

Dans ce programme, la pile contient les coordonnees des reines placees. 
Comme il y a deux valeurs a empiler (x et y), il peut devenir interessant de definir 
deux piles, px et py. Dans ce cas, notons que np peut servir pour indiquer la hauteur 
de chaque pile. Par ailleurs, il nous faut remplacer l'appel de la procedure libre par 
une boucle qui decide si la case (x,y) est couverte par une reine deja placee (meme 
ligne, colonne ou diagonale). On obtient le programme 6.4. 
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DEBUT VAR px, py: TAB [1..8] DE entier; 
np: entier INIT 0; 
x, y: entier INIT 1; 
libre: bool; 
i: entier; 
TANTQUE np < 8 

FAIRE SI x > 8 ALORS x:=1 ; y:=y+1 FINSI; 
Sly<9 

ALORS libre:=VRAI; i:=0; 
TANTQUE libre ET knp 
FAIRE i:=i+l ; 

libre:=NON(x=px[i] OU y=py[i] OU abs(x-px[i])=abs(y-py[i])) 
FAIT; 
SI libre 

ALORS np:=np+l ; px[np]:=x; py[np]:=y 
FINSI 

SINON x:=px[np]; y:=py[np]; np:=np-1 
FINSI; 
x:=x+1 
FAIT 



Programme 6.4. Les huit reines, version 2 



6.2.1. Une version amtliorie 

La solution ci-dessus marche, mais on peut l'ameliorer. Profitons du fait que 

nous Savons qu'il y a une et une seule reine par colonne. II n'est plus necessaire de 
gander une pile des x (programme 6.5). 
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DEBUT VAR py: TAB [1..8] DE entier; 
x, y: entier INIT 1; 
libre: bool; 
i: entier; 
TANTQUE x < 9 
FAIRE SI y < 9 
ALORS libre:=VRAI; i:=l ; 
TANTQUE libre ET j<x 

FAIRE libre:= NON( y=py[i] OU abs(x-i)=abs(y-py[i])) 

i:=i+1 
FAIT; 

SI libre ALORS py[x]:=y; x:=x+1 ; y:=0 FINSI 
SINON x:=x-l ; y:=py[x]; 
FINSI; 
y:=y+l 
FAIT 
FIN 

Programme 6.5. Suppression d'une pile 

Dans cette version du programme, la variable np disparait (sa valeur est la 
meme que celle de x-1). La variable x est la coordonnee en x de la nouvelle reine a 
placer, y etant l'autre coordonnee de la prochaine case a essayer. La marche arriere 
devient necessaire si y arrive a 9, c'est-a-dire s'il n'y a plus de cases a essayer pour la 
valeur courante de x. Le test libre ne comporte plus de test sur x, car les x sont tous 
differents par construction, x servant egalement de pointeur dans la pile des y (il 
indique la premiere place libre dans cette pile). 

6.2.2. Une deuxieme approche 

La dernfere version ci-dessus ressemble un peu a un compteur octal, qui egrene 
tous les nombres en octal de 0 a 77777777. Pour faciliter la programmation, nous 
avons numerate les cases de l'echiquier de 0 a 7 au lieu de 1 a 8. Considerons le 
programme 6.6. 
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DEBUT VAR comp: TAB [0..7] DE entier INIT 0; 
VAR couvert: bool INIT VRAI; 

i, j: entier; 
TANTQUE couvert 
FAIRE % avancer le compteur % 

i:=7; 

TANTQUE comp[i]=7 FAIRE comp[i]:=0; i:=i-1 FAIT; 
comp[i]:=comp[i]+1 ; 

% tester la position creee % 
i:=1 ; couvert:=FAUX; 
TANTQUE NON couvert ET i<8 
FAIRE j:=0; 

TANTQUE NON couvert ET j<i 
FAIRE couvert:=comp[i]=comp[j] 

OU abs(i-comp[i])=abs(j-comp[j]); 

j:=j+l 

FAIT; 
i:=i+l 

FAIT 

FAIT 

FIN 

Programme 6.6. Avec un compteur octal 

Dans ce programme, le compteur (comp) est initialise & la valeur zero. Le 
booleen couvert indique que deux reines "se voient". Son initialisation a VRAI est 
une evidence : dans la position de depart, toutes les reines "se voient". La boucle 
externe tourne tant que la demiere position examinee n'est pas la bonne (aucune paire 
de reines en prise mutuelle). 

Le contenu de la boucle est en deux parties : l'avancement du compteur et 
l'examen de la nouvelle position creee. Les 7 £ droite du compteur deviennent des 0, 
et le prochain chiffre a gauche augmente de 1. Puis, pour chaque reine i (1...7), on 
teste si la reine j (0...i-l) est compatible. Les tests s'arretent a la premiere 
incompatibilite. 

Le nombre de positions & essayer avant de decouvrir toutes les solutions est de 
8 8 , c'est-a-dire 2^, ou 16 Mega dans le sens de l'informatique (1 Mega = 
1024*1024). Considerer un peu plus de 16 millions de cas prend un certain temps. 
Notons que la marche arriere diminue considerablement ce chiffre. En effet, les 
premieres positions essayees dans 1'algorithme du compteur sont les suivantes : 
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00000000 

00000001 

00000002 

et ainsi de suite. Avec la marche arriere, on decouvre immediate me nt que tout 
nombre commencant par 0 0... est a eliminer. La marche arriere donne un 
algorithme qui considere, dans l'ordre arithmetique : 

00000000 
01000000 
02000000 
02100000 

et ainsi de suite. Compter le nombre de cas essayes est laisse" comme exercice. 

Cette ligne de reflexion pourrait nous amener a considerer, a la place du 
compteur octal, les permutations sur 8 objets. On sait en effet que tOUS les x, 
comme tous les y, doivent etre differents. Le nombre de cas a considerer au total est 
de factoriel(8), c'est-a-dire 40320 (& comparer avec les plus de 16 millions pour le 
compteur octal). En faisant de la marche arri&re sur les permutations, il deviendrait 
possible de considerer encore moins de cas. 

Le programme de depart denombre les permutations successives, en utilisant 
un des programmes du chapitre 5 (programme 6.7). 

DEBUT VAR i: entier INIT 1, trouve: bool INIT FAUX; 
TANTQUE NON trouve 
FAIRE trouv6:=libre(perm(i)) 
i:=i+1 

FAIT 
FIN 

Programme 6.7. Avec despermutations 



En y ajoutant de la marche arriere, on pourrait arriver au programme 6.8. 
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DEBUT VAR perm: TAB [1..8] DE entier INIT 0; 
disp: TAB [1..8] DE tool INIT VRAI; 
n, i: entier INIT 0; 
TANTQUE n<8 

FAIRE i:=i+l ; 

TANTQUE i<9 ETPUIS NON disp[i] FAIRE i:=i+l FAIT; 
Sli<9 

ALORS couvert:=FAUX; j:=0; 

TANTQUE ]<n ETPUIS NON couvert 

FAIRE j:=j+1 ; couvert:= n+1-j = abs(i-perm[j]) 

FAIT; 

SI NON couvert 

ALORS n:=n+l; perm[n]:=i; disp[i]:=FAUX; i:=0 
FINSI 

SINON i:=perm[n]; disp[i]:=VRAI; perm[n]:=0; n:=n-1 
FINSI 

FAIT 

FIN 

Programme 6.8. Avec permutations et marche arriere 



Dans ce programme, les permutations & essayer sont engendrees en perm. 
L'616ment disp[i] du tableau bool^en disp indique si l'objet i est disponible, c'est-k- 
dire s'il ne figure pas deja dans la permutation : disp[i] est VRAI si i ne figure pas 
dans perm[l..n]. n indique le nombre de reines deja plac6es. A chaque moment, la 
valeur de i est Findex du dernier objet essaye comme candidat a la place perm[n+l], 
qui est la prochaine reine a placer. 



6.3. Exercices 

Sur la souris et le fromage 

1. Ecrire le programme qui ajoute l'heuristique, c'est-a-dire qui prend comme 
direction pr6f6r6e celle du fromage. 

2. Ecrire le programme de numerotation successive qui trouve un chemin de 
longueur minimale. 

3. II existe combien de chemins differents allant de la case [1, 11 & la case [8, 8] sur 
un echiquier vide ? Une meme case ne peut pas figurer deux fois sur un meme 
chemin. On 6crira un programme qui compte le nombre de chemins. 
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Sur les huit reines 

4. On modifiera le programme donne au §6.2.1 afin de trouver toutes les solutions. 

5. Ecrire la version recursive du programme du §6.2.1. 

6. Ecrire un programme qui compte le nombre de cas essayes en marche arrifere sur le 
compteur octal (voir §6.2.2). 

7. Faire la meme chose pour le programme base sur les permutations, 

D'autres applications de la marche arriere 

8. Extrait partiel du concours de jeux mathematiques, 1987. 

On considere un triangle de cinq lignes : 
X 

x x 

XXX 
X X X X 
X X X X X 

Dans chacune des 15 positions, il faut placer un des 15 premiers entiers en 
respectant la condition suivante : Tender dans chaque position, sauf dans la demi&re 
ligne, a comme valeur la difference entre les valeurs de ses successeurs immediats 
(on COnsid&re le triangle comme un treillis binaire). Un debut de solution pomrait 
etre : 

7 

103 
11 14 



Chaque entier ne figure qu'une seule fois dans le triangle. On ecrira un 
programme utilisant la marche arriere pour trouver une solution au probleme. 
Combien y a-t-il de solutions en tout ? 



123 



Chapitre 7 



Transformation de programmes 



La transformation de programmes, par des techniques automatiques ou 
manuelles, devrait interesser l'informaticien professionnel pour plusieurs raisons : 

• La transformation automatique de specifications en programmes exeeutables 
peut diminuer le cotit de la production de programmes, tout en ameliorant leur 

fiabilit6. 

• On peut d£couvrir des versions plus efficaces de certains programmes. 

- En restant proche de la specification d'un programme, on en facilite la 
portability. 

• Le programmeur am61iore sa maitrise des outils, ce qui nous interesse pour la 
pedagogie. 

Le besoin d'industrialiser la production des logiciels a provoque l'introduction 
de nombreux mots-cles nouveaux. Ces neologismes, tels que programmation 
structured, genie logiciel, . . . . ne represented pas seulement un phenomene de mode, 
mais recouvrent aussi une realite technique. Afin de s'assurer de la qualite des 
produits industriels, on a souvent recours a la specification formelle. Celle-ci sert a 
la demonstration eventuelle d'un programme, a sa validation, voire me me a sa 
production. C'est pour ce dernier aspect que les techniques de transformation 
deviennent utiles. Une specification est une formulation mathematique, surtout 
algebrique, d'un probleme. Dans les methodes de demonstration de programmes, on 
dispose d'une specification et d'un programme, puis on demontre que le programme 
est une mise en ceuvre de la specification. Plutot que demontrer F equivalence de ces 
deux textes, on est arrive naturellement a se poser la question de savoir s'il est 
possible de deriver le programme a partir de la specification. 
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Dans certains cas, la reformulation d'un programme peut s'averer plus efficace 
que la version d'origine. Nous avons vu l'effet de certaines transformations dans le 
cadre de programmes r6cursifs dans un chapitre precedent. 

La transformation d'un programme n'est possible que si le programme est ecrit 
dans une forme propre. L'utilisation d'astuces dependant d'un compilateur particulier 
est a proscrire. Notons que cela revient a dire les memes choses que Ton retrouve 
dans les ouvrages sur la portabilite des logiciels [Brown 19771. 

Enfin, meme si les techniques de transformation n'avaient aucune application 
pratique directe, nous les enseignerions quand meme a nos etudiants. Effectivement, 
en examinant differentes versions d'un meme programme, ils sdparent mieux ce qui 
est intrinseque au probleme de ce qui depend du langage, du compilateur, voire de 
l'ordinateur particulier qu'ils utilisent a un instant-donne. Les paragraphes suivants 
donnent quelques exemples de transformations. Pour aller plus loin, on peut 
consulter utilement des ouvrages tels que [Burstall 1975], [Darlington 1973]. 

7.1. Revenons sur le mode d'un vecteur 

Ce problfeme a forme le sujet du premier vrai programme de cet ouvrage ($2.1). 
II a ete le sujet de plusieurs etudes [Griffiths 19761, [Arsac 1984], car on peut 
l'ameliorer tres simplement, cependant l'amelioration reste difficile a trouver. Nous 
avons decouvert la modification en explorant les techniques de transformation. La 
version de depart de cet algorithme etait celle du programme 7.1. 

DEBUT DONNEES n: entier; 

v: TABLEAU [1 „n] DE entier; 
VAR i, Icour, Imax, indm: entier INIT 1; 
TANTQUE i<n 
FAIRE j:=i+1 ; 

SI v[i]=v[i-1] 

ALORS lcour:=lcour+1; 
SI IcouNmax 

ALORS lmax:=lcour; indm<-i 

FINSI 
SINON lcour:=1 
FINSI 

FAIT 
FIN 

Programme 7.1. Encore le mode d'un vecteur 
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C'est en essayant de rediger la specification formelle de l'algorithme que nous 
avons trouve une amelioration. Essayons de denombrer les cas qui se prdsentent On 
COnsidfcre que le mode d'un vecteur est un doublet : la longueur de la chaine la plus 
longue et la valeur qui parait dans cette chaine. On obtient les cas suivants : 

1. n-l: mode(v[l..n]) = (1, v[1]) 

2. n>1 , v[n]*v[n-1]: mode(v[l ..n]) = mode(v[l ..n-1]) 

3. n>l, v[n]=v[n-1], long(mode(v[1 ..n]))>long(mode(v[1 ..n-1])): 

mode(v[1 „n]) = (long(mode(v[l ..n-1]))+1, v[n]) 

4. n>1, v[n]=v[n-1], long(mode(v[1 ..n]))<long(mode(v[1 ..n-1])): 

mode(v[l..n]) = mode(v[l..n-1]) 

L'amelioration est venu en essayant de formaliser le test qui decide si la chaine 
courante est plus longue que le mode trouve jusqu'alors, c'est-a-dire le test suivant : 

long(mode(v[l ..n]))>long{mode(v[1 ..n-1])) 

Evidemment, on peut laisser le test sous cette forme, mais on cherche 
naturellement un test qui decide si la nouvelle chaine est plus longue que l'ancienne. 
C'est en posant la question de cette maniere que nQUS avons trouve l'idee (pourtant 
simple) de COnsidfrer le test suivant : soit lm la longueur du mode de v[l..n-l]. 
Alors, la chaine dont v[n] fait partie est plus longue que lm si et seulement si 

v[n] = v[n-lm] 



Dans la transformation de la specification ci-dessus, on remplacera le test sur 
les longueurs par ce nouveau test. Par ailleurs, on note que les cas 2 et 4 donne la 
meme solution (le mode est inchange). Mais, si v[n]#v(n-l), il est certain que 
v[n]#v[n-hn]. On en deduit que la comparaison de v[n] avec son predecesseur n'est 
pas necessaire. La specification devient : 

1. n=l: mode(v[l..n])=(l, v[1]) 

3. n>1, v[n]=v[n-long(mode(v[1 ..n-1]))]: 

mode(v[l ..n]) = (long(mode(v[1 ..n-1]))+1, v[n]) 

4. n>1 , v[n]*v[n-long(mode(v[1..n-1]))]: 

mode(v(1 ..n)) = mode(v(1 ..n-1)) 

Le test v[n]=v[n-l] a disparu. En transformant cette specification vers un 
programme iteratif, on obtient le programme 7.2. 
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DEBUT DONNEES n: entier, 

v: TAB [1 ..n] DE entier; 
VAR i, indm, Im: entier INIT 1; 
TANTQUE i<n 
FAIRE i:=i+1 ; 

SI v[i]=v[i-lm] % regard en arriere 

ALORS indm:=i; lm:=lm+1 % nouveau max 

FINSI 

FAIT 

FIN 

Programme 7.2. Le mode ameliore 

On a gagne une variable et un test par element du tableau. Cette optimisation, 
bien que simple, n'est pas evidente au depart 

7.2. La multiplication des gitans 

L'idee de ce programme nous est venu de P.L. Bauer, qui l'a incorpore dans son 
cours de programmation [Bauer 19761. Elle vient de l'observation de la methode de 
multiplication utilisee par des gens sur des marches de l'Europe centrale. Montrons 
cette methode sur un exemple, 25 * 40 : 

25 40 
12 80 

6 160 

3 320 

1 640 

On a divise les nombres de la colonne gauche par deux, en multipliant par deux 
ceux de la colonne de droite en meme temps. L'operation est repetee autant de fois 
que possible, en laissant tomber les moities resultant de la division par deux d'un 
nombre impair. Par la suite, on supprime, dans la colonne de droite, les valeurs 
correspondant £ un nombre pair dans la colonne de gauche. Les nombres restant £ 
droite sont additionnes : 

40 + 320 + 640= 1000 

Ce qui est le r&ultat escompte. 

Les operations ci-dessus ne sont pas le "tour de passe-passe" que Ton pourrait 
imaginer. En fait, 6 * 160 (deuxieme ligne) est strictement egal a 3 * 320. Mais, 
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parce que 3 est un nombre impair, 3 * 320 = 1 * 640 + 1 * 320. Done il faut garder 
les valeurs en face des nombres impairs. En ecrivant ceci dans la forme d'une 
procedure, on obtient le programme 7.3. 

mult(a,b):SI a=0 

ALORS 0 
SINON SI impair(a) 

ALORS b + mult(a-1, b) 

SINON mult(a/2, b*2) 

FINSI 

FINSI 

Programme 7.3. Multiplication avec recursivite 

Pour enlever la recursivite, on introduit une variable res pour le cumul des 
valeurs retenues (programme 7.4). 

multp(a,b,res):TANTQUE a>0 
FAIRE SI impair(a) 

ALORS res:=res+b; a:=a-1 

SINON a:=a/2;b:=b*2 

FINSI 
FAIT 

Programme 74. Suppression de la rtcursivitt 



Comme il s'agit dentiers, representes sous forme binaire dans l'ordinateur, on 
peut simplifier en faisant appel a des decalages (decg pour decaler dune place vers la 
gauche, deed a droite). Par ailleurs, si a est impair, a-1 est pair, done le SINON peut 
disparaitre (programme 7.5). 

multp(a,b,res):TANTQUE a>0 
FAIRE SI impair(a) 
ALORS res:=res+b; a:=a-1 
FINSI; 

decd(a); decg(b) 
FAIT 

Programme 7.5. Version avec decalages 
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Etant donne que a est ddcate & droite, a:=a-l n'est pas ndcessaire, car le 1 que 

Ton enfeve est le dernier bit a droite, supprime de toute facon par le decalage qui va 
suivre. Cela nous permet de deduire ralgorithme de multiplication reellement mis en 
ceuvre par les circuits dans les ordinateurs (figure 7.1). 



100101... 



B 



res 

Figure 7.1. Multiplication de deux entiers 



On decale a vers la droite jusqu'a Tarrivde d'un 1 dans le bit de droite. On decale 
d'autant b vers la gauche. On ajoute la nouvelle valeur de b dans le registre resultat. 
La triple operation est repetee pour chaque bit de a ayant la valeur 1. 

L'introduction d'une variable permettant de supprimer une recursivite sans 
recourir a Line pile est une technique habituelle. On le fait sans reflechir dans 
beaucoup de cas. Considerons le programme bateau de calcul d'un factoriel : 

fact(n): SI n=1 

ALORS 1 

SINON n . fact(n-1) 

FINSI 

Aucun programmeur n'hesite a ecrire la boucle suivante 
fact:=l ; 

TANTQUE n>1 FAIRE fact:=fact * n; n:=n-1 FAIT 

Mais il hesiterait beaucoup plus sur les fonctions suivantes : 
divfact(n): SI n=l ALORS 1 SINON n / divfact(n-l) FINSI 
div2(n): SI n=l ALORS 1 SINON div2(n-1) / n FINSI 
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On utilise des connaissances intuitives des proprietes telles que la 
commutativite, la distributivite, . . . des operateurs. Une etude des conditions 
n&essaires et suffisantes pour permettre des transformations ne peut qu'ameliorer 
notre maitrise de l'outil de programmation. 

Bauer signale que le programme de multiplication peut se generaliser. Le meme 
schema de programme permet la minimalisation du nombre de multiplications 
n&essaires pour le calcul d'une puissance entiere (voir exercice en fin de chapitre). 

7.3. Exercices 

1. Utiliser le programme du 57.2 comme schema pour le calcul de a b . 

2. Enlever la recursivite des procedures divfact et div2 du $7.2. 
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Quelques structures 
de donnees particulieres 



Dans ce chapitre, nous introduisons plusieurs formes nouvelles d'arbres : les 
arbres ordonnes ("ordered trees"), les arbres equilibres ("balanced trees") et les B- 
aibres ("B-trees") 

8.1. Les arbres ordonnes 

Pour qu'un arbre soit ordonne, il faut attacher une valeur a chaque noeud. Dans 
la pratique, comme par exemple dans des systemes de gestion de fichiers ou de bases 
de donnees, cette valeur est souvent une cle, du type decrit pour l'adressage disperse 
(voir chapitre 2) ou d'un style approchant. Ordonner un arbre revient a faire en sorte 
que, si la valeur attachee au noeud n est val[n], chaque noeud dans le sous-arbre dont 
la racine est le successeur gauche de n a une valeur associee vg tel que vg<val[n], les 
noeuds a droite possedant des valeurs superieures. Cette definition implique qu'une 
valeur ne peut figurer qu'une seule fois dans l'arbre. La figure 8.1 montre un exemple 
d' arbre binaire ordonne (voir page suivante). 

L'arbre de la figure aurait pu etre produit par l'algorithme naturel de 
construction si les elements considered s'etaient presented dans l'un des ordres 
suivants : 



7 10 4 2 8 3 1 
7 4 2 110 8 3 
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1 



Figure 8.1. Un arbre binaire ordonne 

Ces ordres sont des ordres partiels de parcours de l'arbre, ou l'ordre partiel doit 
respecter la regie suivante : un noeud devance, dans l'ordonnancement, ses 

successeurs. 

L'algorithme naturel de construction de l'arbre est aussi celui qui sert & 
rechercher la presence d'un objet dans l'arbre (programme 8.1). 



DEBUT DONNEES succg, succd, val: TABLEAU [1 ..taille] DE entier; 
v, racine, pi: entier; 
VAR n: entier; 
n:=racine; 
TANTQUE v*val[n] 
FAIRE SI v<val[n] 
ALORS SI succg[n]>0 
ALORS n:=succg[n] 
SINON succg[n]:=pl; n:=pl;pl:=pl+1; 

succg[n]:=0; succd[n]:=0; val[n]:=v 
FINSI 

SINON SI succd[n]>0 
ALORS n:=succd[n] 
SINON succd[n]:=pl; n:=pl; pl:=pl+1 ; 

succg[n]:=0; succd[n]:=0; val[n]:=v 
FINSI 
FINSI 
FAIT 
FIN 

Programme 8.1. Arbre ordonni 
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L'arbre d6fini par les trois vecteurs succg, succd et val est fourni en entree, 
avec sa racine, la valeur v & trouver ou a placer et 1' index pi de la premiere case libre 
dans chaque tableau. Si v est deja dans l'arbre, n se positionne sur le noeud qui porte 
la valeur v, sinon on cree un nouveau noeud & l'endroit adequat de l'arbre. 

Avec un arbre ordonnd, la recherche d'un objet se fait plus vite qu'avec les 
algorithmes precedents (§4.3.3), qui examinaient tous les nceuds de l'arbre. 
L'ordonnancement permet de suivre toujours le bon chemin, c'est-a-dire celui sur 
lequel est, ou sera, le noeud possedant la valeur recherchee. Le nombre de tests depend 
ainsi de la profondeur de l'arbre. 

8.2. Les arbres equilibres 

L'algorithme du paragraphe precedent ameliore la vitesse de recherche d'un 
objet, mais le nombre de noeuds a examiner n'est pas encore minimise. En effet, 
considerons l'arbre (figure 8.2) qui resulterait d'une arrivee d'objets dans l'ordre 
suivant : 

1 2347810 



1 





Figure 8.2. Arbre ordonne inefficace 
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Les recherches sont minimisees seulement si l'arbre est equilibre, ou du moins 
si la profondeur maximale de l'arbre est minimisee. La figure 8.3 montre l'arbre de 
l'exemple dans sa forme equilibree. 



4 




Figure 8.3. Equilibrage de l'arbre de la figure 8.1 

Soient ng le nombre de noeuds dans le sous-arbre gauche d'un nceud donne, nd 
le nombre de noeuds dans le sous-arbre droit. Alors un arbre est equilibre si et 
seulement si abs(ng-nd)<2. Si, pour tout noeud, ng=nd, l'arbre est strictement 
equilibre, comme celui de la figure 8.3. Cela ne peut arriver que si le nombre de 
noeuds le permet (nombre total de noeuds = 2' ■ 1). On voit que l'algorithme de 
recherche dans un arbre ordonne et equilibre a une efficacite d'ordre o(log2(n)), tandis 
que l'algorithme classique est d'ordre o(n), avec n le nombre de noeuds. 

Ajoutons maintenant deux noeuds, de valeurs 6 et 11, a l'arbre precedent. On 
obtient l'arbre de la figure 8.4. 



4 




Figure 8.4. Avec deux noeuds en plus 
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Cet arbre n'est plus equilibre, car le noeud de valeur 4 a 3 noeuds dans son SOUS- 
arbre gauche et 5 a droite. Pour le r^quilibrer, il y a plusieurs possibilites, dont 
celle de la figure 8.5. 

6 




4 11 
Figure 8.5. Reequilibrage de l'arbre de la figure 84 

Notons que l'arbre de la figure 8.4 est en fait aussi bon que celui de la figure 
8.5 en ce qui concerne l'efficacite des recherches, car leurs profondeurs sont les 
memes. On peut definir la profondeur d'un arbre comme etant la longueur du plus 
long chemin menant de la racine a une feuille. L'efficacite des recherches se mesure 
en calculant le nombre de tests necessaires pour atteindre chaque objet de l'arbre, ou 
pour en placer un nouveau. C'est dans ce sens que les arbres des figures 8.4 et 8.5 
sont d'efficacit£ 6gale (voir paragraphe ci-dessous sur les B-arbres). 

8.2.1. Algorithmes de manipulation d'arbres equilibres 

La creation d'un arbre equilibre a partir d'un ensemble d'objets ordonnes n'est 
pas difficile. Commencons avec un vecteur v qui contient les valeurs, dans l'ordre 
numerique. On peut considerer le programme 8.2. 
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DEBUT DONNEES n: entier; 
v: TABLEAU [1 ..n] DE entier; 
VAR val, sg, sd: TABLEAU [1 ..n] DE entier; 
pi: entier INIT1; 
racine: entier; 
PROCEDURE arb(i, j: entier) RETOURNER entier; 
VAR x, k: entier; 
SI i=j 

ALORS val[pl]:=v[i]; sg[pl]:=0; sd[pl]:=0; pl:=pl+1 ; 

RETOURNER (pl-1) 
SINON SI j>i 
ALORS k:=(i+j)/2; x:=pl; pl:=pl+1 ; val[x]:=v[k]; 
sg[x]:=arb(i,k-1); sd[x]:=arb(k+1 J); 
RETOURNER x 
SINON RETOURNER 0 
FINSI 
FINSI 
FINPROC; 
racine:=arb(1 ,n) 
FIN 

Programme 8.2. Arbre equilibre 

La proc6dure arb procede par dichotomie, en retournant l'index de la racine de 
l'arbre qu'elle cree. Elle cr6e un arbre equilibre avec les valeurs v[i]..v[j]. Si i=j, 
l'unique valeur est inseree dans un nouveau noeud, qui est en fait une feuille. Si i<j, 
l'index median k est calcule. Un noeud est cree pour v[k]. Les successeurs de ce noeud 
sont calcul& par des appels de arb avec les valeurs a gauche de k, puis les valeurs a 
droite. On retourne l'index du noeud cr66. Enfin, si i>j, aucun element n'est a inserer, 
done la proeddure retourne la valeur 0 (successeur vide). 

Le probleme se corse quand on dispose d'un arbre equilibre' et que Ton veut y 
insurer un nouvel objet. L'operation est intrinsequement couteuse, les algorithmes 
directs 6tant particulierement lourds. 

8.3. Les B-arbres 

La minimisation du nombre de noeuds a examiner pendant le parcours d'un 
arbre est particulierement interessante pour la recherche d'un fichier a travers le 
catalogue tenu par le systeme d' exploitation, ou pour les acces a des objets dans des 
bases de donnees. En effet, dans ces cas, chaque examen d'un noeud peut necessiter 
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une lecture sur support externe (en general un disque). La multiplication des entrees- 
sorties est un facteur determinant pour l'efficacite du systeme global, SUrtOUt dans un 
contexte conversationnel, ou les temps de reponse ne doivent pas depasser certaines 
limites. Ce paragraphe reprend des notions exposees dans [Miranda 19871, qui traite 
de tous les aspects des bases de donnees. 

Les B-arbres [Bayer 19711, [Bayer 19721 apportent une solution efficace ail 

probleme, sans pour autant pretendre a une efficacite parfaite. lis representent un 

point d'equilibre, evitant les algorithmes couteux de restructuration que necessitent 
les arbres binaires equilibres. 

Un B-arbre n'est jamais binaire. Un noeud possede au plus 2*d+l successeurs 
et 2*d etiquettes, oil d est la densite du B-arbre. Ce sont des arbres ordonnes et 
equilibres. Dans le vocabulaire des B-arbres, le terme "6quilibr6" se definit de la 
maniere suivante : la longueur d'un chemin de la racine jusqu'a une feuille est la 
meme, quelque soit la feuille. Nous presentons ici les B-arbres les plus simples, 
avec d=L e'est-a-dire que chaque noeud comporte au plus deux etiquettes (les valeurs 
du paragraphe precedent, appelees cles dans le langage des bases de donnees) et trois 
successeurs. La figure 8.6 montre un tel B-arbre, complet a deux niveaux (8 
6tiquettes pour 4 noeuds). 




Figure 8.6. Un B-arbre complet avec d=l 



Par ailleurs, on applique la regie qu'un B-arbre doit etre rempli a au moins 
50%, ce qui revient a dire, pour d=l, que chaque nceud comporte au moins une 
etiquette. Un nceud qui n'est pas line feuille a au moins deux successeurs. Pour 
comprendre la facon de construire un B-arbre, ou de rechercher un objet, prenons 
l'exemple de l'arrivee d'objets dans l'ordre le plus difficile, soit des objets comportant 
des etiquettes 1, 2, 3, . . .Les etiquettes 1 et 2 se logent dans la racine (figure 8.7). 

Figure 8.7. Apres lecture des 1 et 2 

A la lecture du trois, il faut introduire un nouveau niveau, car la racine est 
pleine (figure 8.8). 
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Figure 8.8. Apr£s lecture du 3 



Le 4 peut s'ajouter dans la feuille qui contient deja le 3. Mais l'arrivee du 
5 provoque h nouveau un bouleversement, car on voudrait l'ajouter dans ce meme 
noeud. Comme il n'y a pas la place pour une troisieme etiquette, le 5 provoque 
"l'eclatement" du noeud comportant le 3 et le 4. L'eclatement consiste en la 
"remontee" de l'etiquette intermediaire (dans ce cas le 4), en creant deux nouveaux 
noeuds pour les etiquettes decentralisees (ici 3 et 5). Le resultat se trouve en figure 
8.9. 




[ MM J 




1/1 i l/WI 1/1 * [AAA \ A" V\AA 

Figure 8.9. Apres lecture du 5 



On voit que le 4 peut se loger au niveau au-dessus, car le noeud comportant 
l'etiquette 2 n'est pas plein. A son tour, le six se logera dans la case qui contient deja 
le 5, car il reste un espace. L'anivee du 7 va encore provoquer des eclatements (figure 
8.10). 



L V I 1/1/1 



L h 1 J Ad 
MiWM \ 1/1 s [AAA 



Figure 8.10. Apris lecture du 7 
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Le 7 a fait eclater le noeud qui contenait 5 et 6. II n'y avait pas de place pour le 
6 dans le noeud au-dessus, qui contenait deja 2 et 4. A son tour, ce noeud eclate, 
provoquant l'arrivee dune nouvelle racine. Cela augmente la profondeur de l'arbre. 

Notons que la profondeur augmente toujours par l'adjonction d'une nouvelle 
racine. Cette propri6t6 assure la continuit6 de la definition de B-arbre, oil chaque 
feuille est a la meme distance de la racine. Les noeuds sont egalement remplis & 50% 
dans le plus mauvais des cas, car chaque noeud cree recoit toujours une etiquette. La 
figure 8.10 montre les limites de l'algorithme simple, car l'arbre comporte trois 
niveaux pour un ensemble d'etiquettes qui n'en necessite que deux. Cela resulte du 
choix de donnees de depart, ou nous avons considere l'ordre le plus mauvais. Le 
resultat est neanmoins tout & fait acceptable, surtout en sachant qu'un ordre aleatoire 
d'anivee des etiquettes mene a un bien meilleur taux de remplissage des nceuds, 

Dans les bases de donnees, ou dans le stockage de fichiers en general, les 
Ctiquettes sont des cles dans le sens de l'adressage dispersee (voir chapitre 2). Sauf 
accident, ceci assure une distribution correcte de valeurs. Si l'efficacite etait vraiment 
critique, on aurait pu eviter l'introduction d'un nouveau noeud dans la figure 8.10. 
Rappelons que l'etiquette 7 fait eclater la case qui contenait auparavant le 5 et le 6. 
La montee du 6 fait eclater la racine, qui comportait le 2 et le 4. La remontee du 4 
cree une nouvelle racine. En fait, il y a une place pour le 4 a cote du 3, comme dans 
la figure 8.11. 



I/1 1 \A/\/\ \A*\A* l/l 




Figure 8.11. Reorganisation pour eviter d'introduire un niveau 



II a egalement fallu reorganiser le trio 5, 6, 7. On note que retrouver une 
organisation optimale est loin d'etre trivial. En general, on accepte l'algorithme de 
base, qui est "raisonnable". Des evaluations des couts des differentes operations se 
trouvent en [Bayer 1972], [Knuth 19731, [Alio 19741. 

Le programme 8.3 recherche une valeur dans un B-arbre, en l'introduisant dans 
l'arbre si elle n'y est pas encore. 
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DONNEES sg, sc, sd, valg, vald, pred: TAB [1 ..taille] DE entier; 

V, racine: entier; 
DEBUTVAR n. pos. orph, nn: entier; 
fini, trouve: bool; 
n:=racine; trouv6:= v=valg[n] OU v=vald[n); 
TANTQUE sg[n]>0 ET NON trouve 
FAIRE SI v<valg[n] 
ALORS n:=sg[n] 

SINON SI v>vald[n] ET vald[n]>0 

ALORS n:=sd[n] 
SINON n:=sc[n] 
FINSI 

FINSI; 

trouve:= v=valg[n] OU v=vald[n] 

FAIT; 

SI NON trouve 
ALORS orph:=0; fini:=FAUX; 
1ANTQUE NON fini 
FAIRE SI vald[n]=0 

ALORS fini>VRAI; 
SI n=0 

ALORS nn:=nouveau_noeud; sg[nn]:=racine;sc[nn]:=orph; 

sd[nn]:=0; pred[nn]:=0; valg[nn]:=v; vald[nn]:=0; 

pred[racine];=nn; pred[orph]:=nn; racine:=nn 
SINON SI v>valg[n] 

ALORS vald[n]:=v; sd[n]:=orph 

SINON vald[n]:=valg(nj; valg[n]:=v; sd[n]:=sc[n]; 
sc|n]:=orph 

FINSI 

FINSI 

SINON nn:=nouveau nggud; 
SI v<valg[n] 

ALORS 6change(v,valg[n]); pos:=1 ;valg[nn]:=vald[n] 

SINON SI v>vald[n] ET vald[nl>0 

ALORS valg[nn]:=v; v:=vald[n]; pos:=2 
SINON valg[nn]:=vald[n]; pos:=3 
FINSI 

FINSI; 

vald[n]:=0; pred[nn]:=pred[n]; 

CAS pos 

DANS 1: sc[nn]:=sd[n];sg[nn]:=sc[n];sc[n]:=orph, 

2: sg[nnj:=sd[n]; sc[nn]:=orph, 

3: sg[nn]:=orph; sc[nn]:=sd[n] 
FINCAS; 

sd[n]:=0; sd[nn] :=0 ; n :=pred[n] ; o rph : = n n 

FINSI 

FAIT 

FINSI 

RN 



Programme 8.3. B-arbres 
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La complexity de cet algorithme montre le travail a faire pour mettre en route 
ce genre de technique. Le programme qui restructure un B-arbre pour lui donner une 
representation optimale est encore plus difficile. Un traitement complet de ces 
problemes depasse la portee de cet ouvrage. 

8.4. Exercices 

1. Ecrire un programme qui decide si un arbre binaire donne est ordonne. 

2. Ecrire un programme qui decide si un arbre binaire donne est equilibre. 
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CHAPITRE 2 

Question 1 

Voir chapitre 7. 

Question 2 

Supposons que nous arrivons dans la situation suivant : 
objet=15, bas=centre=7, haut=8, t[7]=10, t[8]=20 

Alors, l'objet n'est pas trouve et t[centre]<objet. La prochaine valeur de centre 
sera entier((7+8)/2), c'est-Mire 7. L'algorithme boucle. 

CHAPITRE 3 

Question 1 

DONNEES n: entier; 

t: TABLEAU [0..n] DE entier; 
PRECOND n>0; 
DEBUT VAR i, j, temp: entier; 

MONDE i: t[l ..i] est trie; 

i:=1 ; t[0]:= ■ maxint 

TANTQUE i<n 

FAIRE i:=i+l ; 

MONDE j est la position du trou, 
temp:=t[ij; j:=i; 

TANTQUE j>1 ET t[j-1]>temp 
FAIRE t[j]:=t[j-1]; j:=H 
FAIT; 

t[j]:=temp 

FAIT 

FIN 

POSTCOND t[ 1 ..n] est trie 



Algorithmique et programmation 



Le probleme vient du fait que quand j=l, t[j-l] n'existe pas. On cree un t[0], 
fictif, en l'initialisant au plus petit entier negatif. En fait, cette initialisation n'est 
pas necessaire, car si t[0] est r6ftsrenc6, il n'est plus vrai que j> 1. 

CHAPITRE 4 

Question 1 

PROC trouver(nom: chaine(8), n: entier) -> (bool, entier): 
VAR b: bool, x: entier; 
SI n=0 

ALORS RETOURNER (FAUX, 0) 
SINON SI val[n]=nom 

ALORS RETOURNER (VRAI, n) 

SINON (b, x):= trouver(nom, sg[n]); 
SI b 

ALORS RETOURNER (b, x) 

sinon retourner trouver(nom, sd[n]) 

FINSI 

FINSI 

FINSI 

FINPROC 

PROC parent(p, enf: chaine(8), n: entier) -> (bool, entier): 
VAR b: bool, x: entier; 
SI n=0 

ALORS RETOURNER (FAUX, 0) 

FINSI; 

SI val[n]=p 

ALORS SI val[sg[n]]=enf OU val[sd[n]]=enf 
ALORS RETOURNER (VRAI, n) 
FINSI 

FINSI; 

(b, x):= trouver(p, enf, sg[n]); 
SI b 

ALORS RETOURNER (b,x) 

SINON RETOURNER parent(p, enf, sd[n]) 

FINSI 

FINPROC 
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Algorithmique et programmation 

PROC cousin_germain(a, b: chaine(8), n: entier) -> (bool, entier, entier): 
SI n=0 

ALORS RETOURNER (FAUX, 0) 

SINON SI ((val[sg[sg[n]]=a OU val[sd[sg[n]]=a) 

ET (val[sg[sd[n]]=b OU val[sd[sd[n]]=b)) 
OU ((val[sg[sg[n]]=b OU val[sd[sg[n]]=b) 
ET (val[sg[sd[n]]=a OU val[sd[sd[n]]=a)) 
ALORS RETOURNER VRAI 
FINSI; 

SI cousin_germain(a, b, sg[n]) 

ALORS RETOURNER VRAI 

SINON RETOURNER cousin_germain(a, b, sd[n]) 

FINSI 

FINSI 

FINPROC 

PROC aieuKa, enf, n) -> bool: 
VAR b: bool, x, y: entier; 
(b, x):=trouver(a, n); 
SI b 

ALORS (b, y):=trouver(enf, x); RETOURNER b 
FINSI 
FINPROC 

REMARQUE. — Cette version suppose qu'un nom ne peut figurer qu'une fois dans 
l'arbre. Elle n'assure pas l'impression, qui pose le probleme particulier de l'ordre 
d'impression demande. Supposons une procedure trouver-bis, qui comporte un ordre 
d' impression. Alors les resultats seraient imprimes avec les noms des enfants avant 
ceux de leurs parents. Le plus simple est d'empiler le nom h l'interieur de la 
procedure (dans une pile globale) et d'imprimer la pile a la fin du programme. 



Question 3 

Un tel arbre comporte (2 n -l) noeuds et 2 n_1 feuilles. 



CHAPITRE 5 

Question 1 

Le nombre de deplacements est 2 M ■ 1 (voir §5.1.1) 
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Algortthmique et programmation 

2 s4 - 1 secondes = 18446744073709551615 secondes - 

= 307445734561825860 minutes 15 secondes 
= 5124095576030431 heures 15 secondes 
s 213087315667934 jours 15 heures +... 
= 5 838 008 664 843 annees 239 jours +... 

(avec l'annee a 365 jours, c'est-Mire sans prendre en compte les annees bissextiles). 

En arrondissant a 5 800 milliards d'annees, on voit que la fin du monde n'est 
pas pour demain. Le resultat est un bel exemple de l'explosion combinatoire. 

Question 3 

L'ordre est inverse a chaque changement de disque, c'est-a-dire que dj fait des 
cercles dans le sens oppose a celui de dj et, plus gdneralement, dj tourne dans le 
oppose a dj.j . 



Question 4 

DEBUT DONNEES i: entier; 
VAR n: entier; 
n:=1; 

TANTQUE 2 DIVISE i FAIRE i:=i/2; n:=n+1 FAIT 

FIN 

i est Findex du coup, n l'index du disque qui bouge a coup i. Tous les coups 
impairs sont du disque dj. Pour les impairs, on divise i par deux, en enlevant le 

disque 1 (passage au monde des 2) . . . Voir le calcul du pivot dans le probleme des 
permutations. 



Question 5 

Le pivot marche toujouis de droite a gauche et on genere le i-ieme permutation 
de n objets 
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Algortthmique et programmation 

perm: TABLEAU (1 :n) DE integer INIT 1; 
nn:=n; 

TANTQU E nn> 1 

FAIRE (p,q):=i DIVPAR nn; 

SI q=0 ALORS p:=p-1 ; q:=nn FINSI; 

j:=n+1; 

TANTQUE q>1 

FAIRE j:=j-1; 

TANTQUE perm(j)>1 FAIRE j:=j-1 FAIT; 
q:=q-1 

FAIT, 

perm(j):=nn; nn:=nn-1; i:=p+1 

FAIT 



Question 6 

Ordre alphabetique, toujous pour tirer le i-i&me permutation. C'est la meme 
chose en placant les objets a l'envers 

TANTQUE n>0 

FAIRE n:=n-1 ; (p,q):=i DIVPAR fad(n); 

SI q=0 ALORS p:=p-1; q:=fact(n) FINSI; 
j:=1; 

TANTQUE perm(j)>0 FAIRE j:=j+l FAIT; 

TANTQUE j<p 

FAIRE j:=j+1 ; 

TANTQUE perm(j)>0 FAIRE j:=j+l FAIT 

FAIT; 

perm(j):=p+1 ; i:=q 

FAIT 



Question 8 

Oui. Chacun des fact(n) permutations est differente des autres (la relation avec 
la permutation d'origine est bi-univoque). Comme elles sont toutes differentes et il y 
a le compte . . . 
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Algorithmique et programmation 



CHAPITRE 6 
Question 1 

Version avec heuristique : 

PROCEDURE prochain(dir) RETOURNER entier; 
RETOURNER(9 dir=0 

ALORS dirpref(x,y,xf,yf) 
SINON SI dir=4 

ALORS 1 

sinon dir+1 
FINSI 

FINSI) 

FINPROC; 

PROCEDURE dirpref(x,y,xf,yf) RETOURNER entier; 
RETOURNER(SI xcxf ET y<=yf 
ALORS 1 

SINON SI y<yf ET x>=xf 
ALORS 2 

SINON SI x>xf ALORS 3 SINON 4 FINSI 
FINSI 
FINSI) 

FINPROC; 

x:=xs; y:=ys; marque(x,y):=true; dir:=0; empiler(O); 
TANTQUE xoxf OU yoyf 
FAIRE SI prochain(dir)<>dirpref(x,y,xf I yf) 
ALORS dir:=prochain(dir); 

xn:= CAS dir DANS (x,x+1 ,x,x-l) FCAS; 

yn:= CAS dir DANS (y+1 ,y,y-1 ,y) FCAS; 

SI NON barriere(x,y,dir)ETPUIS NON marque(xn,yn) 

ALORS x:=xn; y:=yn;marque(x,y):=vrai; 
empiler(dir); dir:=0 

FINSI 

SINON desempiler(dir); 
SI dir=0 

ALORS -- il n'y a pas de chemin 
FINSI; 

X:= CAS dir DANS (x,x-l ,x,x+1)FINCAS; 
y:= CAS dir DANS (y-l ,y,y+1 ,y) FINCAS 

FINSI 

FAIT 
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Algorithmique et programmation 



Question 8 

DEBUT tab: TAB [1 ..5, 1..5] DE entier INIT 0; 
occupe: TAB [0..16] DE bool INIT FAUX; 
i: entier INIT 1; j: entier INIT 0; 
p, q, val: entier; marche: bool; 
TANTQUE i<6 
FAIRE j:=j+l ; 

TANTQUE occupeO] ET j<16 FAIRE j:=j+l FAIT; 

SI j<16 

ALORS tab[5,i]:=j; occupe[j]:=VRAI; p:=5; q:=i; marche:=VRAI; 
TANTQUE p>1 ET q>1 ET marche 
FAI RE val:=abs(tab[p,q]-tab[p,q-1 ]) ; 
SI occupe[val] 
ALORS marche:=FAUX; 
TANTQUE pc6 

FAIRE occupe[tab[p,q]]:=FAUX; tab[p,q]:=0; 
p:=p+1 ; q:=q+l 

FAIT 

SINON p:=p-1 ; q:=q-1 ; tab[p,q]:=val; OCCUpe[val]:=VRAI 
FINSI 

FAIT; 

SI marche ALORS i:=i+l ; j:=0 FINSI 
SINON i:=i-1 ; x:=5; y:=i; j:=tab[5,i]; 
TANTQUE y>0 

FAIRE occupe[tab[x,y]]:=FAUX; tab[x,y]:=0; x:=x-l; y:=y-1 
FAIT 

FINSI 

FAIT 

FIN 



CHAPITRE 7 

Question 1 

puiss(a,b): SI b=0 

ALORS 1 

SINON SI impair(b) 

ALORS a * puiss(a, b-l) 
SINON puiss(a*a, b/2) 
FINSI 

FINSI 
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Algortthmique et programmation 

puiss(a,b,res): res:=1 ; 

TANTQUE b>0 
FAIRE SI impair(b) 

ALORS res:=res*a;b:=b-1 

FINSI; 

a:=a*a; b:=b/2 

FAIT 



Question 2 
div2(n) = 1/fact(n) 

divfact(n) = n*(n-2)*(n-4r.„ 

(n-l) * (n-3) * . . . 

= SI pair(n) 

alors 2 n * (fact(n/2)) 2 / fact(n) 

SINON fact(n) / 2 n " 1 / (fact((n-1 )/2)) 2 
FINSI 

DEBUT DONNEES n: entier; 
VAR div2: entier INIT1 ; 

TANTQUE n>1 FAIRE div2:=div2/n; n:=n-1 FAIT 

FIN 

DEBUT DONNEES n: entier; 

VAR haut, bas, divfact: entier; 
haut:=1 ; bas:=1 ; 
TANTQUE n>1 

FAIRE haut:=haut*n; bas:=bas*(n-1); n:=n-2 
FAIT; 

divfact:=haut/bas 

FIN 
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Algortthmique et programmation 



CHAPITRE 8 
Question 1 

DEBUT DONNEES taille, racine: entier; 

sg, sd, val: TAB (1 ..taille) DE entier; 
PROCEDURE ord(n, min, max: entier) RETOURNE bool; 
Sln=0 

ALORS RETOURNER VRAI 
SINON SI val(n)>max OU val(n)<min 

ALORS RETOURNER FAUX 

SINON RETOURNER(ord(sg(n), min, val(n)-1) 
ET ord(sd(n), val(n)+1 , max)) 

FINSI 

FINSI 

FINPROC; 

imprimer(ord(racine, • maxint, maxint)) 



Question 2 

DEBUT DONNEES sg, sd, val: TABLEAU (1 ..taille) DE entier; 
racine: entier; 
VAR x: entier; 

PROC equ(n, poids: entier) RETOURNE bool; 
VAR eqg, eqd: bool; 

pg, pd: entier; 
SI n=0 

ALORS pOlds:=0; RETOURNER VRAI 
I SINON 6qg:=equ(sg(n), pg ) ; §qd:=§qu(sd(n), p d ) ; 

poids:=pg+pd+l ; 

RETOURNER (eqg ET eqd ET abs(pg-pd)<2) 

FINSI 

FINPROC; 

imprimer(equ(racine, x)) 
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k\gorithmique\ et programmation est Issu 
d'un enseignement delivre depuis vingt-cinq ans 
U des etudiants de 1 e 1 et 2 e l cycles universitaires 
et dans le cadre de la formation permanente. 
II s'adresse done aussi bien aux etudiants 
qui desirent se specialiser en informatique 
qu'auxj professionnels de I'entreprise qui 
souhaiteraient acquerir des elements de rigueur 
en programmation. 

A partir des exemples ciassiques (par exemple 
les tours d'Hanol), la description et I'examen 
des differentes methodes de construction 
d'algorithmes sont proposes. Des exercices avec 
leur solution concluent chaque chapitre et 
resumentj I'essentiel des points abordes. 
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